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ABSTRACT 

Fundamental physical and chemical processes, and measurement 
techniques on the D region are reviewed. Design considerations about 
a partial-reflection system are made, and the main characteristics of 
the partial -ref lection system at the University of Illinois are pre- 
sented. The nature of the partial reflections are discussed, parti- 
cularly reflections produced by gradients in electron density and by 
random fluctuations in a locally homogeneous random medium. Possible 
reasons for disagreement between partial reflections and rocket measure- 
ments are discussed. Some suggestions are made to improve partial- 
reflection data reduction, including the use of only maximums of the 
reflections and deconvolution of the data. The results of partial- 
reflection measurements at Wallops Island, Virginia during the 1971- 
1972 winter are presented and compared to rocket measurements. 
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1. THE IONOSPHERIC D REGION 
1 . 1 General Characteristics 

The ionospheric region between SO and 90 km of altitude is known as 
the D region. Although the lowest part of the ionosphere, it is the least 
well understood, due to the complexity of the processes in effect and to 
the difficulty in making measurements. 

The D region is characterized as a region of weakly ionized plasma 
with complex ionic structure. It is the only atmospheric region where 
both positive and negative ions are present in significant concentrations. 
Above 90 km, negative ions are quickly detached by ultraviolet solar 
radiation. 

8 9 3 

The electron densities are very small, of the order of 10 to 10 m 
during the day, but because of collisions with the relatively dense 
neutral atmosphere, a strong attenuation is produced in HF electromagnetic 
waves propagating in this region. 

Being primarily governed by solar radiation, the ionization of the D 
region almost disappears at night. The seasonal behavior of the 
ionization constitutes an exception to the solar control, and during some 
winter days very high electron densities are observed. Such abnormal 
behavior is part of the so-called winter anomaly, and will be discussed 
later. 

The temperature structure shows a negative gradient with height, in 
an altitude region of 50 to 80 or 85 km, which is known as the mesosphere. 
The altitude of 85 km corresponds to the coldest level in the entire atmo- 
sphere (about 180°K in midsummer), and is known as the mesopause. Above 
the mesopause the gradient of temperature is positive, as a consequence of 
downward heat conduction from the thermosphere. During the winter the 
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temperature profile of the mesosphere shows strong fluctuations, and the 
minimum near 85 km is not clearly identified. 

The mass of air that constitutes the D region is in turbulent motion, 
and the relative composition of the: major neutral constituents is the same 
as that at ground level. Minor and ionized constituents do not follow 
this rule, due to the presence of sources and sinks. 

In the remainder of this chapter, the principal aspects of the 
chemistry and dynamics of the D region, and the techniques of measurement 
of the more important parameters and constituents will be discussed. 

1 . 2 Scnwoes of Ionization 

The D region is formed by the most penetrating radiations. The fol- 
lowing sources are considered as important on its formation, and will be 
discussed below: 

1) Solar Lyman-a radiation 

2) Solar X-rays of wavelengths less than 10 A 

3) Galactic cosmic rays 

1.2.1 Solar Lyman-a.. The strong solar radiation of Lyman-a (1216 A) 
can penetrate the D region, due to the small absorption cross section of 
the atmosphere above 100 km of altitude, for this wavelength. 

Lyman-a radiation can ionize only constituents with ionization 
potential less than 10 eV, in particular nitric oxide (ionization 
potential of 9.25 eV) . Nicolet and Aikin [1960] estimated that the 
ionization of NO by Lyman-a is the major source of ionization in the 
altitude range of 65 to 85 km. Recent rocket measurements of NO 
concentration [ Meira ., 1971] confirmed the importance of Lyman-a in the 
formation of the D region. Calculated production functions of NO, using 
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the NO concentrations as measured by Meira [1971], are shown in Figure 1.1 
[Aikiriy 1972]. The role of Lyman-a as the principal ionizing source 
between 65 to 85 km during quiescent solar conditions is confirmed in this 
calculation. 

1.2.2 Solar X-rays of wavelength less than 10 A. All atmospheric 
constituents can be ionized by X-rays. X-ray photons with wavelengths 
greater than 10 A are absorbed above 90 km, but photons with wavelengths 
less than 10 A are able to penetrate the D region, that, in this way, can 
be considered as a tail of the E region. 

The intensity of the solar X-ray radiation is extremely variable. 
Under quiet sun conditions, the intensity through the region of wave- 
lengths less than 100 A remains practically constant over a period of 24 
hours. Over longer periods of time, however, considerable variations are 
observed for wavelengths less than 10 A, and the intensity may change by a 
factor of 10 to 100 over a period of a few days {Mandel 's tarns 1965]. 

During solar flares, a change in solar flux and spectral composition can 
be observed. The spectrum shifts towards shorter wavelengths, and sharp 
increases in the flux for A <. 5 A are frequently observed. During a class 

I flare, the radiation intensity in the wavelength range of 1 to 8 A can 

2 1 

reach values of the order of 10 erg cm s . 

In Figure 1.1 a calculation of the ion-pair production function is 
shown for the solar fluxes and solar zenith angles listed in Table 1.1 
[. Aikin 3 1972]. The curve X-ray I corresponds to quiescent conditions at 
low solar activity. During this time the contribution of X-rays to the 
ionization of the D region is negligible. Curve X-ray II corresponds to 
conditions of an X-ray enhancement event. Curve X-ray III shows the 



ALTITUDE (km) 



q (cm -3 sec -1 ) 

Figure 1.1 Ion-electron production functions produced by Lyman a (^q) . X-rays 

(X-ray I, X-ray II and X-ray III), cosmic rays and ultraviolet radia- 
tion of wavelength less than 1118 & (qQ .) . (F rom Aikin [1972]). 



TABLE 1.1 


X-ray fluxes and solar zenith angles corresponding 
to the electron-density profiles of Figure 1.1. 
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effect of X-rays during a class I flare. On such occasions X-rays are the 
dominant source of ionization down to 65 km. 

1.2.3 Galactic cosmic rays. Galactic cosmic radiation (GCR) con- 
stitutes the most important source of ionization below 65 km [Nicolet and 
Aikin , 1960; Webber^ 1962]. Calculated ion-electron production functions 
produced by galactic cosmic rays are shown in Figure 1.1 for maximum and 
minimum of solar activity, and a latitude of 50°. 

1.2.4 Solar ultraviolet radiation of wavelengths less than 1118 A. 
Wavelengths between 1216 ft and 1026.5 ft penetrate the D region in small 
intensities, so that the ionization of minor constituents by such radia- 
tion does not constitute an important source of ionization. The major 
constituents in ground state cannot be ionized by such radiation. If 
excited metastable species are present in sufficient concentrations, 
however, their ionization can contribute appreciably to the formation of 


the D region. 

Ionization of 0 ( 1 A ) by such mechanism has been suggested as a 
^ Q 

source of ionization of the D region. CL(*/\ ) has an ionization threshold 

4 d 

at a wavelength of 118 ft, and its ionization cross section is estimated as 
-18 2 

3 x 10 cm . First calculations of the ion-electron production function 
resulted in a source as important as the ionization of NO by Lyman-a. In 


such calculations, however, atmospheric absorption by CO^ was not con- 
sidered. Ionization rates including C0 9 absorption [Huffman et al. t 
1971] are of one order of magnitude smaller than the ionization produced 


by Lyman-a. A calculation of the ionization production function of 
O^C^A ), including the absorption by C0 9 is shown in Figure 1.1. 
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1 . 3 D -Region Chemistry 

The gases that constitute the atmosphere are chemically inert at low 
altitudes. In the ionosphere, however, they are excited and dissociated 
by solar radiation, and become extremely active. A multitude of reactions 
take place. The study of such reactions and the resulting distribution of 
neutral and ionized species is the objective of ionospheric chemistry. 

As it is impossible, at the present state of knowledge to obtain a 
global model of the D region, involving all chemical reactions, normally 
particular models are developed, appropriate to a group of constituents, 
in a certain range of altitudes. It is usual to divide the £>-region 
chemistry into positive-ion chemistry, negative-ion chemistry and neutral 
chemistry. Some aspects of each will be presented below. Electron loss 
processes encompasses both positive- and negative-ion chemistry, and will 
be presented in a separate section. 

1.3.1 Positive-ion chemistry , The principal aspects and main 
problems related to the positive-ion chemistry of the D region have been 
presented in review papers by Donahue [1972] and Thomas [1974]. 

Mass spectrometric measurements of positive-ion composition of the D 

region have been successfully made since 1963 [Narcisi and Bailey t 1965]. 

The results obtained by Narcisi and Bailey are shown in Figure 1.2. 

+ + + + 

Above 82 km, 32 , (0 ? ), and 30 (NO ) are the dominant ions. Below 
82 km, water cluster ions of mass 37 + (H + • (H^O) and 19 + (H + *(H 7 0)) are 
the most predominant. It is possible that heavier hydrated ions of the 
form H + * (I^O) are present, but they are dissociated by collisions with 
the rocket, and are not observed in the measurements. 28 + is 

observed in small quantities, since the ions formed by X-ray ionization of 



ION DENSITY — IONS/CC 


Figure 1.2 Positive-ion composition of the D region as measured 

by Na.rci.si and Bailey [196S] . 






9 


. + + + 

N 7 disappear very rapidly by the reaction + ^ + ^2‘ * n 

turn dissociatively recombines with electrons or with negative ions. 

The origin of the hydrated ions is one of the leading problems in 
l?-region chemistry. Fehsenfetd and Ferguson [1969] suggested a scheme of 
reactions which is based on the formation of the cluster N0 + *H o 0 in a 
three-body reaction. This scheme was later supplemented by another 
starting with 0^ + instead of N0 + [Ferguson and Fehsenfetd , 1969]. The 
complete reaction scheme is shown in Figure 1.3 [Donahue, 1972], 

The reaction 0 ? + + 0^ + M -»■ 0^ + + M, the initial step in the proposed 
formation of the hydrated ions is a fast reaction in the lower ionosphere, 
and provides an effective way of changing 0^ ions in hydrated ions. 

Above 80 km in daytime the reaction 0^ + + 0 -+ 0^ + + 0^ is very effective 
in reducing 0^ + ions, and can explain the disappearance of hydrated ions 
near 82 km. 

Recent calculations of the ionization rate of O-^A ) have shown that 

2 g 

this ionization is much smaller below 80 km than believed before. As a 
result, the 0^ + source will not be enough to explain the measured ion 
concentrations below 77 km. 

The reaction scheme starting with N0 + also presents some problems 
[Donahue, 1972]. In the first place, the first created ion is H^O* • (H^O) ? , 
and not H^O^fH^O), that has been observed to be the dominant ion; second, 
the conversion of N0 + to the hydrates via the first three-body reactions is 
too slow, relative to the dissociative recombination of the cluster ions, 
and finally, if the NO measurements of Meira [1971] are correct, and if 
the ionization rate of NO is of about 10 cm s , as commonly accepted, 

+ X _ *7 

NO would be the dominant ion at 80 km, with density close to 3 x 10 cm . 
This situation does not correspond to the facts. 




Og + H 2 0 + M — * 02- H 2 0 + M 
2- 2.8 x ICT 2 ® cc/sec M = N 2 

Figure 1.3 Reaction scheme for conversion of 0 2 + and N0+ to hydrates. Three- 
body reaction rate constants in units of 10 - 28 cm & sec" 1 . Two- 
body rate constants in units of 10“^ cm-^ sec" 1 (From Donahue [1972]) 
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This problem would be solved if the ionization rate of NO by Lyman-a 

at 80 km were not as high as believed. Donahue [1972] suggests as 

possible explanations for a low ionization rate that the absorption of NO 

+ 

by a Lyman-a photon would not lead directly to NO , but rather to an 
excited state of NO, the ionization rate increasing by collisions at low 
altitudes, or that there would be a variation in the cross section of NO 
within the width of 1 A of the Lyman-a line, the structure being a 
function of pressure. 

1.3.2 Negative- ^on chemistry. Measurements made up to the present 
time of negative-ion composition of the D region have yielded conflicting 
results [Narcisi et al, 3 1972a; Arnold et al. 3 1971]. Consequently the 
knowledge of the negative-ion chemistry is based on laboratory measure- 
ments and theoretical models. A review relating the most recent 
progress has been given by Thomas [1974] . 

1.3.3 Electron loss processes . The predominant electron loss source 
above 80 km is dissociative electron-ion recombination with NO + and . 

At lower altitudes three-body reactions become important, and electron 
attachment to 0^ can be realized with N^ or 0 ? as the third body in the 
reaction. 

The effective recombination coefficient, a can be measured at any 

eft’ 

occasion when a sudden variation in electron production is observed, as 
during solar flares [Montbriand and Belrose 3 1972] and solar eclipses 
[Sechrist 3 1970]. 

According to the results obtained by Sechrist [1970], above 80 km 
lies between 10 ^ and 10 . The rapid variation with height of & e ££ 

in this region corresponds to the gradient of electron density observed 
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near 85 km. Below 75 km, a ff increases sharply from 4 x 10" 5 to 10 3 cm 3 

s" 1 . Between 78 and 85 km, a g£f appears to be nearly constant with a 

5 3-1 

value of about 4 x 10 cm s . The loss mechanisms in this altitude 
range are not well understood, but electron recombination with hydrated 
ions is probably the dominant electron loss process [Beidj 1970] . 

1.3.4 Neutral chemistry . 0-region neutral chemistry is related to 
the chemical processes involving minor neutral constituents. All 
atmospheric constituents in the D region, with exception of 0 and N ? fall 
in the minor category. They include rare gases, metallic atoms such as 
Na, Ca, Al, Mg, Ni, Cr, Fe, and a number of molecules formed from the 
elements nitrogen, oxygen, hydrogen, and carbon. 

To the present time, measurements of minor neutral constituents have 
been made by using particular techniques that are able to measure only one 
type of constituent. For atomic oxygen, for example, no reliable 
technique of measurement has been developed until now. As a result, our 
knowledge of the concentration of minor constituents is very 
unsatisfactory. 

The chemistry of the most important minor netural constituents has 
been described in review papers by Strobel [1972] and Thomas [1974]. 

1.4 Measurement Techniques 

The complexity of the parameters involved, the relatively high 
density of the atmosphere, and the small concentration of some of the most 
important constituents are some of the factors that contribute to the dif- 
ficulty in performing 0-region measurements. Satellites cannot be used, 
and measurements have to be made by using rocket or ground-based 
experiments . 
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Rockets have the disadvantage of giving only instantaneous informa- 
tion of the parameter involved, and are not well suited to time variation 
or synoptic measurements. Most D-region ground-based experiments, as will 
be described below, do not furnish the accuracy that would be desirable. 

Because of the complexity of the processes involved, Z)-region 
measurements must be carefully coordinated, to take maximum advantage of 
the data obtained to the best theoretical interpretation, and in planning 
the next generation of experiments \Bowhill> 1969]. In the next sections 
the following types of measurements will be discussed: neutral structure 

and composition, ion density and composition, and electron density and 
collision frequency. 

1.4.1 Neutral structure and composition. Measurements related to 
the neutral structure of the V region include measurements of density, 
neutral temperature, winds, minor neutral constituents and turbulence. 

- Density measurements : Below the turbopause, measurements of gas 

density give the concentration of each major constituent, since the 
composition is nearly the same as at the ground. 

Density profiles have been obtained by using falling spheres and 
pressure gauge experiments. 

In the falling sphere experiment [Bartman et al. s 1956] a collapsed 
sphere is launched to a high altitude and then ejected and inflated after 
the rocket power has been exhausted and atmospheric drag has reached a 
tolerable magnitude. 

Using falling sphere experiments Faucher et al. [1963, 1967] have 
performed several density measurements between the heights of 20 to 135 
km; they estimated an error of about 4 percent at 110 km, increasing to 50 


percent at 135 km. 
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Pressure gauges aboard rockets were used by La Gow et al . [1959] in 
the determination of air density to an altitude up to 210 km. 

_ wind s and neutral temperature : The most useful method of measuring 

winds and neutral temperature is the rocket-grenade method. The experiment 
is made by carrying aloft and ejecting grenades from rockets. 

The fundamental parameters measured are the positions in space, the 
total travel time of sound waves to the ground, and angles of arrival of 
successive sound waves at an array of microphones located on the ground. 
From these parameters, wind velocities and temperature can be determined. 

The grenade experiment has been used successfully several times, and 
a large amount of data obtained [Stroud et at. t I960]. 

The accuracy of the measurement is limited by errors in measurements 
of the sound arrival times, and the resulting error below 75 km is 
generally less than ± 15° in wind direction. For higher altitudes, the 
intensity of the sound at the ground is so weak that the errors sometimes 
increase by a factor of 10, 

- Minor neutral constituents : No technique has been developed to the 

present time that is able to measure different minor neutral constituents 
simultaneously. Mass spectroscopy has not been able to detect minor con- 
stituents in the D region. The information we have is the result of 
measurements made by rather specialized techniques, usually only appli- 
cable to one particular constituent, or by indirect means. Most measure- 
ments have been made by using absorption spectroscopy and glow emissions. 

Absorption spectroscopy has been used successfully in measurements 

V 

of ozone. 

Diurnal distributions of ozone were measured by Johnson et at. 

[1952], Weeks and Smith [1968], Reed [1968],#ays et at, [1972] and 
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Miller and Ryder [1973]; nocturnal distributions were obtained by Carver 
et al. [1966] . 

Dayglow emissions have been used in measurements of sodium and nitric 
oxide. 

Sodium profiles were obtained by Donahue and Meier [1967] by 
measuring the brightness of the sodium V lines as a function of altitude, 
using rochet-borne photometers. 

Ground-based measurements of sodium dayglow brightness with Zeeman 
photometers have given information about the total concentration and time 
variations of atmospheric sodium [Blamont and Donahue ^ 1964] . 

Barth [1964], Pearce [1969], and Meira [1971] obtained NO concentra- 
tions in rocket experiments, by measuring the dayglow emission of NO in 
the gamma bands, using scanning ultraviolet spectrometers aboard rockets. 

- Turbulence : Turbulence parameters have been measured by release of 

chemicals in the atmosphere [ Zimmerman, 1965] and from reflections of 
radio waves from meteor trails [Greenhou and Neufeld> 1959]. Sodium and 
other chemicals have been released from rockets. The sodium trail is 
observed at twilight due to resonant scattering of 5890 A sunlight. Some 
other chemicals have been used that produce a chemiluminescent trail. 
Turbulence parameters are obtained by taking pictures of the trail and 
analyzing its spreading as a function of time. The spread of the trail is 
due to molecular diffusion, wind shear and turbulence. There are some 
problems in interpreting the results, since the pictures taken contain 
only two-dimensional information, and the passage of the rocket introduces 
a pertubation in the atmosphere, and can alter the turbulent regime, 

1.4.2 Ion composition and density. Measurements of positive-ion 


composition have been successfully performed several times, by using 



rocket-borne mass spectrometers [Naroisi and Bailey > 1965; Naroisi et at. t 
1972a, b; Krankowsky et at., 1972]. Measurements of negative-ion 
composition in the D region, that have been carried out by Naroisi et at. 
[1972a] and Arnold et at. [1971], produced conflicting results; such 
measurements are very difficult to be made, due to the interactions of the 
ions with the structure of the rocket and spectrometers. 

Ion densities have been measured by different kinds of probes. 

Naroisi and Bailey [1965] used a spherical electrostatic analyzer, that 
consists of a spherical collector surrounded by a concentric wire mesh 
grid. When the collector has a negative polarization, and the grid is at 
the rocket potential, the collector current is proportional to the 
positive-ion density. Bale et al. [1968] used parachute-borne blunt 
probes in measurements of ion concentrations. 

1.4.3 Electron density and collision frequency measurements. The 
techniques of measurement of electron density and collision frequency can 
be classified into rocket and ground-based techniques. 

- Rocket techniques : The first rocket measurements of electron 

density of the D region were obtained as early as 1947 [Seddon, 1953], 
utilizing the Doppler effect. Since then different techniques have been 
developed and combined, so that rocket measurements became the most 
effective way of measuring electron density and collision frequency in the 
lower ionosphere. 

Rocket techniques utilize effects introduced by the ionosphere on 
propagation of radio waves, or current probes. The radio propagation 
experiments are based on Doppler effect, Faraday rotation, pulse delay and 
absorption, and on phase variations. Several measurement systems have 
been developed. 
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The University of Illinois, in cooperation with GCA Corporation, 
developed a measurement system incorporating both differential absorption 
and Faraday rotation, and a Langmuir probe. Since 1964 several successful 
measurements have been performed using this system. Bowhill [1965] has 
described the propagation experiment. The Langmuir probe, that will be 
discussed below, adds height resolution to the measurement. 

- Langmuir probe : The Langmuir probe is a direct current instrument, 

that works by collecting charged particles in a plasma, by using a 
collector a few volts positive (in electron-density measurements) or 
negative (in positive-ion measurements). The probe current can give, in 
this way, an indication of electron density and of ion number density. 

Its usefulness at altitudes below 100 km, however, is severely limited by 
the lack of an adequate theory to describe the collection of electrons and 
ions by a space vehicle when the mean free path of the particles is not 
much larger than the Debye length. Such a problem has been solved by 
using Langmuir probes simultaneously with radio propagation rocket 
measurement, as described before. The Langmuir probe in such simultaneous 
measurements has the advantage of increasing the height resolution to an 
order of 10 m. 

- Ground-based experiments : Ground-based electron density and 

collision frequency measurements have been obtained from the following: 

VLF propagation, incoherent scattering, cross modulation and partial 
reflection. 

Only the partial-reflection technique will be discussed here. 
Information about the other techniques can be obtained on review papers 
[Fejer, 1970; Seahris 1974], 



_ Partial reflection: The partial-reflection technique is based on 


the measurement of weak reflections produced by small irregularities in 
the D region. If an HF pulse is sent to the ionosphere, a series of small 
reflections are observed between the heights of 50 and 90 km. If pulses 
with polarization corresponding to the ordinary and extraordinary modes of 
propagation are transmitted, from the ratios of the amplitudes, or from 
the phase of the reflections of the two modes, electron densities can be 
calculated, if a collision frequency profile is assumed. 

Weak reflections produced by irregularities in the ionosphere were 
observed by Dieminger [1952], by examining conventional ionograms. Such 
reflections were used for the first time to calculate electron densities 
by Gardner and PaMsey [1953]. The technique of Gardner and Pawsey is 
based on the measurement of the amplitude of the reflections, and for this 
reason sometimes called a differential absorption technique. Belrose and 
Burke [1964] introduced the generalized magnetoionic theory in the expres- 
sions for calculation of electron densities, and developed the 
experimental techniques. A theory based on the work of Belrose and Burke 
will be presented below. 

If the reflection produced at a height h is due to a small dis- 
continuity in the refractive index, the coefficient of reflection is given 
by: 


P. 

o } x 

where n ^ and are the refractive index immediately above and below the 
irregularity, and the indices o and x refer to the ordinary and extra- 
ordinary modes of propagation. In expression (1.1) it was assumed that 


2 1 
n 2 +n l 


( 1 . 1 ) 
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the ionosphere is horizontally stratified, and that the propagation is 
quasi- longitudinal, so that the wave equation is decoupled for the 
ordinary and extraordinary inodes, and the classical Fresnel expression for 
the coefficient of reflection can be used. 

As the reflection observed are very small, it can be written: 

= n, and = dn, resulting: 


dn 
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( 1 . 2 ) 


The expressions for the refractive indices are given by the 
generalized Appleton-Hartree magnetoionic formulas [Sen andWyller, 1960], 
for quasi- longitudinal approximation: 
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where i<> = 2 ttj, f being the operating frequency; 

\ 

Wj. = 2irf^ cos0,> being the gyrofrequency , and 6 the angle between 
the earth's magnetic field and the path of propagation; 

W o is the plasma frequency; 

C^/2 (a:) , (r) are integrals which have been tabulated 

by Buvke and Hava [1963]; 
a is the velocity of light in free space; and 
v is the collision frequency of monoenergetic electrons. 

Assuming that the irregularities are produced only by fluctuations in 


electron density, and not in collision frequency, it results from (1.2) 
and (1.3). 
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( 1 . 4 ) 


where |f? | and [i? | are the amplitudes of the ordinary and extraordinary 
coefficients of reflection, respectively. 

The intensity of the signal reflected at a height h, when received 
at the ground is given by: 
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and the ratio of the intensities of the extraordinary to ordinary signals 
is given by: 

( 1 . 6 ) 
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If reflections from two different heights h ^ and ^ are measured, it 
can be written: 
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where 
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the indices 1 and 2 referring to the heights h ^ and h ? , respectively. 

Taking the distance h„-h. sufficiently small, the factors K and K 

Z. J. X ( 

can be assumed constants in the interval h ? -h^, resulting: 


A 1 noyv = A ln( W ' 2iK x~ K o )Ah 


( 1 . 8 ) 


where A h = 

The value of K 0 * x is given by: 
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where N is the electron density 0 

From expressions (1.8) and (1.9), it results 
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Expression (1.10) can be used in the determination of the average 
electron density between the heights h . and h*. The quantities |Fj/|fl 
and F q are obtained from equations (1.4) and (1.9), if a collision 
frequency profile is assumed. 
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Some comments are necessary about the signal processing techniques 
and the theory used. 

The received signals are produced from reflections of transmitted 
pulses with width between 25 and 50 y sec. Considerations of bandwidth 
and height resolution limits the lower and upper values of the pulse 
width. The pulse repetition rate normally used ranges from 0.5 to 60 
sec 1 , depending on the processing capabilities of the system. Average 
values of the scattering cross sections at height intervals of approxi- 
mately 1 km must be obtained from such signals. The average scattering 
cross sections are determined in the following ways: 

1) By measuring the amplitudes of the received signals corresponding 
to a given height in every sample, and taking the average power over all 
samples : 
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( 1 . 11 ) 


where M is the number of samples, (a J is the amplitude of the received 

( S J Ob 171 

signal, from a given height, in the m-th sample, and A is the average 

O £C 

value of , Several ways of reducing the influence of noise are used, 
and some of them will be discussed in Chapter 2, 

As the pulse width is finite, reflections from a height h are really 
produced in a height interval of ± Wc'/ 4, W being the pulse width, and <?' 
the velocity of light in the medium. The average value of the signals 
measured at a height h , in this case, will be representative of the 
scattering cross section at this height if the reflections are equally 
distributed over the interval ± Wo 1 / 4, during the sampling period, and if 
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the scattering cross section per unit volume does not change very rapidly 
with height. The consequences of a breakdown of the above assumptions 
will be discussed in Chapter 4. 

2) By measuring only the peaks of the received signals, and taking 
the average of the measured peaks at each height. In this procedure it is 
assumed that the peaks of the signals correspond to reflections at a given 
height, and are not the result of interference of signals reflected from 
different heights. In practice, this method offers the disadvantage that 
normally only a few reflections are observed from certain height 
intervals, principally between 76 and 80 km. As a result long sampling 
and processing times are necessary. The influence in the electron-density 
profiles of using one or the other of the above methods will be discussed 
in Chapter 4. 

The theory of scattering used in the determination of the ratio Iff 1/ 
\R q \ is a second point that deserves some comments. 

The ratio |l? \/\E | was obtained by assuming that the reflections 
cc o 

are produced by irregularities over height intervals much smaller than one 
wavelength, and that the atmosphere is horizontally stratified, at least 
over a region comparable to the first Fresnel zone. If the reflections 
are produced by irregularities distributed continually over the height 
interval, if there are a great number of discrete irregularities in the 
volume occupied by the pulse, or if the atmosphere is not horizontally 
stratified, a new model of reflectors must be introduced, or the use of 
equation Cl<4) must be justified. 

BetvoQe and Burke [1964] showed that reflections produced by irregu- 
larities distributed over a volume are related to reflections produced by 
a small discontinuity by a geometric factor, that is the same for both the 


24 


ordinary and extraordinary modes, and in consequence the ratio 
should be the same for both models of reflectors. This result, however, 
depends on assuming that the collision frequency and correlation function 
of the irregularities remain constant over the reflection volume, and that 
the atmosphere is horizontally stratified. 

The effect of irregularities in the refractive index produced by 
irregularities in collision frequency was considered by Piggott and Thrane 
[1966]. As shown by Gallet [1955], however, if the irregularities are 
caused by adiabatic turbulent mixing on the gradient of electron density, 
the pressure fluctuations, and consequently the collision frequency 
fluctuations, are orders of magnitude smaller than the electron-density 
fluctuations . 

Flood [1968] extended the theory of partial reflections by 
considering a reflector model consisting of irregularities continuously 
distributed over a homogeneous background medium. It was assumed an 
ionosphere horizontally stratified over an area at least equal to the 
first Fresnel zone, and quasi- longitudinal propagation. 

If the above assumptions are valid, the wave equation is decoupled 

for the ordinary and extraordinary modes of propagation and the scattering 

cross section can be obtained by using the same techniques developed for 

scattering from an isotropic medium [Tatarskiy 1961; Booker 3 1959]. The 

2 2 

expression obtained by Flood for the ratio A /A ~ is given by: 
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where y = y q - w+w^/v, x is the pulse width, e is the charge of 

the electron, m is the mass of the electron, and a the velocity of light 
in free space. Figure 1.4 shows electron-density profiles obtained from 
the same A /A profile, by applying the theories of Belrose and Burke 
[1964] and Flood [1968]. The difference between the two profiles is at 
most 100 percent. As pointed out by Holt [1969], if the differential 
absorption inside the scattering volume is not too great. Flood's theory* 
reduces to Belrose and Burke's theory. If the differential absorption 
inside the scattering volume is great and changes with altitude, the two 
theories will give results considerably different. Such a situation 
occurs normally above 80 km. 

Cohen [1971] developed an expression for the scattering cross section 
for a reflector model including the following characteristics: 

a) There is a finite number of irregularities distributed 
in a random way; 

b) The irregularities are uncorrelated; 

c) The ionosphere is horizontally stratified. 

Assuming that the probability of occurrence of M different irregu- 


larities follow a Poisson distribution, Cohen arrived at the following 
expression for A ^ 

^ OjX 



e^BNCoh/2] 
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Figure 1.4 Electron-density profiles obtained using Flood [1968] 
theory (- 0 -- 0 -) and Belrose and Burke [1964] theory 
( ) . (From Flood [1968]). 
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where the variation of electron density of the irregularities is supposed 
to have a height dependence of the form A N(h) - $N(h)g(h) being the 
constant g(K) a geometric factor; y is the average number of irregularities 
per unit length, which was assumed to be a constant independent of height; 
|aj is a geometrical factor; t is the pulse duration, and the other symbols 
are as defined before. 

In Chapter 4 the problem of comparison of different scattering 
theories will be discussed again. 

The real structure of the irregularities is not known, and thus it is 
not possible to define the more appropriate model. The problem of deter- 
mination of the structure of the irregularities will be discussed in 
Chapter 4. 

If the ionosphere is not horizontally stratified, there will be a 
coupling between the ordinary and extraordinary modes. If a pulse with 
one mode of polarization is transmitted, it will be depolarized, and 
reflections containing both modes will result. The wave equation in such 
a case cannot be separated, and its solution is very difficult. An 
approximate solution for the scattering cross section of a magnetoionic 
medium containing continuously distributed irregularities over a 
homogeneous background was developed by Simonieh and Yeh [1971] . In the 
model of Simonieh and Yeh the collision frequency was assumed as zero, and 
so the results cannot be applied directly to partial reflections. The 
wave equation was solved by using an approximate technique developed by 
Lighthill [I960]. The calculated scattering cross sections show that for 
the range of parameters encountered in partial-reflection experiments, the 
results are practically the same of Flood's [1968] theory. 



28 


Coupling between the ordinary and extraordinary modes in a 
horizontally stratified medium where the propagation is not longitudinal 
was discussed by Connolly and Tanenbaum [1972] . 

The results of Connolly and Tgnenbaum [1972] show that for a hori- 

10 -3 

zontally stratified ionosphere, and for an electron density of 10 m > 
the coupling coefficient is always less than 3 percent. This is not large 
enough to be considered in partial -reflection calculations, since the 
error introduced is much smaller than the imprecision of the technique. 

Until now only electron-density determinations based on measurements 
of amplitude of the reflections were discussed. Von Biel et at. [1970] 
developed a method of determining electron densities from the differential 
phase between the ordinary and extraordinary reflections, using the same 
model of reflector and the same approximations of Flood [1968] , and 
utilizing measurements of amplitude of the reflections at different 
directions to determine indirectly the phase. The electron-density 
profiles obtained are very similar to the profiles measured by Flood’s 
differential absorption technique, as should be expected, since both are 
based on the same theory. 

Austin [1971] developed a method of determination of electron 
densities by measuring directly the phase difference between the ordinary 
and extraordinary reflections. 

1.5 Comparison of Electron- Density Profiles Obtained by Different 

Techniques . 

A comparison of electron-density profiles obtained by different 
techniques has been given by Seohrist [1974]. Differences between partial 
reflection and rocket measurements are of particular importance for this 
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work. Such differences occur principally for heights greater than 75 km, 
where partial-reflection profiles show a tendency to produce a valley in 
electron density that is not observed in rocket profiles. Above 80 km 
partial-reflection profiles are normally too low if compared to rocket 
profiles. Such differences will be discussed in Chapter 4. 

1.6 The Winter Anomaly: Theories and Experiments. 

The exceptionally high absorption of HF radio waves that is observed 
in some winter days at middle latitudes is known as the winter anomaly. 

It was first observed by Appleton [1937] and later confirmed by 
several workers [Appleton and Piggott s 1948, 1954; Dieminger , 1952]. 

The causes of the winter anomaly remain unexplained until now. 

Several theories have been suggested, however, to explain its existence. 
These theories are basically related to: 

a) An increase in electron production due to an increased 
concentration of NO; 

b) Decreased electron loss by a change in the dissociative 
recombination process; 

c) Increased electron production due to an increase in 
precipitated energetic electrons. 

The theories related to items a and b above assume a meteorological 
origin to the winter anomaly. Studies of Bossolasco and Elena [1963] and 
Gregory [1965] have demonstrated that the increases in P-region absorption 
are correlated with increases in the temperature at the 10 millibar (30 
km) level of the stratosphere. 

Geisler and Dickinson [1968] suggested that the winter anomaly would 
be produced by enhancements of NO concentrations, and that such enhance- 
ments would be caused by vertical transport resultant of planetary waves. 
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A second possible meteorological cause of the winter anomaly is 
related to variations in electron loss rates. Reid [1970] and Sealwist 
[1970] arrived at the conclusion that the principal electron loss 
mechanism around 80 km is probably recombination with molecular or 
hydrated ions. Seehriet [1970] suggested that changes in mesospheric 
water vapor content could alter the concentration of hydrated ions, and 
consequently the electron loss rate. Reid [1970] suggested that besides 
the 37 + ions observed by Nareisi and Bailey [1965] there are heavier 
hydrated ions in the mesosphere, as 55 + , 73 + , etc., that are not observed 
in rocket measurements because they are fragmented by the passage of the 
rocket. More complex ions would have higher recombination coefficients, 
since they have more degrees of freedom. Increases in mesospheric 
temperature would be sufficient to break up the larger ions, reducing the 
electron loss rate. 

Maehlum [1967] showed correlations between satellite measurements of 

precipitated energetic electrons and days of high absorption, and 

advocated the theory that the winter anomaly is produced by precipitated 

electrons. Gettev and Seohnist [1971] calculated the electron production 

produced by precipitating electrons necessary to explain the observed 

electron densities during days of high absorption and concluded that a 

very localized production rate of electrons would be required, and that 

3 - 1 

the production rate at 80 km should be over 100 cm s . That is not 
probable . 

Measurements of the winter anomaly have been in most cases obtained 
only by absorption measurements. Belvose [1966] obtained electron-density 
profiles during the winter by means of a partial-reflection experiment. A 
rocket electron-density profile during a day of high absorption was 



obtained by Meahtly and Smith [19681. The complexity of the factors 
involved in the winter anomaly, however, makes necessary simultaneous 
measurements of several parameters. A coordinated ground-based and rocke 
experiment was carried out in January and February of 1967 at Wallops 
Island, Virginia [Seohrist et at. , 1969] . The experiment included ground 
based ionosonde and radio-wave absorption measurements to determine the 
days of high and low absorption, and rocket measurements of temperature, 
winds and electron density. The electron-density profiles obtained by 
Sechrist et al . [1969] are shown in Figure 1.5. 

Figure 1.6 shows the rocket-grenade temperature and winds for the 
normal day (February 3, 1967) and for the anomalous day (January 31, 
1967). 

The electron-density profile for the anomalous winter day was very 
similar to electron-density profiles obtained during summer for the same 
zenith angle. The temperature profiles on normal and anomalous days were 
significantly different, with a warming occurring on the anomalous day 
above 70 km. The results of the experiment supported a meteorological 
origin to the winter anomaly. 

r 

During winter of 1970-1971, a coordinated rocket experiment was per- 
formed by the members of the Ionosphere Research Laboratory of 
Pennsylvania State University and the Atmospheric Sciences Laboratory of 
White Sands Missile Range, New Mexico [Mitchell et al., 1972]. The 
experiment intended to measure electron and positive- ion densities, 
neutral temperature, winds, ozone concentrations and neutral air 
densities on three different days, two of them preferably anomalous, and 
one a control day. The measurements were performed on January 22, 1971 
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Figure 1,5 Rocket electron-density profiles comparing summer values 
with normal and anomalous winter electron densities. 
(From Seohvist et al. [1969]). 
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Figure 1.6 Rocket -grenade temperatures and winds for a 

normal (February 3, 1967) and anomalous winter 
day (January 31, 1967) at Wallops Island. 

(From Seohrist et at. [1969]). 
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and January 26, 1971, anomalous days, and on February 1, 1971, the control 
day. Mitchell et at. [1972] suggested that the anomalous absorption on 
January 22 was probably related to meteorological phenomena, and 
apparently was characterized by a decrease in the electron loss rate, and 
the other anomalous day was possibly related to a solar disturbance, and 
was accompanied by an increase in detachment. 

During the 1971-1972 winter the University of Illinois conducted 
another coordinated experiment for measurements of the winter anomaly at 
Wallops Island. The description of the experiment and analysis of some of 
the results obtained will be made in Chapter 5. 

1.7 Statement of the Problem 

As discussed in the sections above, for a better understanding of the 
^-region processes it is necessary to develop new techniques of 
measurements, and to improve the accuracy of the existing techniques. One 
of the most important parameters to be measured is the electron density, 
in particular its time and spatial variations, and the small-scale struc- 
ture of electron-density profiles. The partial-reflection technique, 
being a ground-based technique, and relatively inexpensive, is 
particularly suitable to this type of measurement. The accuracy and 
limitations of the partial-reflection technique, however, are not very 
well known at the present time. The real structure of the irregularities 
producing the reflections is not known. The range of validity of the 
scattering theories and signal processing methods used, for example, must 
be examined critically. The appearance of minimums in electron density 
between 75 and 80 km is another point that deserves attention. 

It is the purpose of the present work to make an evaluation of the 
partial-reflection technique, to analyze the range of validity of the 
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theories used, to extend the scattering theory to a locally homogeneous 
random medium, to verify the possible reason for disagreement between 
rocket and partial-reflection profiles, to improve the signal analysis 
techniques, and to present and discuss results of partial-reflection 
measurements made during the winter 1971-1972. 

In Chapter 2 the specifications and signal processing methods for a 
partial-reflection system will be discussed and the system in use by the 
University of Illinois will be discussed in Chapter 3. 

In Chapter 4 a theory of scattering in a locally homogeneous random 
medium, will be developed. A medium like this is observed above 80 km, 
where the large gradients in electron density make the D region depart 
considerably from a homogeneous medium, as is normally assumed in partial- 
reflection measurements. Comparisons of the developed theory with the 
classical theories used in partial reflection will be made. The possible 
influence of the scattering theories in the electron-density profiles will 
be discussed. The influence of a strong dependence on height of the 
scattering cross section will be studied. 

In Chapter 5 results of partial-reflection measurements obtained 
during the winter 1971-1972 will be presented and discussed. 
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2. DESIGN CONSIDERATIONS FOR A PARTIAL- REFLECTION SYSTEM 

In this chapter will be discussed the specifications of the trans- 
mitter, receiver, antenna system and signal processing methods, adequate 
to partial-reflection measurements. The nature of the received signals 
will be presented initially, and based on its characteristics, the desired 
specifications of the system will be analyzed. 

2.1 The Nature of the Received Signals 

Partial-reflection signals are the result of small reflections, 
produced at heights between 50 and 90 km altitude. The incident signals 
are HF pulses, of w r idth between 10 and 50 Usee, with polarization corre- 
sponding to the ordinary and extraordinary modes of propagation. The 
reflection coefficient of the ionosphere is very small in this height 
range, going from approximately 10" 6 at the lower altitudes to 10~* 
near 90 km. Measured values of the reflection coefficients will be 
presented in Chapter 5. The number of reflections observed on each 
sample is normally small between 60 and 75 km, and increases above 
this height. As a result, isolated reflections can be observed below 
75 km, but the amplitude of the signals increases steadily with height 
above approximately 80 km. Shown in Figure 2.1 is a sample of the 
received signals for the ordinary and extraordinary modes of propagation. 

The small reflection coefficient below 75 km, and the strong absorption 
for the extraordinary mode above 80 km, make the received signals very 
weak, and some signal processing system able to detect signals with a 
signal -to-noise ratio close or below 1 must be used, in order to obtain a 
better utilization of the received data, principally at low altitudes. 

The dynamic range of the amplitudes of the signals is very large, 
the signals being, on the average, at 80 km, approximately 20 times greater 



Figure 2.1 Sample of received signal. The upper trace 
corresponds to ordinary polarization and the 
lower trace corresponds to extraordinary 
polarization. The vertical center line of 
the graticule is the 60 km mark, with height 
increasing to the right; each centimeter repre 
sents 15 km [Pirmat and Botihi.ll, 1968]. 
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than at 60 km. Besides that, the fading of the signals produces a large 
fluctuation of amplitudes at the same height. As a result, amplitude 
fluctuations as large as 50 dB must be detected. Such fluctuations impose 

some restrictions on the specification of the receiver, that will be 

> 

discussed in Section 2.3. A sample of the variations of amplitude of the 
ordinary and extraordinary signals for reflections from 82 km altitude is 
shown in Figure 2.2. It is observed that the signals present a fading 
with a correlation time of approximately 4 sec. 

The correlation coefficients, between the A and A signals at the same 

O' 

height, and between A signals at heights 1.5 and 3.0 km apart, are shown 
in Table 2.1, for the data obtained on January 28, 1972, at 12:00 hours, 
local time, at Wallops Island, Va. The average A and A amplitudes, the 
noise amplitude, and the electron densities, obtained from the same data 
used in Table 2.1 are shown in Table 2.2. The noise amplitude is the result 
of averaging the power of noise from samples collected between the heights 
of 40 and 44,5 km, where partial reflections are not present. It is observed 
thatiat low altitudes the A and A signals are very well correlated. The 

O X 

* >9 • * • /• f 

correlation decreases with altitude, and is very poor for altitudes above 
80 km. Such low correlation coefficients cannot be explained by the pre- 
sence of noise, since the signal-to-noise ratio at 81 km, for example, is 
approximately 5, as shown in Table 2,2. Austin [1971] suggested that the 
low correlation coefficients at higher altitudes could be the result of the 
different group velocities of the ordinary and extraordinary modes. The 
height of reflection of the signals is determined by assuming that both 
signals propagate with the velocity of light in free space. As a result, 
at higher altitudes, reflections considered from the same height, are 





TABLE 2.1 

Correlation coefficients between A 0 and A x , at the same height, and 
between A 0 signals separated by 1.5 and 3.0 km, for the data obtained 
on January 28, 1972, at 12:00 hours, local time, at Wallops Island., 
Virginia. 
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Correlation coefficient 
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A (h) -A (fc- 1.5km) 


A Ah') -A (h- 3km) 





TABLE 2.2 

A Qy A x and electron densities for the data obtained on January 
28, 1972, at 12:00 hours, local time, at Wallops Island, Va. 

No sample received was rejected or reached saturation of the 
receiver. 

Height 

Average 

Average 

Electron density 




„ -3. 

(km) 

^ o 
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63 

10.0 

15.7 

1.19xl0 8 

64.5 

10.2 

16.3 

1.20xl0 8 

66 

9.9 

16.3 

: 4.90xl0 ? 

67.5 

10.6 . 

18.2 

1.32x10 s 

69 

13.3 

22.8 

; 1.44xl0 8 

70.5 

17.2 

28.5 

1.25xl0 8 

72 

19.5 

31.6 

O 

1.3.4x10 

73.5 

19.7 

30.9 

1.56xl0 8 

75 

18.8 

27.7 

3.08xl0 8 

76.5 

17.8 

22.4 

6.15xl0 8 

78 

19.3 

17.4 

6.41xl0 8 

79.5 

22.5 

15.2 

5.22xl0 8 

81 

25.6 

14.2 

l,31xl0 9 

82.5 

36.1 

13.1 

2.99xl0 9 . 

84 

74 . 0 

12.5 

. 2 . 10x10° 

85.5 

123.4 

13.7 

5.85xl0 8 

87 

88.5 
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164.3 
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182.4 

: i 
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19.4 

20.2 
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Average noise 
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really produced from different heights. The error in the determination of 
height, however, is much less than 1.5 km, as will be shown in Chapter 4; 
Table 2.2 shows that above 80 km, the A o signals for heights separated by 
1.5 km are very well correlated. So the justification of Austin for the 
low correlation coefficients between A and A ^ does not seem reasonable. 

Von Biel et al, 3 [1970] calculated the correlation coefficient between the 
ordinary and extraordinary reflections that would be produced by volume 
scattering, instead of a sharp reflector, and showed that in this case the 
correlation coefficient can assume very small values. A second indication 
of reflections produced by volume scattering is the good correlation co- 
efficient between A q signals separated by 1.5 or. 3.0 km, as shown in- 
Table 2.2. 

2.2 Sources of Errors in the Determination of the Electron-Density 

Profiles. 

The sources of error in partial -reflection measurements will be dis- 
cussed in this section in order to specify the system in a way of minimizing 
such errors. Errors directly produced by an inadequate choice of the model 
of reflector were shown in Chapter 1 and will not be discussed further. 

The expression for the determination of the average electron density 
between heights hj and is given by equation (1.10) that will be repeated 
below: 


io g [cvV2 /( Wi ] - 

- V 


( 2 . 1 ) 


where 
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R 


a)-a) r 2 r w-Wjr 2 
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and 



0)±{J) r 


(2.3) 


Errors are produced by iraprecisions in the measurement of the ratio 
(j4 J A XJ 04 / A )., in the determination of (F /F )_/•(/? /R ) 1 , that depends 
on the choice of the model of reflector and on the. previous knowledge of 
the collision frequency profile, and in the choice of the collision fre- 
quency profile, that alters both the factors n / {RJR^) ^ and (F^-F^) . 

Imprecisions in the measurement of (4 /A ) ? /(A l A ), are the results of 

CO O ^ CC O 1 

errors introduced by noise, or by wrong interpretation of the received data. 

I 

This second class of errors will be discussed in. Chapter 4. 

r 

The percentage of error produced by the factor (F / R ),/(F / R ). can 

4 ^ i y Ld C / 4- 

be expressed by: 


M _ A R 

N ~ [log(F) - log(4) ]F 


where A - (VV2 /( VVl and R = ' ( W2 /C W 1* 

It is observed that the percentage of error increases as the factor 

log(F) - log(4) decreases. Such a situation occurs at low altitudes where 
the electron density is small. Consequently, electron-density calculations 
become very imprecise at lower altitudes, principally below 60 km. If the 
frequency of the HF pulses is increased, the denominator in equation (2.4) 
decreases, and the percentage of error M/ff increases. The tolerable 
error &N/N, assuming a, given error AF/F, imposes a restriction on the maxi- 
mum frequency to be utilized for measurements over a given range of altitudes. 
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The lower limit of the frequency to be utilized is imposed by the attenua- 
tion of the extraordinary signal, that will make its amplitude to reach the 
minimum value detectable by the receiving system utilized. As the attenua- 
tion of the extraordinary signal affects principally reflections from 
higher altitudes, and errors due to imprecisions in the determination of R 
are more pronounced at lower altitudes, a given frequency of operation will 
be suitable for measurements only over a given range of altitudes. For a 
complete measurement of the D region, from 50 to 90 km, more than one fre- 
quency should be used, for better results. 

To estimate the value of the error produced by imprecisions in the 
determination of the ratio [R /R ) J [R /#)., the following calculation was 

CC 0 ** CC O -L . 

performed. The electron densities and collision frequencies shown in 
Table 2.3 were taken as a reference. They are representative of rocket 
profiles during conditions of maximum and minimum of solar activity. It 
was assumed that reflections were produced by sharp reflectors spaced by 
2.0 km, and the resultant R JR q profile was calculated, for frequencies of 
2.66 and 5.0 MHz. The corresponding profiles were obtained from the 

collision frequencies listed in Table 2.3. From the (A fA ) and (R JR \ 
profiles, the electron densities were calculated, reproducing, evidently, 
the data of Table 2.3. After that, the ratio [R/R)„J{R R ). was 

35 O £ 30 O X 

changed by ± 2, 4 and 10%, and the error introduced in electron density 
was verified. The results are shown in Figures 2.3 to 2.6. As expected, 
the error increases at lower altitudes, and for measurements below 70 km, 

(i? /R )_/(i? /R ) , must be determined with a precision of at least 2%. As 
can be observed in Figure 2.6, the frequency of 5.0 MHz is impracticable 
for partial -reflection measurements, even at heights of 90 km, for condi- 
tions of minimum solar activity. 
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TABLE 2.3 

Electron densities and collision frequencies used as a reference for the 
calculation of errors in electron-density measurements produced by impre- 
cisions in the determination of the ratio (R X /R Q ) 2 / 1 > as shown in 

Figures 2.3 to 2.6. The electron densities of column (a) correspond to 
conditions of maximum of solar activity, and column (b) to minimum of solar 
activity. 


Height 

-3 

Electron density (m ) 

Col. Freq^ 

(km) 

' (a) 

(b) 

-1 

sec 

60 

2 . 5xl0 7 

1 . OxlO 7 

2.5xl0 7 

62 

3.5xl0 7 

1 . OxlO 7 

1.7xl0 7 

64 

4 .OxlO 7 

2. OxlO 7 
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66 

6 . 0x1 0 7 

3. OxlO 7 

. 7.5xl0 6 
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3. OxlO 7 

5.5xl0 6 

70 

' 2.0xlO S 

4 4 OxlO 7 

4.5xl0 6 
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5. OxlO 7 
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6. OxlO 7 

2 . 5x1 0 6 
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6. OxlO 8 
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1.0x10 
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80 
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i 

9. OxlO 5 

82 ' 

1 .OxlO 9 

1 . OxlO 8 

6.0xl0 5 

84 

4. OxlO 9 

3. OxlO 8 

4.5xl0 5 

86 

5. OxlO 9 

6. OxlO 8 

3. 5x1 0 5 

88 

8. OxlO 9 

3. OxlO 9 

2.5xl0 5 

90 

l.OxlO 10 

6. OxlO 9 

1.5xl0 5 
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Error in electron density produced by errors of ±2, 

4 and 10% in the ratio (R x /R 0 ) / (R x /Ro ) . , for the 
frequency of 2.66 MHz, electron densities and colli- 
sion frequencies listed on Table 2.3 (column a), 
corresponding to conditions of maximum solar activity. 
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Error in electron density produced by errors of ±2, 4 
and 10% in the ratio {R x /R ( ^ 2 /(/?,,/# 0 ) for a frequency 

of 2.66 MHz, electron densities and collision frequencies 
listed on Table 2.3 (column b) , corresponding to conditions 
of minimum of solar activity. 
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Figure 2.6 Error in electron density produced by errors of ±2, 4 and 
10% in the ratio (V^o)2^Wl» for a frequency of 
5.0 MHz, electron densities and collision frequencies 
listed on Table 2.3 (column b) , corresponding to conditions 
of minimum of solar activity. 
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The error produced by the ratio {AjA^^ is given by: 


LN _ -M 

N ~ [log(S) - log(4)]4 


(2.5) 


This expression is similar to equation (2.4), that gives the error produced 
by (i? /R ) /(i? /R ),, so that the same discussion made above applies to the 
present case. It should be noted, that for a frequency of operation of 
2.66 MHz, and at lower heights, the error in the determination of 
(4 /A ) / {A /A ), has to be less than 2%. Such tolerance will impose 
restrictions on the signal processing method to be used, as will be dis- 
cussed in Section 2,6. 

An expression for the error introduced by imprecisions on the 
determination of the collision frequency, v, would be too involved, since 
v appears in the factors (i? /R)J(R/R), , and (F-F ) . The influence of 
the value of v in the calculation of electron densities can be better 
observed from an analysis of Figure 2.7. In Figure 2.7 a plot of electron 
density as a function of A is shown, using as a parameter the values of 
collision frequency at heights aiK * that are assumed to be separated 
by 1.5 km. If a given value of A is considered, the influence of the values 
of collision frequency on the calculation of electron densities can be 
obtained by taking the values of electron density corresponding to different 
curves, each one corresponding to a pair of values of v, at heights h ^ and 

It is observed, that even if the collision frequency is changed by a 
large amount (from curves a to b, for example, the collision frequencies 
change from 2.4 x 10^ and 1.8 x 10^ s ^ to 1.7 x 10^ and 1.08 x 10^ s *), 
the corresponding variation of electron density is not very large. As a 



(Ax/AoJa/Ux/Ao)! 



Figure 2.7 Electron density as a function of the ratio {A x /A 0 )2/ 

for different values of collision frequency. 
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consequence, the choice o£ the collision frequency model to be used is not 
very critical. 

2.3 The Receiver 

The most stringent specifications for a receiver to be used in partial 
reflection are the dynamic range and the linearity of the receiver. Satu- 
ration of the signals can produce serious errors in the measurement, and as 
the amplitude of the reflections cover a very wide range, the receiver to 
be used must be linear over a range of at least 40 dB, with a tolerance in 
linearity of at most 10% „ Such specification imposes some limitations on 
the design of the intermediate stage, responsible for non-linearities of 
large input signals, and on the detector, responsible for non-linearities 
of small input signals. The minimum input signal to be detected is limited 
by atmospheric noise, and not by the receiver noise, and as a consequence, 
the sensitivity is not an important specification. 

The bandwidth is given by the expression 

bW = l/W (2.6) 

where W is the width of the transmitted pulse. Such width is limited on the 
lower side by the increase of the necessary bandwidth, with deterioration of 
the signal-to-noise ratio, and increase in interference produced by spurious 
frequencies. On the upper side, the width of the pulse is limited by the 
height resolution of the measurement. A pulse of 50 psec, for example, 
occupies a height range of 7.5 km, Pulse widths used in the existing partial - 
reflection systems range from 25 to 50 psec, corresponding to bandwidths of 
40 to 20 kHz, For a better height resolution, 25 psec or less should be 
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used all the times that the signal-to-noise ratio is good enough. A pulse 
width of 10 psec, corresponding to a bandwidth of 100 kHz should be tried, 
principally for a comparative analysis of electron-density profiles obtained 
with different pulse widths. A more flexible partial-reflection receiver 
should be able to have its bandwidth adjusted in order to operate with pulse 
widths of 10, 25 and 50 psec. Phase distortion must be kept in a low level, 
in order to avoid distortion of the pulse, so that single tuned circuits are 
recommended in the design of the IF stage. 

The frequencies utilized on partial-reflection measurements are very 
close to the broadcast and commercial services bands. To avoid interference 
from such frequencies, a passband filter must be added to the front end of 
the receiver, with a bandwidth of approximately 400 kHz, and an attenuation 
of at least 50 dB for frequencies 200 kHz outside the bandwidth. 

The rejection of image and spurious frequencies of the receiver must be 
of at least 90 dB. 

2.4 The Transmitter 

The specifications to be discussed in the transmitter are the power 
output and the pulse repetition rate. 

The necessary power output can be determined from the following 
expression: 


P.-P+L , in dB 

trs 


(2.7) 


where P^ is the output power, is the received power and is the system 
loss, that is given by: 


L 

s 



+ L 
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+ L h 


+ R. 


( 2 . 8 ) 
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where L is the free space loss, given by: 

L = 20 log{^- s } C 2 - 9 ) 

p a 

f being the frequency, o the* velocity of light, s the ray path, is the 
ionospheric absorption, G_^ and G^_ the gains of the receiving and transmit- 
ting antennas, L ^ are the ohmic losses in the transmission lines, and R \ 
the reflection coefficient >f the ionosphere. 

The necessary received power P^ depends on the minimum signal-to- 
noise ratio that the system is able to process, on the atmospheric noise at 
the receiver location, and on the bandwidth of the receiver. • 

The power of noise is given by [Davies, 1965]: 

P = F + B - 204 dB/w (2.10) 

n a 

where F ^ is the external noise power available from a noiseless antenna, 
expressed in decibels above K2‘ , where K is Boltzmann's constant and 
T = 288.39 °K is a reference temperature, which is the noise generated 
in a unit bandwidth by a thermal source in a temperature f t and 
B = 10 log(b) , b being the bandwidth. The values of F ' are given in the 
literature (CCIR Repoi'tNc. 322). 

Based on the above expressions, the power necessary to a transmitter 
to be located at Urbana, Illinois will be calculated below. 

The following parameters will be considered: 

Receiving and transmitting antenna gains: 20 dB, each one; 

Minimum signal-to-noise ratio acceptable: 1; 
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Ohmic loss in the transmission lines: 3 dB; 

Frequency of operation: 2.66 MHz; 

bandwidth: 40 kHz (25 ysec pulses). 

Since the smallest reflections are produced in the height range of 60 
to 70 km, a height of 65 km will be used in the calculations, resulting in 
a free space loss L = 92 dB, 

The noise factor F , for winter conditions, during the morning, is 

a 

23 dB [Davies j 1965], Assuming that a reflection coefficient of the iono- 
sphere of 3 x 10~6 must be detected, a transmitter power of at least 100 kw 
must be used. 

The pulse repetition rate ideally should be equal to the inverse of 
the correlation time of the noise, that is equal to the inverse of the band- 
width, assuming a white noise at the input. Two practical reasons, however, 
limit this ratio: 

1) As the repetition rate increases, the average power to be delivered by 
the transmitter increases; 

2) Between the reception of the ordinary signal and the transmission of 
the extraordinary signal, a switching device must change the polariza- 
tion of the antennas from one mode to the other; the time of operation 
of such a device is a limiting factor in the pulse repetition rate. 

Taking into consideration the above restrictions, a reasonable pulse 
repetition rate can be considered as 50 double pulses per second. Such a 
rate, for a transmitter with a power output of 300 kw, and a pulse width 
of 25 ysec, implies an average power output of 750 w. 

2.5 The Antenna System 

The directivity of the antennas is an important factor in the design 
of a partial-reflection system, in order to avoid a strong contribution of 
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oblique reflections to the received signals. Assuming that all the trans- 
mitted power is uniformly distributed over the main lobe of the antenna, 
the reflections received from a height h will be produced by a volume with 
an average width of 2?jtg(6/2), where 6 is the width of the main lobe, and 
a height cW/2, where a is the velocity of light, and W the pulse width. 

The indetermination produced by the pulse width, in height, is equal to 
oW( 2, and the indetermination in height produced by the finite angle of 
radiation of the antenna is given by h [1-cos (9/2) ] , as can be observed in 
Figure 2.8. If the main lobe of the antenna is of ‘20°, for reflections 
from 80 km altitude and a pulse width of 25 usee the scattering volume 
will have dimensions of 2.75 km of height and 28 km of width. The inde- 
termination in height produced by the lobe of the antenna will be 1.3 km, 
smaller than the height of the scattering volume, that depends only on the 
width of the transmitted pulse. An antenna lobe of 20°, in this case, is 
enough for partial-reflection measurements, and will keep contributions of 
oblique reflections at a low level. 

, An attenuation of secondary lobes of at least 20 dB must be obtained, 
in order to reduce the reflections produced by such lobes by at least 40 dB, 
if both the receiving and transmitting antennas have the same characteristics. 

The antennas must operate with the polarizations corresponding to the 
ordinary and extraordinary modes of propagation. A switch device must be 
incorporated to the antenna system in order to change the polarization from 
one mode to the other. Since the polarizations used are circular, an 
array of cross dipoles is the most convenient type of antenna for partial 
reflection. 

In the design of the antennas, the same antenna can be used for trans- 
mission and reception, or independent antennas can be used. The first 
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solution implies the use of a decoupling device between the receiver and 
the transmitter. An example of such a device is the transmission line 
bridge shown in Figure 2.9 [Westcott, 1948]. The use of a single antenna, 
although more economical, gives less flexibility to the system. If two 
antennas are used, there is the possibility of transmitting one mode of 
polarization and receiving a different mode. Such type of operation is 
convenient for phase measurements [Von Biel et al. 3 1970], and for measure- 
ments of rejection of the undesired mode of propagation of the antennas. 

2.6 The Signal-Processing System 

The relatively high correlation time of the reflections (4 to 5 sec, 
approximately), permits the use of a series of signal processing methods, 
in order to improve the effective signal-to-noise ratio of the received 
signals. Due to the limitations of the partial-reflection systems in 
operation, however, the advantages of such methods have not been fully ex- 
ploited to the present time. Most of the processing techniques used in 
Thomson-scatter experiments, for example, could be adapted to partial- 
reflection measurements. In the system in operation at the University of 
Illinois, some techniques have been introduced, by using the correlation 
function of the signals [D a R. Ward and S. A . Bowhill , private communica- 
tion], and integration of the signals, as in the present work. The inte- 
gration method will be discussed in Chapter 5. 
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3. THE PARTIAL- REFLECTION SYSTEM AT THE UNIVERSITY OF ILLINOIS 

3.1 Charaeteristies of the Equipment 

The partial-reflection system in operation by the University of 
Illinois, at Urbana, at the time the present work was made, is described in 
detail by Henry [1966] and by Pimat and Bowhill [1968], and had the fol- 
lowing characteristics: 

Transmitter: 

Peak power: - 50 kW 

Frequency: - 2.66 MHz 

Pulse width: - adjustable between 20 and 50 ysec 

Output impedance: - 50 ohm, unbalanced 
Receiver: 

Noise figure: - 3 dB, maximum 

Selectivity: - 25 kHz between -6 dB points, as shown in Figure 3.1 

Linearity: - according to Figure 3.2 

Recovery time: - 200 msec for the receiver to drop into noise after 
0.1 volts RMS at the signal frequency applied at the input is 
removed 

RF input impedance: - 50 ohm, unbalanced. 

As can be observed in Figure 3,2, the linearity of the receiver is not 
good enough for partial-reflection applications. To compensate for the 
effects of non-linearity, two subroutines, VALUE and LINAP, were included 
in the computer program that calculated electron density. Such subroutines 
are listed in Appendix II, The subroutine VALUE is a list of the values of 
input signals at the receiver, TU(I), the corresponding values of output 
signals, TUO(I), and the ratios of increments at the input and output. 








S(I) . The subroutine LINAP uses the data listed in VALUE to transform 
output signals into the corresponding input signals, that are used in the 
calculations of electron densities, eliminating in this way the effects of 
non-linearity of the receiver. The subroutine LINAP employs a linear inter- 
polation between two of the values tabulated in VALUE to perform the 
linearization. 

The dynamic range of the receiver is not enough for measurements of 
electz’on densities between 60 and 90 km, saturation generally occurring for 
reflections above 84 km. To overcome such a problem, a programmed attenua- 
tor, adjustable between 0 and 30 dB, was introduced at the input of the 
receiver. This attenuator operates only on alternating frames, and when 
it is used, the output data consists of two interconnected sets of data, 
one with the attenuator on, and the other with the attenuator off. The in- 
convenience of such a method is that only one-half of the measured samples 
are used in the calculation of electron density at each height. The com- 
puter program PROAT calculates electron densities using the data obtained 
as described above. 

Using the above method of data collection it was possible to obtain 
electron-density profiles between the heights of 70 and 90 km, without 
saturation of any sample. Figure 3.3 shows two electron-density profiles 
obtained in this way. 

The antenna system is composed of two independent arrays for trans- 
mission and reception. Each array consists of two sub-arrays of 30 
parallel half-wave dipoles. The sub-arrays are prependicular to each 
other. A schematic diagram of the two antennas is shown in Figure 3.4. 
Switching between the ordinary and extraordinary modes of polarization is 
controlled by. the polarization reversal control, operated by a pulser. 



HEIGHT (km) 


OCT. 9, 1971 
8:50 LOCALTIME 

OCT. 17, 1971 
8:25 LOCAL TIME 


60 L 
10 7 



ELECTRON DENSITY (m*3) 


Figure 3,3 Electron-density profiles obtained with a programmed attenuator at the input of the 

receiver, at the University of Illinois, Urbana, Illinois, on October 9, 1971 at 8.50h 
local time (a), and October 17, 1971 at 8:25h local time (b). 





Figure 3.4 Antenna layout (From Henry [1966]) 
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For the operation of the partial-reflection system at Wallops Island, 
Virginia, the receiver and transmitter were installed in a van, and a new 
antenna was constructed. A single antenna was used, consisting of two 
perpendicular sub-arrays of four short dipoles. Isolation between trans- 
and receiver was obtained by means of a transmission line bridge, 
as shown in Figure 2.9. At the input of the receiver, as an additional 
protection, it was placed the diode circuit shown in Figure 3.5. This 
diode circuit is operated by a pulser, and opens the input of the re- 
ceiver all the time, with exception of the time interval between 20 and 
300 ysec after a pulse is transmitted. 

The measured gain of the antenna used at Wallops Island [G. W. 

Henry, Jr., private communication] was approximately 16 dB. 

3.2 The Signal Processing System 

In the system operated by the University of Illinois, the output of 
the receiver is coupled to an analog-to-digital converter, that for each 
reflection takes 21 samples of signal, from 60 to 90 km, with a height 
interval of 1.5 km, and four samples of noise, corresponding to the height 
interval between 40 and 44.5 km, where reflected signals are absent. The 
sampled signals are fed to a PDP-15 digital computer, and stored on mag- 
netic tapes. The sampling and storage systems are described in detail by 
Birley and Sechrist [1971], 

The data are processed by a computer program, PROAX, that is listed in 
Appendix II, The processing system is based on taking the average power of 
signal at each height, and subtracting it from the average power of noise, 
that is obtained from the samples measured in the interval corresponding to 
the heights between 40 and 44.5 km. In this way, the average signal -4^ ^ (h) 



FROM 

TRANSMISSION 



FROM PULSER 


Figure 3.5 Diode circuit used to protect the receiver from overload. 


ON 
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from a height h is given by: 


A (70 = 
o } x 


i - 1 

n n_ 

N 


n 


(3.1) 


where N is the number of samples taken, a Q x 00 n the amplitude of the 
ordinary (extraordinary) signal received from a height h, in the n-th 
sample, n is the average amplitude of noise in the n-th sample, for 

the ordinary (extraordinary) frame, and is given by: 


n 

o 3 cc 3 n 

a corresponding to the m-th sample of the four samples of noise 

o,x 3 n 3 m 

taken during each frame. 

It was observed that besides the Gaussian noise, an impulsive noise, 
probably produced by nearby power transmission lines, appeared on the 
signals o An example of such a noise is shown in Figure 3.6* 

The integration procedure expressed by equation (3.1) is effective in 
eliminating a Gaussian noise, but not an impulsive noise. 

The impulsive noise observed is characterized by a rapid fading, and 
normally does not appear on successive frames at the same height. This 
characteristic was used to minimize its influence on the averaging of the 
signals, in the following way: each sample from a given height, a 0 ^ x ^ n 

is compared with the sample from the same height in the frame immediately 
before, a If the difference - a 0yX ^ n -l is Sweater 





Figure 3.6 Sample of received impulsive noise. 

The horizontal scale is of 0.2 msec 
per graticule. 
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than a specified value, that is an input of the processing program, for the 
height h, h - 1,5 km, or h + 1,5 km, the sample is rejected, since such a 
rapid fading corresponds probably to the presence of impulsive noise. 

Too noisy frames, produced principally by lightning, are also rejected. 
If the average noise on two consecutive frames, one ordinary and the other 
extraordinary, exceeds a given value, both frames, for all the 21 heights, 
are rejected. 

The signals used in the above computations are linearized signals, 
obtained from the subroutines VALUE and LINAP, described before. 

Once the average A and A are obtained, electron densities are calcu- 
lated by calling the subroutine CALC, that is listed in Appendix II. 

The number of samples necessary to the calculation of an electron- 


density profile can be obtained from the expression of the statistical 


uncertainty associated with the total power A ( h ) at a given height, 

that is given by: 



2 

where n is the average power of noise, and N the number of samples taken. 
Using as the minimum signal-to-noise acceptable 1, and an uncertainty in 

2 o 

the total power A (h) = 0.04 A“ (K ) , it results N = 5000 samples. 

O j O 

3.3 Results Obtained 

Using the system and the processing method described above, it was 
possible to obtain reliable electron-density measurements down to a height 
of 63 km, taking approximately 5000 samples for each profile. Figure 3.7 
shows two electron-density profiles obtained in this way, on August 8, 1971, 



AUGUST 8, 1971 
9:20-9:50 LOCALTIME 

AUGUST 8, 1971 
11:00-11:30 LOCAL TIME 


ELECTRON DENSITY (m“3) 


Figure 3*7 Electron-density profiles at low altitudes obtained at Urbana, Illinois 
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between 9:20 and 9:50, local time, and on August 8, 1971 between 11:00 and 
11:30, local time, at Urbana, Illinois. In Table 3.1 the corresponding 
values of average noise, number of samples taken, A (h) and A (h) are shown. 
In Chapter 5 a more detailed analysis of the results obtained will be made. 
3.4 Suggested Improvements in the Urbana System 

The main limitation of the system in operation at Urbana, is the poor 
linearity of the receiver. The subroutines VALUE and LINAP that correct 
the non-linearity are relatively slow, and most of the processing time is 
taken by them. In consequence, the processing of an electron-density 
profile with 5000 samples takes approximately 30 minutes of computer time. 

A more linear receiver, besides saving processing time, and permitting an 
increase in the number of collected samples, would permit obtaining profiles 
from 60 to 80 km without the use of a programmed attenuator, if the dynamic 
range of the receiver is increased to at least 40 dB. The necessary modi- 
fications on the receiver implies a new design of the IF and detector 
stages. 

A second important improvement in the system would be an increase of 
the pulse repetition rate. In the system used now, the pulse repetition 
rate is limited by the average output power of the transmitter in two 
double-pulses per second. If the power output of the transmitter were 
increased, the repetition rate could be increased up to 20 double-pulses 
per second, this second limitation being imposed by the recovery time of 
the switching device that changes the polarization of the antennas from 
the ordinary to the extraordinary mode. A system with a speed of 20 
double-pulses per second would be able to take 12,000 samples in 10 minutes, 
improving the precision of the measured profiles, principally at low 


altitudes. 



TABLE 3.1 


Values of A x , A Qi A X /A Q average noise and number of samples taken, 
corresponding to the electron- density profiles of Figure 3.3. 


Height 

August 8, 71 9:20-9:50 h 

August 8, 

71 "10:00-11:30 h 

km 

A 

o 

A 

X 

A /A 
x o 

A 

o 

A 

X 

A JA 

X 0 

60 

17.5 

. 28.0 

— 

1.60 

30.2 

40.0 

1.33 

61.5 

20.1 

32.2 

1.60 

47.2 

66.9 

1.42 

63 

24.6 

39.7 

1.61 

75.1 

114.6 

1.53 

64.5 

31.4 

53.6 

1.71 

115.2 

180.5 

1.57 

66 

48.6 

80.6 

1.66 

159.8 

237.7 

1.48 

67.5 

71.9 

112.2 

1.56 

196.5 

270.4 

1.38 

69 

108.4 

147.1 

1.36 

223.0 

283.6 

1.27 

70.5 

146.4 

172.5 

1.18 

224.0 

284.4 

1.17 

72 

167.7 

174.3 

1.04 

244.9 

257.8 

1.05 

73.5 

179.9 

164.4 

0.91 

226.1 

200.1 

0.89 

75 

174.1 

: 131.7 

0.76 

218.5 

148.0 

0.68 

76.5 

165.7 

94.3 

0.57 

220.7 

111.7 

0.51 

78 

165.9 

71.4 

0.43 

223.7 

90.7 

i 

0.41 

1 

79.5 

200.0 

60.2 

0.30 

231.6 

79.5 

0.34 

81 

257.3 

57.8 

0.22 

258.6 

75.5 

i 

0.29 


Average noise: 

- 21.8 

Average noise: - 

27.8 


4824 

samples taken 

5040 samples taken 


51 frames rejected 

82 frames rejected 
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Finally, an increase of the peak power output of the transmitter to 
300 kw would permit an increase in the sensitivity of the system. 
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4. PARTIAL REFLECTIONS IN THE PRESENCE OF GRADIENTS OF ELECTRON DENSITY 

Observation of partial reflection and rocket electron-density 
profiles above 75 km shows the following characteristics: 

1) Partial-reflection measurements frequently present valleys in 
the electron-density profiles between 75 and 85 km, as can be verified on 
the profiles of Figures 4,1 to 4.3, obtained at Urbana, and on the results 
published by several experimenters [Belrose and Burke , 1964; Von Biel 

et al. 3 1970] . 

2) The A q profile in this height range shows a strong increase with 
altitude of the scattering cross section per unit volume. Figures 4.4 

to 4,6 show the A profiles corresponding to the electron-density profiles 
of Figures 4.1 to 4.3. 

3) Above 80 km partial reflections produce electron densities that are 
too low if compared to rocket measurements. 

4) Rocket profiles above 75 km are characterised by small changes with 
height up to 80 km where sharp gradients in electron density are observed 

[ Bowhill , 1969], 

Such differences between rocket and partial-reflection measurements 
will be investigated in this chapter. It will be shown that the gradients 
in electron density observed in rocket profiles can be responsible for the 
increase in the profiles near 80 km, and that in a region where the 
scattering cross section is a strong function of altitude, classical methods 
of partial- reflection data processing can introduce errors great enough to 
explain the observed differences. 

The gradients in electron density can give a two-fold contribution, to 


partial reflections: 



HEIGHT (km) 


MAY 7, 1971 

9:00 -9:25 h LOCAL TIME 
SAMPLES TAKEN -756 


ELECTRON DENSITY (m^) 


Figure 4.1 Electron -density profile obtained at the University of Illinois 
Urbana, 111., on May-7-71, at 9:00h; local time. Valleys of 
electron density are observed near 71 and 80 km. 
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Figure 4.2 Electron-density profile obtained at the University of 
Illinois, Urbana, 111., on August-11-71, at 7:30h, local 
time. A valley of electron density is observed near 80 km. 



HEIGHT (km) 


78 



Figure 4.3 Electron-density profile obtained at the University of 
Illinois , Urbana, 111. on August-11-71, at 8:05h, local 
time. Valleys of electron density are observed near 
70 and 80 km. 





HEIGHT (km) 


80 



Figure 4.5 A 0 profile corresponding to the electron- density profile of 
Figure 4.2. 
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1) The gradient itself produces reflections, due to changes in the 
refractive index; 

2) If there is turbulence, eddies with different electron densities 
will be mixed, and reflections will result from the random medium produced 
in this way. 

If such reflections are strong enough, they can be responsible for the 
changes in the scattering cross section near 80 km. 

The theories currently used in partial-reflection calculations do not 
take into account a medium as the one described above. Belrose and Burke 
[1964] assumed reflections produced by a sharp irregularity, and Flood 
[1968] assumed a random medium with a homogeneous background. 

The influence of the height dependence of the scattering cross section 
on the values of the calculated electron densities will be discussed in 
Section 4.1. It will be shown that the calculated electron-density pro- 
files depend on the width of the transmitted pulse, and that "fictitious 
valleys", very similar to that observed in the measurements, can appear 
on the profiles. The dependence of partial-reflection profiles on the 
pulse width has been already discussed by Coyne and Belrose [1973]. In 
Section 4.2 the reflection coefficients produced by gradients in electron 
density will be calculated, for gradients of the magnitude of that observed 
on rocket profiles. In Section 4.3 the theory of scattering in a random 
medium will be extended to the case of electron-density fluctuations in 
a locally homogeneous background medium. Results will be compared with 
the theories of Belrose and Burke [1964], Flood [1968] and Cohen [1971]. 
Possible influences of the correlation function of the irregularities on 

the ratio Itf /R I will be verified. The relative contribution of random 
1 x o' 
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irregularities and of gradients in electron density to partial reflections 
will be discussed. In Section 4,4 the nature of the irregularities pro- 
ducing partial reflections will be analyzed and in Section 4.5 some sug- 
gestions to overcome the influence of the finite pulse width on partial- 
reflection calculations, and experimental results, will be presented. 

4.1 Partial He flections in a Region where the Scattering Cross Section 

is a Function of Altitude, 

Observation of profiles, as shown in Figures 4.4 to 4.6 shows that 
the scattering cross section per unit height changes sharply with heights 
above 78 km. Such sharp increase in scattering cross section has been 
observed in practically all partial-reflection measurements [Gardner and 
Pawsey, 1953; Belrose and Burke t 1964]. 

The changes in scattering cross section with height are really 
greater than the changes observed in the profiles, since A q is the 
result of integration of reflections produced over a height interval cW/ 2, 
where c is the velocity of light in the medium, and W is the pulse width. 
To study the influence of changes in the scattering cross section on the 
calculated electron-density profiles, two models of scattering cross 
section per unit volume as a function of height will be considered: 

1) The scattering cross section per unit volume has a constant 

value cf. below a height h . a 0 between the heights h and h, , and o, above 
1 a L a d 5 

the height h^. Discontinuities are observed at h a and h as shown in 
Figure 4 e 7a; 

2) The scattering cross section per unit volume is a constant 
below a height h&> changes exponentially to a value at a height h and 
has a constant value o 3 above as shown in Figure 4.7b. 



HEIGHT (km) 


SCATTERING CROSS SECTION 
PER UNIT VOLUME (er) 


Figure 4.7 Scattering cross-section per unit volume 
profiles used on the calculations of 
Section 4.1. 
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For both models it will be assumed that the electron density and col- 
lision frequency, N and v, are constant over the region of interest. The 
A I A profiles that would result from a plane wave pulse propagating ver- 
tically in such a medium will be obtained, and, from such profiles, electron 
densities will be calculated using the theory of Belrose and Burke [1964] . 
The real electron density N will be compared to the calculated electron 
densities. 


4.1.1 Partial reflections in a region where the scattering cross 
section per unit volume changes in steps . If a plane wave pulse is inci- 
dent on the ionosphere, the amplitude of the reflected pulse will be given 
by 


where A~(h^) 

4n 


a 


D 

W 

h l 

o(h) 

c 

For the 

with h, - h 
b a 

-For h, < h 

l — a 


A Olj) = D 


i ^ He 

*i + r 


a(h) • exp 


. Wc 

h - ~ 


{-4a dh}dh 


(4.1) 


= reflected power 
= incident power 
= attenuation coefficient 

= average width of the volume V occupied by the pulse 
= width of the incident pulse 
= height of the center of the pulse 
= scattering cross section per unit volume 
= velocity of light in the medium. 

scattering cross section per unit volume profile of Figure 4.7a 
<_ f/e/4, it results 
- Wc/4 
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If both characteristic modes of propagation are transmitted, the ratio 
of the amplitudes of the reflected signals, A I A , can be obtained from 

SC o 

equations (4.2) to (4.6), using for a the following expressions, for the 
ordinary and extraordinary modes, respectively 


r N w+io r 

5 <3 s-t r xj -v 

CL ' ' C P- ( n ~ ) 

o 4 e mcv 5/2 v 


(4.7) 


N u~u> r 

d e n r jun 

a - C r , 0 ( J 

or 4 c 5/2 v - 


(4.8) 


The quasi- longitudinal approximation being assumed. 

The ratio of the scattering cross sections, c^/o , is given by 

CC 0 

Belrose and Burke [1964] 


uij. u>— ] 2 
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co— - ui t- ) 2 

(— k 


aj+(i) r oa+o)- 'j 2 r - 


2 5/2 ^ v 


u>+w r ) 2 
(— £•) 


(4.9) 


In the case considered here, the collision frequency is constant, and 

the ratio a /a will be constant. 
x o 

Figures 4.8 and 4.9 show the A and A X /A Q profiles for a pulse width 

9 -3 

of 50 ysec, an electron density of 10 m , a collision frequency of 
7.7 x 10 5 s” 1 , typical of heights near 80 km, and the following scattering 
cross sections per unit volume: 



HEIGHT (km) 


ha ~3f 


hn-6 


a ctj. = ct 3 = 0.01 <T 2 
b (Tj = cr 3 = 0. 5 erg 
C cr 3 = a-g =10001 



Aq (arbitrary units) 


Figure 4.8 Calculated A profiles, for a region where the electron density is of 10 m" , the 
collision frequency 7.7 x 10 5 s'i, and the scattering cross-section per unit volume 
profile is that shown on Figure 4.7a,with oj = 03 - 0.01 02 » ( a ) * a l = a 3 = a 2» 

and o„ = o„ = 100 a 
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1) 

C lo = 

o - = 0.01 

3c 

°2 c 

■ 

km 

2) 

°lo = 

a, = 0.5 

3o 

°2 0 

, h b -h a * 3 

km 

3) 

°2o = 

a 3o =100 

0 

• h - h a * 5 

km 


The index o refers to the ordinary mode. 

Electron-density profiles calculated from the AjA Q profiles of 
Figure 4.9 by using Belrose and Burke's [1964] theory are shown in 
Figure 4.10. It must be noticed that the real electron density is con- 
stant and equal to 10 m . As it is observed, variations of the scattering 
cross section with height, and the finite pulse width, can produce valleys 

in a region where the real electron density is constant. The depth of the 

_2 

valley can be as great as 2.5 x 10 “ of the real value, as in Figure 4.10a; 
even if the scattering cross section changes only by a factor of 2, as in 
Figure 4.10b (and in this case changes only by a factor of 1.25), the 
depth of the valley will be of 0.7 of the real value of the electron 
density. 

A comparison of Figure 4.10 with Figures 4.1 to 4.3 shows a great 
similarity between experimental results and model calculations. In 

g 

Figure 4.3, for example, the electron density changes from 6 x 10 at 

75.7 km to 1.6 x 10^ at 80.2 km; the average A changes from 52 to 75 km 
to 140 at 81 km. Such variations are comparable to the theoretical 
results of Figures 4,10c and 4.8c. 

To analyze the influence of changes of the pulse width in partial- 

reflections measurements, electron-density calculations were made for a 

9 -3 

region where the real electron density is 10 m , the collision frequency 

7.7 x 10^ s’ 1 , using pulses of 50 and 25 usee, for the following changes 
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in backs cattering cross section 

a) = °3 = °° 01 a 2 

^ 10 

b ) a 2 = °3 = 3" °1 

The results are shown in Figure 4.11. Figure 4.12 shows an experi- 
mental measurement, made by Belrose [1971], using the same pulse widths. 
Experimental and theoretical results agree very well, again. In both 
cases the depth of the valley is practically independent of the pulse width, 
but for the 25 psec pulse, the valley appears at higher altitudes, and is 
narrower. 

The above calculations show that the knowledge of A q profiles is 
necessary to the interpretation of partial-reflection results. Every time 
the average A q changes rapidly with height, valleys observed in electron- 
density profiles cannot be considered as real, but may be a consequence of 
changes in the scattering cross section and of the finite pulse width. 

4.1.2 The scattering cross section per unit volume changes exponen- 
tially as a function of height. The scattering cross section per unit 
volume to be considered is shown in Figure 4.7b and is given by: 

a(h) - Oj h < h a 

ff(h) = exp (6h) >_ h h a 

a Qi) = cfj h >_ hy 

The reflected power, as obtained from equation (4.1) is given by: 

-If h < h a - fcfc/4 

4 2 (h x ) a 

sinh(affc) exp(-4ah^) 

A • 
in 


(4.10) 



HEIGHT (Km) 



Figure 4.11 Electron-density profiles calculated using Belrose and Burke r s theory, for a region 

where the real electron density is of 10 ^ the collision frequency 7.7 x 1(P sec““^, 

the scattering cross section per unit volume profile is that shown on Figure 4.7, with 
a l - a 3 = 0.01 02 , pulse widths of SO usee, (a), and 25 usee, (2?), and with 02 * 03 *= 
Ailoi, pulse widths of 50 usee, (e) , and 25 usee, {d) . 
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-If h b - Wo/ 4 >_ fe ?_h a - Wo/ 4 

{h^) a.. r, 

— 2 55 4S" {exp(-4a(?: 1 - p ) ] - exp(-4a7z a )} 

A • 
in 


* (g ~ - " 4a) ex P^ “ + - exp[(g - 4a)fc a ]} 

-If h a + Wo /A >_ >_ - f/c/4 

Av a i i 

—2 = 4 ^ {exp[-4a(^ 1 -|%]-expC-4oi^ a )}+ -^^y {exp [ (e-4ot)ft fc }-exp[ (g- 

^in 


+ 4“ (exp(-4o* b ) - exp [-4a (ftp p)]} 

-If h ^ + Wo / 4 >_ + 

exp[(g-4a)7i & ) - exp[(S-4a) (?i 1 

+ ^ {exp(-4a - exp[-4a(h 1 - p)]) 
-If h l >_ h b + Jito/4 


4 2 (fcp 


m 



4 2 (^ 1 ) a 3 

— ^ ^ sin?z(cxftfc) exp(-4ah 1 ) 


AT 

in 


C4.ll) 


4a) k a ]} 


(4.12) 


(4.13) 


(4.14) 
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Equations (4.10) to (4.14) can be used in the calculation of electron 

densities as would be measured from partial reflections, using the same 

procedure as in Section 4.1.1. Figure 4.13 shows the electron-density 

9 -3 

profiles obtained for the case where the real electron density is 10 m , 
the collision frequency 7.7 x 10 5 s" 1 , the pulse width 50 usee, and the 
following scattering cross section per unit volume profiles: 

1) o lo = 1, 3 0 = 1.63 x 10~\ o 3o = exp(-B^) = 100 

h =0 h. - 3 km 

a b 

2) = ~ x 3 a 3 o ~ * 

h =0, h. ~ 3 km 

a b 

3) a lo = = 7,66 x 10 " 4 ’ °3o = 1 

h - 0, h-, = 3 km 

a b 

The results are essentially the same as obtained in Section 4.1.1, 
and they show that deep valleys in electron density can be produced even 
in the case where the scattering cross section per unit volume changes 
continually with height. 

4.2 Reflections Produced by Gradients in Electron Density. 

As discussed before, the scattering cross section of the ionosphere 
increases sharply at altitudes near 80 km. Such altitudes are charac- 
terized by gradients in the electron-density profiles, as can be observed 
in rocket measurements. On this item the reflection coefficient produced 
by gradients will be calculated, in order to verify if they have enough 
magnitude to explain the observed partial reflections. Initially it will 
be assumed that the electron density changes linearly with altitude, and 
a full-wave calculation will be applied to the determination of the 



HEIGHT (km) 



Figure 4.13 Electron- density profiles calculated using Belrose and Burke's [1964] theory, for a 
region where the real electron density is 10® the collision frequency 7.7 x 10“* 

sec"l, and the scattering cross-section profile is that shown on Figure 4.7b, with 
h a -hb - 3 km, and with a \ Q = 6 o = 1.63 x ID -3 , 03 ^ = exp(-ghfc) , $q = 7 *66 x 10 ' 4 
= 1 , (c) • 
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reflection coefficients. After that, an approximate solution will be 
applied to different geometries of electron-density profiles. 

4.2.1 Reflections from gradients with a linear variation in electron 
density. The model of reflector utilized in this item is shown in Figure 
4,14. The electron density is constant above a height h 2 and below a 
height and increases linearly with height between and h^ The col- 
lision frequency is assumed constant, the ionosphere horizontally stratified 
and the propagation quasi-longitudinal. For such a simple model it is pos- 
sible to find an analytical solution for the reflection coefficients, 
avoiding numerical calculations. 

As the propagation is longitudinal, and the ionosphere horizontally 
stratified, the solution of the wave equation can be written as [Budden 3 
1960] 

1) If h ± h 1 

A plane wave solution: 


E ' ~ exp{ - i (h-h^)} + exp{£ (4. 15a) 

H ol ~ ^ n l o ex P^~ ^ ^ n l o R o ex P^ ^O n lo^”^l^ (4.15b) 

1 = exp{- i exp{i k^n^Oi-h^)} (4.16a) 

H xl = 1 n lx exp{- i ' 1 n lx R X exp{i (4.16b) 


where fc n is the wave number in free space, R , is the reflection coef- 

U O 

ficient, ^ is the refractive index for heights equal or below h^, the 
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indices o and x referring to the ordinary and extraordinary modes of propa- 
gation, respectively. 

2) If h >_ h 2 

A plane wave solution: 


E 


o 3 


H 


o 3 


E 


x3 


H 


x3 


- exp{- k^i 2 )} 

(4.17a) 

= i n 2o C Q exp U k Q n 2x (h-h 2 -)} 

(4.17b) 

= C x exp{- i kQn^h-h^} 

(4.18a) 

~ in 2x C x ex P U 0 n 2 x^ 2 )} 

(4.18b) 


where n 2o x is the refractive index for heights equal or above h^ t the 
indices o and x referring to the ordinary and extraordinary modes of 
propagation. 

The solution to the wave equation in a medium with linear variation 
in electron density is given by Budden [1960], and will be extended below 
to include the Sen-Wyllev [1960] expressions for the refractive index. 

The refractive index of the medium for quasi-longitudinal approxima- 
tion, is given by equation (1.3) that will be repeated below 


01 


0)±U) T 


n 


o.x 


= 1 - 


- W 

. 5 o 


0)±01 T 


— C 3/2<V> - { 2 -V- C s,2<-T^ 

QJV 


(4.19) 


The height dependence of the electron density can be expressed by 

N = B(h - ft ) 
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and In being constants. As a result n~ can be written as 
o & o,x 


n = 1 ~ y (h-h ) 
OjX o J 


(4.20) 


where 


rCn>±w j 


'OjX 


(0±U) 


- 5 1 


R j ± p 

— L ' ^ r 


L 03V 


2 3/2 v 2 v 5/2 v 


~ ! 

2„ 

h i 

g 3 

V 

e jn 

_ 

0 


(4.21) 


As a consequence the wave equation has the form 


d 2 E 


dh 


f t£+ 


<r J 


- 0 


(4.22) 


Taking 


t = ( k 2 y ) 1/3 (h-h -l/y) 

OjX 0 0 3 X K o 


(4.23) 


and substituting equation (4.23) into equation (4.22) it results 


d 

d 


2 


E 

o 3 x 

2 

T 

o 3 x 


T E 
o 3 x o 3 x 


(4.24) 


Equation (4.24) is the Airy's equation whose solution A . (x ) and 

'l' O y X 

B .( t ) are discussed in Budden [I960]. 

'V O j X 

Using the solution A.(t ) and B .( t ), the electric and magnetic 

'l' OyX "l* OyX 

fields can be written as 


E 

O t X 


(4.25) 



102 


i * 21/^ 

H = - 1/k {B1 A.(t ) + B2 B . (t )}(fe ) 

' c> o,x i o,x o,x t- o 3 x o 3 x 


(4.26) 


Where 


4.- (X J = 


<3 4.(x ) 

x o.x 


i 1 - o.ar d x 


and 5 . ( x ) = 

-rl '- /O /V» r 


d B.(x ) 
x o.ar 


o,x 


0,ar i t 


o 3 x 


Applying the boundary conditions at h ^ and h^, using equations (4,15) 
to (4.18), (4,25) and (4.26) for the fields, the following expression is 
obtained 


1+1? n. [a A.(xl ) + B. (tl )] 

o^x _ lOjic 1 - OjX x o^x ^ o 3 x 


1 - R 


o 3 x 


y 


f(-^-){a 


a\ (xl ) + b\ (xl )} 

^ v ✓n rY*s < v iO ***** 


O^X 'L OyX' 


Q^X' 


(4.27) 


where 


lo,* 1 ' 3 ' 


- "2o.x V T W } 


o,x 


(4.28) 


Y o x 1/3 ' 

n _ „ 4,(x2^ J - i(-g^) 4,(x2^ J 


OyX V O y X 


OyX' 


ho.x ■ value of T o,x at helght *1 


T 2o,ar * Value ° f \,x at helght V 


The reflection coefficients, R and R can be determined from equa- 

’ o x n 


tions (4.27) and (4.28) 
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Plots of |i? [j | R | and ji? /R \ as a function of - h 1 are shown in 

O X X O £ A 

Figures 4,15 and 4.16, for the case where the electron density changes from 

9-39-3 5-1 

2 0 5 x 10 m to 3 x 10 in , the collision frequency is 2.15 x 10 s , 

the frequency is of 2,66 Mhz, and the distance h 2 - changes from zero to 

8 -3 9 -3 

150 m. If the electron density changes from 6 x 10 m” to 7 x 10 m 
over a distance of 2 km, values typical of rocket profiles near 80 km, for 
the same frequency and collision frequency as above, the results are 
R = -3,7 x 10" 6 + i 3.9 x 1(T 6 and R = -2.4 x 10" 4 - i 1.8 x 10" 4 . 

O X 

The results show that reflections produced by gradients in electron 
density are of the same order of magnitude of measured partial reflections. 
Resonances are observed, for distances close to a multiple of half wave- 
lengths. The situation is the same as in a transmission line with 
characteristic impedance changing as a function of length. The difference 
between the wavelengths of the ordinary and extraordinary modes produces 
minimums and maximums in \R /R |. If enough number of samples are taken, 

minimums and maximums tend to cancel, and the average )/? / R | will be the 

X o 

same as if = 0, that is, a sharp reflector as in Belrose and 

Burke's [1964] theory. The changes of the ratio |l? /R | as a function of 
the geometry of the irregularities show that in partial-reflection measure- 
ments it is necessary to take the average value of A q and A ^ over all samples 
received, without rejecting samples with low signal-to-noise ratio or samples 
that reach the saturation of the receiver. If only a few percent of the 
samples reach the saturation of the receiver, even if the ordinary and the 
corresponding extraordinary samples are rejected, a bias can be introduced 
on the average A and A producing errors in the calculated electron 
densities. Such a situation is exemplified in Figure 4.17. Figure 4.17a 

I 



|Rx|, |Ro| 
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Figure 4.15 ji?^j and |i?^| as a function of the height interval h 7 -h^, 
for reflections produced by a gradient in electron density 
as shown in Figure 4.14, for a frequency of 2.66 MHz, an 
electron density changing from 2.5 x 10^ nr3, and a col- 
lision frequency of 2.15 x 10^ s"-*-. 




Rx | / 1 Ro| 
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ELECTRON DENSITY (m~3) 


Electron densities calculated using all samples received 
(a) and rejecting all samples above 300 mV at the output 
of the receiver (b) . 
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shows an electron-density profile obtained at Wallops Island, Virginia 
using all the samples received, with the receiver gain adjusted to avoid 
saturation of any sample. Figure 4„17b shows the electron-density profile 
calculated from the same data, but rejecting all samples that reached a 
value greater than 300 mV at the output of the receiver, in a simulation 
of a receiver with a saturation level of 300 mV, In the processing of 
the data, if one ordinary sample was rejected because of excessive ampli- 
tude, the corresponding extraordinary sample was also rejected, to avoid 
an additional bias on averaging the data. The number of samples rejected 
at each height is shown on Table 4,1. As it is observed the rejection of 
some samples can introduce serious errors on the calculated electron 
density. 

At 76.5 km, for example, only one sample was rejected. The electron 

8 - 3 

density between 75 and 76.5 km, however, changed from 2.55 x 10 m to 
8 —3 

10 m” . At 81 km, 9 samples were rejected. The electron density between 

8 -3 8 —3 

79.5 and 81 km changed from 7.5 x 10 m to 5.7 x 10 m 

Reflections produced by gradients in electron density could explain 
some of the preferred heights of reflection observed in partial-reflection 
measurements by some experimenters [Gregory, 1961]. 

Preferred heights of reflection have been confirmed on the Urbana 
measurements, as seen in the histograms of occurrence of maxima in the 
A q profiles shown in Figures 4.18 and 4.19, and that correspond to the 
electron-density profiles of Figures 4.1 and 4.3. It is observed that 
heights of occurrence of a maximum of reflections in the histograms are 
associated to valleys on the electron-density profiles. Such association 
is one more indication that the valleys are not real, but produced by varia- 
tions in the scattering cross section with height. 
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TABLE 4.1 

Number of samples rejected at each height, in the calculation 
of the electron-density profile of Figure 4.17b. 


HEIGHT 

(km) 

NUMBER OF SAMPLES 
REJECTED 

67.5 

1 

69.0 

0 

70.5 

0 

72.0 

0 

73.5 

0 

75.0 

0 

76.5 

1 

78.0 

0 

79.5 

0 

81.0 

9 
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NUMBER OF MAXIMUMS IN THE A 0 PROFILE 

Figure 4.18 Histogram of occurrence of maximums in the A 0 profiles for the 
data used in the calculation of the electron density profiles 
of Figure 4.1. It is observed that the valleys in electron 
density of Figure 4.1, near 71 and 80 km, correspond to an 
increase in the number of maximums in A a at the same heights. 
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NUMBER OF MAXIMUMS IN THE A 0 PROFILE 


Figure 4.19 Histogram of occurrence of maximums in the A Q profiles, 
for the data used in the calculation of the electron - 
density profile of Figure 4.3. Valleys in the electron - 
density profile of Figure 4.3, near 80 km corresponds 
to an increase in the number of maximums in A a at the 
same height. 
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4 ,, 2. 2 Reflections produced by gradients in electron density; - ar 
approximate solution 0 The full wave solution presented in Section 4.2.1 is 
possible only for particular geometries of the electron-density profile. 

As partial reflections are very weak, however, a first order approximation 
can be applied to the determination of the reflection coefficients, and 
results obtained for an arbitrary geometry. The solution of the problem is 
presented by Collin [1966] for the case of an isotropic medium, and can be 
extended to a horizontally stratified magnetoionic medium, for quasi- 
longitudinal propagation, as will be shown below. 

The electron-density profile can be approximated by infinitesimal 
steps dN over distances dhj as shown in Figure 4.20, 

Each discontinuity dN will produce a reflection coefficient given by 


dR 


dn 


o, x 2 n 


<?,ar 


o.x 


(4.29) 


If second order reflections are not considered, the total reflection 
coefficient is given by: 


R 


o, x 


-2 ik n h d ,, „ ' 

e c o.x gjjCln n 0 ' X )dh 


'0 


(4.30) 


where L is the length of the reflector. 

If the height dependence of n is known, R can be readily evaluated 

Of CC Of 

from integration of equation (4.30) . 

As examples, the reflection coefficient will be calculated for three 
different geometries. The indices o and ar will be dropped from the equa- 
tions, and the results can be applied to the ordinary or extraordinary 



HEIGHT, h 



REFRACTIVE INDEX, n 


Figure 4.20 Approximation of the electron-density profile by 
infinitesimal changes dN over height interval 
dh, as used on the calculations of Section 4.2.2. 
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modes, depending on the expression used for the refractive index n. 

1) The electron density is an exponential function of height between 
heights zero and and constant outside this range. 

In this case the refractive index is given by: 


n 


n Q exp \ 


ln(n 1 /n Q ) 




(4.31) 


where n q and n ^ are the refractive indices at heights zero and 
respectively. Substituting equation (4.31) into (4.30), it results: 



expC-in^fej) In 



sin (^q^j) 

"Wi 


(4.32) 


where n was assumed as a constant equal to n Q in the exponential factor, in 
performing the integration of equation (4.30), 

A normalized plot of |i?| as a function of the height interval h ^ is 
shown in Figure 4.21. The results are similar to that obtained on 
item 3.2.1. 

2) The refractive index has the following height dependence: 



(4.33a) 



(4.33b) 





115 


Substituting equations (4.33a) and (4.33b) into equation (4.30), and 
performing a straightforward integration it results: 


E = j exp (-ing/Cg^) In 


f ' > 


'sin^fep^) 

V J 


Voh J 


(4.34) 


A normalized plot of |i?| as a function of the height interval h ^ is 
shorn in Figure 4,22, 

Observation of Figures 4,21 and 4.22 shows that the rate of change of 

|/?| as a function of the height interval h ^ depends on the geometry of the 

irregularity. As a consequence the variations of | R /R | as a function of 

cc c? 

h j will also depend on this geometry. 

3) The refractive index changes between heights zero and ir/y as a 
function: 


n - exp [a cos(yli)] 


(4.35) 


and is constant outside this height interval. This profile is characterized 
by having no discontinuity in the refractive index or in the first derivative 
of the refractive index with height. 

If equation (4.34) is substituted into equation (4.30) one gets: 


a exp 

r i} W 1 

y 


S o 

1 


2k o n o 

Y J 

2 ^ | 

- 1 



(4.36) 


1 



2 IRI/In (ni/no) 
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A normalized plot of E as a function of n/y is shown in Figure 4.23. 
The results show that resonances appear even in the case where no discon- 
tinuity in the refractive index or in the first derivative of the refrac- 
tive index are observed. 

In all the three cases shown above the amplitude of the reflection co- 
efficient is proportional to ln(n^/ng) , n ^ and being the refractive 
indices at the bottom and at the top of the irregularity, respectively. 

Rocket electron-density profiles near 65 km show electron- density 
gradients where In (n 1 /n Q ) is of approximately one order of magnitude 
below In observed near 80 km. Partial reflections at 65 km are one 
order of magnitude below partial reflectors near 80 km. 

Such reflections , in this case, can also be explained by gradients in 
electron density, 

Belrose ccnd Burke [1964] considered the electron-density profiles as 
proportional to the A q profiles at low altitude, assuming that the reflec- 
tion coefficient is proportional to the electron density. It should be 
noted, however, that the A ^ profile is the result of integration of 
reflections over a height range equal to one-half of the pulse width, and 
is not proportional to the scattering cross-section profile, being a func- 
tion of the pulse width. Even if the reflection coefficient is proportional 
to the electron density, the gradient of electron density with height should 

be sharper than the corresponding gradient in the A profile. 

o r 

4,3 Scattering from Random Irregularities in a Locally Homogeneous 

Background Medium 

Conti ibutions to partial reflections of reflections produced by random 
irregularities in electron density near 80 km altitude will be discussed in 
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this section. The background medium is characterized by a sharp gradient in 
electron density, and so it departs considerably from a homogeneous medium. 

To analyze such reflections, the theory of scattering will be extended to 
the case of a locally homogeneous background medium. The method of analysis 
will be similar to that used by Tatarski [1961] and applied to partial re- 
flections by Flood [1968], 

4.3.1 Scattering cross section of random irregularities in a locally 
homogeneous background medium. The model of reflector to be considered is 
characterized by weak irregularities superposed to a horizontally strati- 
fied, slowly varying, background medium. The medium will be assumed 
initially as isotropic, and the scattering cross section will be determined. 
The refractive index of the medium can be written as: 

— ' -> -V 

n(h ) = n(h) + n^(l?) , if f? is inside a volume V 

n(h) = n{ 0) , for points outside a volume V 

where n(h ) is the average value of n(h) , and n^(h) « nQi) , for points 
inside V, that will be assumed as having a lower boundary at the plane 
h = 0. 

Taking a time variation given by the factor exp (-iwt) the Maxwell 

\ 

equations have the following form, in a medium without currents or charges: 

v x t = imtl (4.37a) 

V x H = -iue Q n (4.37b) 

V • n E = 0, or ^ t = -2 1 * 7 (In n ) (4.37c) 
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The wave equation can be written as 

V 2 £ + k Q 2 n 2 E + 2 V[E ♦ V(ln w)]*0 


(4.38) 


If a series solution of the type E = E^ + is assumed, 
and a first order approximation used (Born's approximation) the electric 
field is given by 


= $ + 

o 


(4.39) 


Substituting equation (4.39) into equation (4.38) and collecting terms 
of the same order of magnitude it results: 

V 2 ! Q + k 2 n{h) 2 0 (4.40a) 


V 2 ! x + k 2 n(h) 2 * -2 k 2 Q nQOn^R) - 2V[£ Q • V In (n(h))] (4.40b) 


Equation (4.40a) has the well-known W.K.B. solution 


E q - 2 0 n{h)~ 1 ^ 2 exp - 

O 

h 

n[h)dh. 



0 


(4.41) 


Where the incident wave is assumed propagating vertically, and consequently 
the vector A is on the horizontal plane. It should be noted that the use 
of the WKB solution implies in not considering reflections produced by the 
background medium. Such reflections, however, werfe discussed in 


Section 4.2. 
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Substituting the value of ^ given by equation (4.41) into equation 
(4b 40b) the following equation for the scattered field is obtained: 


> 2 t 


+ k 0 n(h) 2 E. = 


2* 2 n(fc) 1/2 * (7?) exp 


f 

rh - 

i k Q 

n[h)dH, 

[ J 

0 


2V . 

[n Vi) exp[i 

■ h ' 

n(h)dh] • v [ln(n(h))]< 


l J 

0 J 


(4.42) 


To solve equation (4.42) the Green's' function will be determined. The 
Green's equation is given by: 

V 2 <7$, j?') + k 2 n{h ) 2 G(R,R') = -4ir6 (#-£*) (4.43) 

The solution of equation (4.43) is given by Liu [1967] and has the 
following expression: 


£(/?> ) = i. 


’ 2o 

2nr sin0 


exp 


*t- TT 

4 J 


exp J i k^rsine sin6 Q - u)(0j] 


/sin8 cos 6 
[q{h) q{h' )] 


1/2 


d 6 


for - « < h < h' 


(4.44) 


! 2 — 2 2 

where r = ]7?-7?' |, q (h) ~ n{h ) - sin 6, 8^ is shown in Figure 4.24, the 


path of integration r is shown in Figure 4.25, and 


w(8) = 


*h 

[n 2 (t) - sin 2 6]^ 2 dx 
**7i' 


(4.45) 



h (HEIGHT) 



(OBSERVATION POINT) 


Figure 4.24 Geometry of the source and observation points 
used in the determination of Green's 
function, equation (4.46). 
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Integration of equation (4.44) can be performed by the saddle-point 
method. The integral can be put in the form: 


G(J?,*») = i?(6) expj>/(e)]de 


Where 


, . . / K 0 f -in ) / sine cosQ 

F 6 " / 2?rr sin6 0 ® Xp ( 4 j q(?t’)] 1/2 


A ~ i T<qv 


/(0) = sine. sin0 o - 


The saddle point is given by: 


%P\e e = = ° 


Taking n(h ) = n(o) + y(h)> where y{K) « n(o). 

It results from equations (4,47c) and (4. SO) 


f ce fl ) 


cos9 sinfi,. 
s 0 


sine cosQ . ysinO cos6 dh 

s s a 1 e s 

r -, rt .2 . 2 . 1/2 C 0 + 2 r . r -, n .2 • 2 a i 3 / 2 

[re(0) -sin 0 ] J [ref 0) -sin 6 ] 

s h S 


(4.49) 


or, as y « 1 
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sin6 cose 

Q S 

COS0 sine n = s ry 7y 

S 0 r~~~ 2 . t A -i 1/2 

[«C0) -Sin 8 0 ] 


cos6 0 = 0 


(4.50) 


X 


X 


Resxilting 


sine = n(0) sine 
s u 


(4.51) 


Performing the saddle-point integration it results 



( 

f > 

■h 

expj 

I* k 0 

[ n(h ) - n( 0) ]dh ' 


l 1 

0 


(4.52) 


Using the Green's function (4.52) the solution for the scattered 
field 2^, in equation (4.41) is given by 


2~ 


-2fc Q n(0) exp (i &Q«(0)i?)4 0 


-t t 


^1 ” 4ttS 


exp{-i J^n(0)m*fl 7 }exp {£ 


rh 


* 0 


n.(fl') n{ h ') 1/2 exp J 

V l 


[n{K) - «(0) ]dh} dV' 


i k 


n ih)dn\ 


n(0) exp(i fe 0 h(0)i?) 


2vi? 


V { nOi ) 1/2 exp[£ k Q 




V 


nWdh}A Q . v ln(n 1 (R''))) 


exp{--£ k^rri ■£'} . exp{i k Q 


J 0 


\nW) - n(0)]dh}dv • 


(4.53) 


where the following approximation was introduced: 

exp(t k 0 2 ?) = exp (2 R) exp (-i m ♦ i? 1 ) 
m being a unit vector in the direction of j$„ 


(4.54) 
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The second term in the second member of equation (4.53), in a first 
order approximation, produces a component of the electric field in the 
direction of r. Such longitudinal component is annulled by the longi- 
tudinal part of the first term of the second member of equation (4.53) 
and does not contribute to the scattered power. This statement is proven 
in Appendix I. As a result it can be written 


^ -2&q n( 0) exp(i k Q n(Q)R ) 

£>1 = J 


(4.55) 


where 


I * 


r h 


n.(^’) n( h') I//2 exp { i k Q n(h)dh} exp { -i & 0 n(0)m * $'} 

v Jo 


exp { i feg 


[n(K) - n{Q)]dh } dV' 


0 


The field can be determined from (4.37a) 


4 - 1 -¥ 

ff, = -A- V x E. , or 

1 t-cuy 1 


_ 2k 2 n{0) rA _ ^ 

H 1 = f 7r4coy V X { iT eX ? {i k 0 M(0)i?}J / 


(4.56) 


2k 3 Q n( 0) 2 


exp(i /c 0 n(0)i?} I • (m x ? Q ) 


4Til?(oy 


(4.57) 
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The Poynting vector is given by 


$- 


F x ft* 
J 1 X a l 


2 k s 0 n( 0) 3 
(4ir) 2 ojuff 


jl J* [1 Q • ? 0 -(M 0 )]m 


and the flux of power density in the direction of m is given by 

* * 2k l J C0) 3 2 2 

S = 5 • m - — y 4 sin x .11* 

(4 it) covd?“ 

—y 

where X is the angle between the direction of the incident field, E 
the direction of observation, r. 

The average power density in the direction of m is given by 


S m " 

m 


2k S Q n( 0) 3 
(4 it) 2 . 


sin x IX* 


where the bar means average. 

For the scattering cross section, a, defined by 

Scattered power per unit solid angle, per unit volume 
Incident power per unit area 

for following expression results 


4kg n(0) 3 

(4 it) 2 


II* . 2 
* -y~ sm X 


i 


(4. 58) 

(4.59) 
, and 

(4.60) 


(4.61) 
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This product II* is given by 


fP= I n(h ) 1/2 n*(h 2 ) 1/2 expi-i ^*(0)^} 

h, h* 


1 2 


exp{i &qH( 0)7?^} exp{2 f ?Cq 


— . -¥ -V 


[n(h 1 ) - nCO)]^} • exp{r ^n*(0)m*i? 2 } 


i^i 


exp{-i k 0 n*( 0)h 2 }exp{-2i fc 0 


[n(ft ? ) - n* (0)3 ^2 -^1^2 (4.62) 


(tf^n* (i? 9 ) can be identified as the correlation function of the refractive 
index, B^^R^) . 

For a general medium equation (4„62) cannot be integrated, unless the 
analytical expression of B^iR^,^^ i s known. 

A slowly varying medium, however, can be assumed as a locally homo- 
geneous random medium [Tatarski, 1961], and £^(i?^,i? 2 ) can written as 

vM 2 > ■ ° 2 

Substituting equation (4,63) into equation (4,62), and introducing the new 
coordinates 


R l +R 2 


n 




(4.63) 


2 


(4.64a) 
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■+ 




(4.64b) 


It results 


II* = 


V, 


h 1/2 


h 1/2 


~lh n *(h s ~ 2^") d (# 5 ) ex P^ [ fle ( n (°)))} 


S e 


exp{-2?c 0 (JOT(n(0)))^*? s }exp{2'i^ { 


VV 2 

o e 


[n(fr) - n(0)]c#z) 


exp{-2i^Q 


•h -h /2 
S e 


[n*(h) - n*(0)]dfo} exp{i k n [2?e(n(0)))fr 

v 6 


e*p{2k[I)«(n(Qj)]fc~}»b (fl )d7„ dF 

v H 62 o g 


(4.65) 


where JcgS = (£q - kgin) , Re and In mean real and imaginary parts, respec- 
itvely. For a homogeneous random medium without losses, n(h ) is real and 
constant, and II* reduces to 


II* = n 


b [R ) exp{£ k^t * n(0)} dV <- 


dV 


J v q h 

S e 


* n V $ $r>n - k n m) 
n u 


(4.66) 


where $^(^Qn - k^n in) is the spectral density of power. 

In this case the scattering cross section reduces to the well-known 
expression: 


a = 


4 -4 

0 n .2 


4 k n n 


(40 • 


sin x $ n C^o n " k Q n ^ 


(4.67) 
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If the background medium is not homogeneous, equation (4.65) can be 
solved if n(h) is known. Assuming a linear variation with height: 


n(h) = 1 - a(h - h Q ) , h >_ 0 


(4.68a) 


n(h ) = 1 + a h 


0 


, h < 0 


(4.68b) 


it results 


ZP = f o 2 (R s ) exp 2 k Q (Im) Qi 2 - $ • } • 

J V 

s - 

h 1/2 __ h 1/2 

n(h Q + n* (h„- ~) b (R ) exp{?: k \b-t i?e(n(0))] 


V 


2(Rea)h ?»„]} exp { (Zfia) —}dV dV 

& O 


e S 


(4.69) 


1/9 ^ \ / 2 • ^ ' p 

The factors n{y^- + h ) , n* (h - ~) and{exp (Ima)—] can be con- 

sidered as a constant in the integration in V^, since b (R ) goes to zero 

& Tl & 

very rapidly. Recognizing j exp {i^.QR}b^(R)dV = 
the power spectrum of rZj(/?), it results: 


II* = 


° 2 e ^P (ft^ 2 


”^Vo [fie n (°)^ 


3/ 


- ^ ^ 0 )} n(^ 5 ) i/2 . n*(hg ) 1/2 
, fe Q [J?e n(0)]b g - 2 0?ga)y }dy g 


(4.70) 
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For backscattering: 


U * = ) v ° Cffg) exp{4fe 0 j [Im n{h)]dh} *^{ 2 ^ [2te k(* 5 )]> 


n(h s ) l/2 n*(h s ) 1/2 dV s 


(4.71) 


n(0) 3 


i »cy 1/2 '>*(y 1/2 <- 2 «y 


exp {4 Jc q [/m(w(fc))]dfc} 4^(2^ [Re n(h )]} dV^ 


(4.72) 


4„3.2 Application to partial reflections. For quasi-longitudinal 
propagation, for the ordinary and extraordinary modes of propagation. 


V> 4 C 3/2 (y o,x ) + *"4 C 5/2^o > x ) ^ 


(4.73) 


Where $ . is the power spectrum of the electron-density fluctuations. 


U> ± 

and y - — . From equation (4 0 72) and (4.73) it results for re- 

o 1 CC V 

flections produced by a pulse of width W 3 centered at a height h : 


k t r 2 


— exp{4lc 0 Im (.n 0iX mw 

lV '\-CW/A 


th^CW/A 


^ N {2k^[Re{n (hmUV 


(4.74) 
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The signals A and 4 measured at ground are given by: 
0 cc 


A = a V exp {2 
o,x o,i c 


h -CP// 4 
o 


kq [Im n(h)]dh} 


(4.75) 


The ratio {A /A ) is given by: 

£u O 


Where 



2 



rh -CP// 4 
JQ 


K 


XyO 


= V I?7! 


(4.76) 


(4.77) 


Using the expression of o given by equation (4.74) it results: 

o 


r ^ 



2 


fe* <1/2^ * r C 5/2^» 

[i'o C 3/2 ( ^ ) + r- C 5/2^>l 


; (h +CP//4 

exp (4 0 k {2k [Ren (fe)]}dV 

y ft -CP// 4 33 0 x 

/ ;_o 

C rh +CP//4 


exp (4 k dh'H J2k [Ren (h)]MF 

V h -cm a 0 “* 0 ° 

^ J o 


rh^-CW/ A 


• exp { 2 


<w^ } 


J 0 


(4.78) 


If the background electron density does not change inside the volume V, 
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equation (4.78) reduces to: 


ro 

X 


oJ 






22 ?5 2 * exp{4 

C 3/2<iV> * r *1/2^ 


h +CW/A 
o 

[K -K )dh } 
} h -CW/4 x ° 


( h Q -CW/A 


exp { 2 


iK x -K o )dh} 


(4.79) 


0 


That is identical to the results obtained by Flood [1968]. 

If the irregularities are uncorrelated, equation (4.78) takes the 

form 



exp 


tvl C 3/2fex> * r c ln^ 

bl C 3/ 2 <V * f- C 5/ 2 fe 0 » 
rh o +CW/ 4 





exp (4 

h 

exp (4 
J V 


'h +CW/4- 
j h°-CW / 4 

"rhTcw/4 

jh o -CW/4 


K dh } dV 
x 


K dh} dV 
o 


(4.80) 

; 


Such expression is similar to that developed by Cohen [1971]. 

As a numerical application, the ratios A I A were calculated at 78.5 

x cp 

5 - 1 

and 81.5 km for a pulse of 20 asec, a collision frequency of 7.7 x 10 s 

and the electron-density profile shown in Figure 4.26. The following power 

2 2 

spectrums were assumed: A Gaussian spectrum (<t>^ = exp(-a k /2), with 

“ft 

a = 0, 50 and 70 m) , and a power law spectrum ($^ = k n - 11/3 and 6) . 

With the values of A / A obtained, the electron density between 78.5 

x o J 

and 81.5 km was calculated using Belrose and Burke's [1964] theory. For 
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Figure 4.26 Electron-density profile used on the calcula- 
tions of Section 4.3. 



135 


comparison, the electron density was calculated assuming the existence of 
sharp reflectors at 78 „5 and 81.5 km, instead of a random medium. The 
results are shown in Table 4.2 e They show that if the power spectrum is a 
Gaussian spectrum, with the parameter as great or equal to 50 m, a serious 
error can be produced in the calculation of electron densities using 
Belrose and Burke's theory* If the power spectrum is a power law spectrum, 
such error will be small. 

4.3*3 Electron-density fluctuations necessary to produce partial 
reflections ♦ An estimate of the fluctuations of electron density necessary 
to produce partial reflections of amplitude comparable to the observed re- 
flections can be obtained from equation (4,67), for altitudes above 80 km, 
since for such heights the losses produce a small effect on the value of 
the scattering cross section for the ordinary mode, and in a first approxi- 
mation, the background medium can be considered as homogeneous. 

If the power spectrum is of the form 

_ 2 

<& n (£) = 2(ir ) 3/>2 L 3 • exp (- L h ’ n i (4.81) 

it results from equations (4„67) and (4.81), for the scattering cross 
section of the ordinary mode: 

.4 -4 -2 

4/C 0 n o , 0 v 3/2 r 3 r *0 1 2 r 2 n ~ ,, _ 

°o = 2~ * L ex P(" 2~ k 0 L * n lo (4.82) 

(4tt) 

Where the angle x in equation (4.67) was taken as 90° (back- scattering), 
and is given by: 



TABLE 4.2 


Electron-density calculations using Belrose and Burke's theory [1964], 
for reflections produced by sharp reflectors and by a random medium, 
for the electron-density profile of Figure 4.26, and a collision 
frequency of 7.7 x 10^ sec"-'-. 



Sharp 

Gaussian 

Gaussian 

Gaussian 

Power 

Power 





Law 

Law 


Reflector 

Spectrum 

Spectrum 

Spectrum 

Spectrum 

Spectrum 



a=0 

a=50m 

a=70m 

n-ll/ 3 

n~6 

ln[(/t x M 0 ) 2 ] 

IntCVVll 


-0.683 

-0.43 

-0.05 

-0.66 

-0.65 

Electron 
density at 
80 km (nr 3) 

i 

i 

l.SxlO 9 

l.lxlO 9 

0.69xl0 9 

• 

0.12xl0 9 

1.03xl0 9 

1.02xl0 9 
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where 



s M 2 

4 ( 1 +?) 2 


(4.83) 



e m<u 

o 


(4.84) 


ilfj being the r.m.s. electron-density fluctuation. 

If the transmitting antenna has a lobe with an angle 6, the trans- 
mitted pulse with W, and the irregularities are distributed over all 
the volume occupied by the pulse, the power back-scattered per unit angle 
solid per unit incident power density, from a height h is given approxi- 
mately by 


* 2 

o 


V 

h 2 °° 


Tftsinj) 



(4.85) 


If R is known, the value of can be determined from equations (4.82) 

and (4.85). A plot of is shown in Figure 4.27, for a pulse of 50 

9 -3 

psec, and antenna lobe of 15°, a background electron density of 10 m , 

a frequency of 2.66 MHz, and a gyrofrequency of 1.5 MHz. 

As will be shown in Chapter 5, the reflection coefficient for the 

ordinary mode, near 80 km, is of the order of 10 - ^, and the electron density 

9 3 

of the order of 10 m . From Figure 4,27 we conclude that, if L ^ 50 m, 
fluctuations of electron density of the order of 0.1 percent will be 
enough to produce such reflections. 
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4.4 The Nature of the Irregularities 

The choice of the best model of reflector for partial reflections is 
impaired by the lack of knowledge of the real structure of the irregulari- 
ties producing reflections, Belrose and Burke [1964] arrived at the con- 
clusion that the reflections are produced by isolated reflections and not 
by a turbulent medium. Such conclusion was based on the fact that 
reflections produced by pulses of different width showed maximums at the 
same height, and more defined minimums for the narrower pulses, and that 
reflections produced by pulses with different frequencies (40 kHz apart) 
showed similar structures. 

Further work, however, [Belrose* 1970] indicated that the reflection 
mechanism is probably produced by a combination of turbulent and isolated 
scatters, Gregory [1961] observed that the reflections are produced 
principally from the heights of 55, 61, 66, 74 and 86 km. Von Biel et al. 
[1970], based on the correlation coefficient between the ordinary and 
extraordinary reflections, suggested that the reflections are produced 
by small irregularities distributed over all the volume occupied by the 
pulse. 

Fraser and Vincent [1970] studied the irregularities of the D region 
by measuring the phase variation of the received signals and the space 
and time correlations of the ground diffraction pattern. Sixty percent 
of all the 70-80 km and 20 percent of the 80-90 km reflections were co- 
herent echoes. During the winter, reflections showing a pronounced 
stratification were observed near 85 km. 

Manson et al. [1969] based on small irregularities in electron - 
density profiles published by Mechtly and Smith [1968] concluded that the 
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isolated irregularities of the rocket profiles are sufficient to explain 
the partial reflections. It should be noted, howeveT, that the irregularis 
ties observed on rocket profiles are at least partially produced by the 
precession of the rocket, which changes the angle of attack of the current 
probe with the ionosphere, producing variations of current that are not 
really produced by irregularities. In the interpretation of the rocket 
data such irregularities cannot be taken into account. 

Von Biel [1971], measuring the amplitude distribution of the received 
signals observed that below 80 km the echo amplitude distributions are pre~ 
dominantly Rayleigh-like in character, suggesting that they are produced 
by a turbulent medium. 

Above 80 km the distribution approximates to a Rice distribution with 
the Rice parameter increasing with height, showing that above 80 km there, 
is a contribution of coherent scattering. Such conclusion supports the 
theory that above 80 km the scatters are at least in part produced by 
gradients in electron density. 

As can be concluded from the available information about the structure 
of the irregularities, that was summarized above, further work must be done 
on the subject. 

4.5 Electron-Density Calculations for a Region with Sharp Gradients in 

Electron Density 

In Sections 4.2 and 4.3 it was shown that the received signals A Q (h) 

and A (h) , in some circumstances, do not represent the scattering cross 

%c 

section at the height h. 

To overcome such a problem, at least two solutions can be tried: 

1) To use only maximums of the a and a profiles in the calcula- 

0 

tion of the electron densities. 
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2) To deconvolute the average A q and A profiles, 

(a is the signal recieved from one sample, and A „ the average of 

O yX OyX 

a o) Both solutions will be discussed below, 

O y CL 

4.5,1 Electron-density calculations using only maxima of the a Q and 
a profilesc If a maximum on the a arid a profiles at a height is con- 

X 0 X 1 

sidered as produced by a reflection from this height, and not the result 
of phase interference of waves reflected from different heights, the A 

O y\C 

profiles obtained by averaging only the signals measured at such maximums 
will be representative of the scattering cross-section profile. A measure- 
ment made by this method, besides being based on the above assumption has 
the inconvenience that only a few maximums are observed on each received 
sample, they are not equally distributed over all heights, and as a con- 
sequence a great number of samples has to be taken to obtain a reasonable 
number of samples at each height. The sampling and processing time in 
this case becomes very large. 

An electron-density profile obtained by the method above described, 
at the partial-reflection system of the University of Illinois is shown 
on Figure 4.28a. Table 4.3 shows the total number of samples and the 
number of maximums observed at each height. The electron-density profile 
obtained from the same samples, but averaging all the signals at 1.5 km 
of interval is shown in Figure 4.28b, for comparison. 

The criteria used to consider a signal at a height h as a maximum 
was the following: 

a 0j>x W > (4.86a) 


a (h) 


a {h- 1.5) 

o.z' 


(4.86b) 



HEIGHT (km) 



Figure 4.28 Electron-density profiles measured at Wallops Island, Virginia, on Feb-14-1972, 
at 12:12 local time, using only maximums of the A 0 and A x signals, (a), or 
averaging all the samples at 1.5 km of height interval, (b) . 
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TABLE 4.3 

Number of maximums observed at each height, in the data used in the 
calculation of the electron-density profiles of Figure 3.28. 


Height (km) 

Number of maximums 
in the A profiles 

Number of maximums 
in the A x profiles 

69 

145 

150 

70.5 

137 

127 

72 

185 

202 

73.5 

170 

190 

75 

170 

196 

76.5 

165 

187 

78 

132 

166 

79.5 

195 

149 

81 

62 

132 

82.5 

i 

191 

141 

84 

263 

182 

85.5 

247 

224 

Total number of 

samples 

1,584 
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Observation of Figure 4.28a and 4.28b shows that the processing method 
can alter significantly the resultant electron-density profile. The 
profile of Figure 4.28a, obtained using only maximums of the a Q and a ^ 
profiles shows a steep gradient above 84 km that is not observed on the 
profile of Figure 4.28b, obtained by the method described in Chapter 2. 

Below 75 km the profile of Figure 4.28a presents lower electron 
densities than the profile of Figure 4.28b. 

The signal processing method using only the maximums of the received 
signals, although time consuming, should be considered in partial- 
reflection systems having a transmitter with a rapid pulse repetition 
rate (order of 10 s ^) and a capability of processing the number of 
samples necessary to get an electron-density profile (order of 15,000) in 
a reasonable time. 

4.5.2 Deconvolution of the A and A -profiles . Each received 

O X 

sample a or a from a given height is given by the expression: 

O %C 


a 


0 3 X 


.00 


0 + OO 

s (t- 70 a(t)dr 

J _oo 


(4.87) 


Where a(t) is the waveform of the transmitted pulse and a(r) is the 
scattering cross section per unit volume. 

If integral equation (4.87) is solved, the scattering cross section 
as a function of height can be obtained. A solution to this equation, and 
consequently the deconvolution of the A and A profiles was given by 
Austin et al. [1969] and Coyne and Belrose [1973]. 

As the reflections from different heights are not correlated, equa- 
tion (4 0 87) can be averaged without taking the phase of g(t) into 


\1 
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consideration, resulting: 


A 


2 

o t x 


f+co 


O (t) s{x-h)dx 

Csj Hi 

J ^CQ 


(4.88) 


Equation (4.88) is a Fredholm integral equation of first kind, and 
admits oscillatory solutions „ Such solutions, although mathematically 
plausible, do not correspond to the physical solution of the problem. To 
avoid the oscillations it is necessary to know with great precision the 
variations of A and A with altitude, from the lowest to the highest 

O X 

level, and to determine precisely the waveform of the transmitted pulse. 

The solutions are particularly sensitive to the values of A (h) 

Of x 

at lower altitudes, and small variations of such values will produce 
solutions with oscillations growing with height* Due to limitations on 
the precision of the measurement of A Qi) at lower heights, where the 
signal-to-noise ratio is low, all the attempts made during the execution 
of the present work to deconvolute the A (i h ) profiles did not produce 

Q y X 

reliable results,. Due to the importance of this subject for partial re- 
flection, and the convenience of further work to improve the deconvolution 
technique, the method used and the results obtained will be discussed 
below. 

Methods of deconvolution of physical data have been discussed by 

several authors [Stone, 1962; Grisson et al, , 1968; Ritchie and Anderson , 

1966] , A review of such methods is made by Raveck [1969] . 

The technique used in decovoluting the A (h) profiles was that 

O 3 X 

given by Grisson et aZ , [1966] „ In this technique, the integral equation 
(4.88) is approximated by the following series: 
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Am 2 = l s i (/z)C i A (4.89) 

Where the subscripts o , x were dropped, the region of integration was sub- 
divided into n sub-regions of width A - h 2 - and the coefficients Cb 
are chosen to yield an approximation of the integral given by Simpson’s 
rule, o(h) has been replaced by a set of expansion coefficients 
To determine o . an error function is defined by 

•t 

E = S x + yS 2 (4.90) 
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S 2 is a positive definite quantity whose magnitude is a direct measure 
of the continuity of the solution {o.K The solution is obtained by mini- 

*ls 

mizing the error E, for a given choice of y, and given by the simultaneous 
system of equations 


a e 

3o . 

-I 


= 0 


i - i j t , . . . n 


(4.94) 


By adjusting the value of y one imposes a continuity requirement on the 

solution and reduces the fluctuations to a certain degree. 

The computer program to deconvolute the A q h ) profile by using the 

above method is in Appendix II (DECOMM) . The results obtained were not 

reliable, and even adjusting y over a wide range, the oscillations appear 

on the solutions, and as the oscillations on the A ( h ) and A (h) profiles 

are not in phase, they produce large fluctuations in the A / A ratio. In 

Figure 4.29 it is shown an electron-density profile obtained by the method 

described above (Figure 4.28a) and by using the program PROAX (Figure 4,29b). 

4.5.3 Correction for the different group velocities of the modes of 

propagation . The different group velocities for the extraordinary and 

ordinary modes of propagation can be an additional source of error in the 

determination of electron-density profiles above 80 km, since in the signal 

processing of partial reflections A (h~) and A {h) are obtained by assuming 

that both signals propagate with the velocity of light in free space up to 

the height h. In Figure 4.30 the group velocities v and v for both 

go gx 

inodes of propagation as a function of electron density are shown. It is 

9 -3 

observed that for N > 10 m the difference between v and v has to be 

- go gx 

taken into consideration* If the wave propagates 1.5 km in a region where 


j 



HEIGHT (km) 



ELECTRON DENSITY (rrr 3 ) 

Figure 4,29 Electron-density profiles obtained by the programs PROAX (a) and DECOMM (b) measured at 
Urbana, Illinois, on October 19, 1971, at 11:00 A.M. local time. 
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N = 5 x 10 9 m“ 3 , for example, the error in the determination of A^(fr) will 

be of 180 m. If the A (h) profile changes considerably over this distance, 

x 

the corresponding A /A profile will be in error. 

X O 

In an attempt to compensate this difference in propagation velocities 
computer program "RETARD” was developed. In this program, that is listed 
in Appendix II the corrected value of at a height h, is taken as 


A (fc+1.5) - A 00 

*_,(*) = + V (?0 3 

xp x gx (1.5 x 10/o ) 


(4.95) 


where o q is the velocity of light in free space. The corrected value of 
A Qi) is obtained by making a linear interpolation between two measured 

X 

values, at heights separated by 1.5 km. 

To calculate V (h) it is necessary to know the electron density at 

gx 

this height, and to determine the electron density it is necessary to 

know A (fc). This problem is solved on the program RETARD by an iteration 

xxx 

technique. The electron density is initially calculated by taking 
A xr Qi) = A^{h')<> With this value of electron density, A^Qi) is calcu- 
lated using equation (4.95). The value of electron density is recalcu- 
lated, and a new A (h) determined, until the value of electron density 
xx 3 

before and after a correction of A (h) shows a difference of less than 

XX 

15 percent, that is considered as satisfactory in the present calculations. 
Electron-density calculations with and without correction for the different 
group velocities are shown in Figure 4.31. The profile of Figure 4.31a 
was calculated using the program PR0AX, described in Chapter 2. The 
profile of Figure 4.31b was calculated using the program RETARD. 



(km) 
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5. PARTIAL-REFLECTION MEASUREMENTS DURING THE 1971-1972 WINTER 

During the 1971-1972 winter, a coordinated rocket program, headed by 
the University of Illinois, was established to study the winter varia- 
bility of the D region through measurements of electron density, ion 
composition, pressure and temperature, at Wallops Island, Virginia. It 
was decided to make measurements on three different days, with the fol- 
lowing characteristics: 

a) a day of low ionization of the D region; 

b) a day of high ionization, with no magnetic disturbances during 
at least 8 days prior to the measurement; 

c) a day of high ionization, during magnetic disturbances. 

The rocket measurements were performed only at the low ionization 
day, on January 31, 1972, since no day attending the specifications of the 
items b and c above were observed during January and February of 1972. 

To determine the conditions required for a given day, several measure- 
ments were made every day by different workers in United States and Canada, 
and the results sent to Wallops Island. 

The measurements performed were the following: vertical absorption at 

1.8 and 3.3 Mhz, oblique absorption at 2.2 Mhz, VLF phase, AjA Q ratio at 
76 km, at Ottawa, Canada and Raleigh, North Carolina, 10 mb temperature, 

A index, particle precipitation, solar radiation fluxes, and electron- 
density profiles by partial reflection at Wallops Island, 

In this chapter will be presented the results of the partial-reflection 
measurements made at Wallops Island, by the University of Illinois. 

The equipment used is described in Chapter 3. Electron densities were 
normally measured from 9 to 13 o'clock, local time, with exception of a few 
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days when interference problems or hardware defects prevented the measure- 
ments. On the average, 10 profiles were taken each day. The profiles were 
obtained from 513 frames, each one corresponding to one ordinary and one 
extraordinary sample. The pulse repetition rate used was 2 double pulses 
per second, and consequently, each electron density corresponds to the 
average over approximately 4 minutes. When strong interference was observed 
during a measurement, it was interrupted and started again when the inter- 
ference level reached a reasonable level. As a result, some profiles were 
obtained over a period of time greater than 4 minutes. The pulse width 
used was of 50 ysec. 

Results obtained during the month of October of 1971, at Urbana, 
Illinois, using the partial-reflection system of the University of Illinois, 
that is described in Chapter 3, are also presented in this chapter. 

5.1 Partial Reflections in October , 1971 

The electron densities at 70.5, 75 and 78 km, measured at Urbana, 
Illinois, for the month of October, 1971, are shorn in Figure 5.1. The 
electron densities presented are the median electron densities obtained 
from measurements performed at solar zenith angles between 65 and 55°. 

The electron density at each height is the result of averaging electron 
densities measured at two adjacent heights, separated by 1.5 km. The 
A / A ratios, at 69, 76.5 and 82.5 km, for a solar zenith angle of approxi- 
mately 60° are shown in Figure 5.2. 

As can be verified from Figure 5.1, some variability, although not 
strong, was observed in this month. The ionization on October 14, for 
example, was approximately 50 percent greater than on October 13. 

On some days, measurements were made during all the morning, including 
sunrise time. One of such measurements, made on October 14, is shown in 
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Figure 5.1 Daily median electron densities for the month of October, 1971 at 
70, 75, and 78 km of altitude, solar zenith angles between 65 and 
50 degrees, measured at Urbana, Illinois. 
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Figure 5„3» For solar zenith angles above 80° no reflections above the 
noise level were observed between 65 and 80 km, but strong reflections were 
measured above this height., Reflections between 65 and 75 km were observed 
only for solar zenith angles below 70°. During sunrise, reflections pro- 
duced by the C layer, at heights between 60 and 64,5 km were observed. One 
example of such reflections is shown in Figure 5.4, where the A profile 
measured on October 14 is plotted, for a solar zenith angle of 87°. Re- 
flections above the noise level are observed at 61.5, 63 and 64.5 km, and 
above 81 km. 

Reflections from the C layer during sunrise is a good way of measuring 

the collision frequency at lower altitudes, since at that altitude the 

ionospheric absorption is very small, one can take A /A = R t /R , and from 

o cc o 

this ratio calculate the collision frequency. For the data obtained on 

October 14, collision frequency of 2.0 x 10 7 was calculated at 63 km. 

5.2 Partial Reflections in December , 1971 

The median electron densities at 70, 75 and 78 km for the month of 

December, 1971 are shown in Figure 5,5. The median was obtained for 

solar zenith angles between 70 and 60°. The A l A ratios, at 69, 76 and 

x o 

82.5 km, for a solar zenith angle of approximately 60° are shown in 
Figure 5.6. 

Two periods of low absorption, on December 8 to 10 and December 15 
to 17, and two periods of high absorption, on December 12 to 13 and 
December 18 to 19 were observed. The variation from low to high absorp- 
tion showed a periodic behavior, with a period of approximately 6 days. 
Figure 5.7 shows the median electron-density profiles during a day of high 
absorption, December 12, and a day of low absorption, December 9. It can 
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Figure 5.4 A Q profile measured on October 14, 1971, at a solar 
zenith angle of 87°. 
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DAY (DECEMBER, 1971) 

Figure 5,5 Daily median electron densities, for the month of December, 
1971, at 70, 75 and 78 km altitude, solar zenith angles 
between 70 and 60°, measured at Wallops Island, Virginia. 



a - 69.0 km 
b - 76.5 km 
c - 82.5 km 
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DAY (DECEMBER, 1971) 

Figure 5.6 A x /A 0 ratios for the month of December, 1971, at 
69, 76.5 and 82.5 km altitude, solar zenith 
angle of approximately 60°. 






Figure 5.7 Median electron densities for a day of low absorption (December 9) and a day of 
high absorption (December 12) . 
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be observed that the high absorption day is characterized by an increase 
of ionization all over the D region, from 69 to 76 km. Above this height 
no measurement was made on December 12, due to saturation of the receiver. 
At 74 km, the electron density on December 12 was approximately four times 
greater than on December 9. 

The daily variation of the ionization was closely correlated to the 
variations of geomagnetic activity, as can be observed on Table 5.1, that 
gives the A p indices. On December 12, 13, 18 and 19, days of high ioniza- 
tion, the A index was of 9, 13, 22 and 11, respectively. On December 8, 

P 

9, 10, 15, 16 and 17, days of low absorption, the A index was 3, 7, 3, 3, 
8 and 67, respectively. December 17, although being a day of low absorp- 
tion and high magnetic activity, was prior to a day of high ionization. 

Such correlation suggests that the variations in ionization during 
December 71 were produced by geomagnetic effects. 

The A I A profile of Figure 5.6 shows that the best indicator of the 
ionization is the A l A ratio at 76.5 km. Above this height, the A JA 
ratio is influenced principally by the strong gradient in electron density 
that appears near 80 km, and at lower altitudes, as 69 km, the differen- 
tial absorption is too small, producing practically no variation in the 

A I A ratio. 
x o 

5.3 Partial Reflections in January and February , 1972 

The median electron densities at 70, 75 and 78 km, for solar zenith 

angles between 70 and 60°, for the months of January and February of 1972, 

are shown in Figures 5.8 and 5.9. The A /A ratios at 69, 76.5 and 

X o 

82.5 km, solar zenith angle of approximately 60° are shown in Figures 5.10 
and 5.11. 
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DAY (JANUARY, 1972) 


Figure S.8 Daily median electron densities at 70, 75 and 78 km 
altitude, solar zenith angle between 70 and 60°, for 
the month of January, 1972. 
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DAY ( FEBRUARY, 1972 ) 

Figure 5.9 Daily median electron densities at 70, 75 and 78 km 
altitude, solar zenith angle between 70 and 60°, 
for the month of February, 1972. 
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Figure 5.10 


A /A ratios for the month 
7o,5 and 82.5 km altitude, 
approximately 60°. 


of January, 1972, at 69, 
solar zenith angle of 
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DAY (FEBRUARY, 1972) 


Figure 5.11 A X /A Q ratios for the month of February, 1972, at 69, 
76.5 and 82.5 km altitude, solar zenith angle of 
approximately 60°. 
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The behavior of the D region during the months of January and 
February of 1972 was quite different from that observed in December, 1971. 
Only one day of high absorption, on January 11, was observed. Low absorp- 
tion days appeared more frequently, as on January 25, 28, 31 and 
February 10. The daily variations of electron densities were much smaller 
than in December. 

The days of higher ionization on January 11, 15, 26 and 29, and on 
February 14, 20 and 26, were also associated to variations of geomagnetic 
activity, as can be observed in Table 5.1. 

Such results indicate that variations of ionization during all the 
winter 1971-1972 were mainly due to geomagnetic effects. 

Measurements of electron densities during the occurrence of solar 
flares were made several times. The flare that produced the most pro- 
nounced effects on the D region occurred on January 19, between 11:34 and 
12:45 hours, local time. The electron-density profiles measured during 
this period are shown in Figure 5.12. The profile at 11:09 was taken 
before the flare. The ionization reached its maximum values on the pro- 
file measured at 11:47. All the D region, down to 63 km, was affected. 

The electron density produced between 63 and 64 km by the X-ray radiation 

8—3 

was of approximately 3 x 10 m . The X-ray solar fluxes during the flare 

are shown in Table 5.2. 1.8 vertical absorption measurements made at 

Greenbelt, Maryland [£. Gnanalingam, private communication] showed a good 

correlation to the AJA q ratio at 76.5 km, measured at Wallops Island as 

can be observed on Figures 5.13 and 5.14, that show the scatter plot of 

the A l A ratio at 76.5 km versus 1.8 MHz vertical absorption, for the 
0 

months of January and February, 1972, respectively. 



TABLE 5.2 


X-ray solar fluxes during the solar flare on January 19 
1972 between 11:34 and 12:45 hours, local time. 
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ELECTRON DENSITY (m -3 ) 

Figure 5.12 Electron-density profiles measured before and during the solar flare that 
occurred on January 19, 1972, between 11:34 and 12:45 hours, local time. 
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1.8 MHz ABSORPTION (dB) 


Figure 5,13 Scatter plot of A x /A 0 ratio at 76,5 km, measured 

at Wallops Island, Virginia versus 1.8 MHz vertical 
absorption measured at Greenbelt, Maryland, for 
the month of January, 1972. 
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1.8 MHz ABSORPTION (dB) 

Figure 5.14 Scatter plot of A x /A ratio at 76.5 km altitude, 

measured at Wallops Island, Virginia, versus 1.8 MHz 
vertical absorption measured at Greenbelt, Maryland, 
for the month of February, 1972. 
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Figure 5.15 shows the electron-density profiles measured on 
January 31, 1972, between 11:40 and 12:30 hours, local time. This day was 
chosen as the low day, and a Nike Apache rocket launched at 12:30 for the 
measurement of electron densities and collision frequencies, by Faraday 
rotation, differential absorption and current probe. The technique 
employed on the rocket measurement is the same as described by Meohtly 
et ale [1967] o The results of the rocket measurement are also shown, on 
the traced line, in Figure 5.15 [S'. A. Meohtly , private communication]. 

As can be observed on Figure 5.15, the profile numbers 1 and 2, taken at 
11:40 and 11:58 hours are lower than the profile numbers 3 and 4 taken 
at 12:12 and 12:30 hours. This difference is the consequence of a small 
flare, that occurred after 12:00 and was strong enough to change the 
ionization of the D region. Comparison of the rocket profile with partial- 
reflection profile number 4 shows good agreement above 79 km, with a dif- 
ference in height between the two profiles of approximately 1.5 km, which 
is less than the precision of the partial-reflection measurements. 

Below 72 km, the profiles will also be in good agreement if the partial- 
reflection profile is lowered by 3 km. Between 72 and 79 km, however, there 
is a complete disagreement. There is no explanation for such disagreement, 
unless the rocket profile suffered the influence of a sudden change in the 
solar X-ray radiation due to the flare cited above, and that this in- 
fluence did not appear completely on the partial-reflection measurement, 

because this measurement was the result of an average over 12 minutes. The 

-3 -2 -1 

Solrad 9 satellite reported a flux of 1.0 x 10 erg. cm s , dropping 

rapidly, for the range of 1-8 A, at 12:35. At 13:10 this flux dropped to 

„ , -4 -2 -1 
4 x 10 erg. cm s 
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Figure 5.15 Electron-density profiles obtained from partial -reflection measurements (full 
line) and rocket measurements (traced line) at Wallops Island, Virginia, on 
January 31, 1972. 
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It should be noticed that no ground measurement made on January 31 
indicated a day of high absorption, including the vertical absorption at 
1.8 MHz. Finally, it is interesting to observe that the electron densities 
obtained by the rocket between 72 and 79 km are very high even for a day 
of high absorption. 

In Figure 5.16, the profile shown in full line was obtained by taking 
at each height the maximum electron density measured by rockets on previous 
days of high absorption [Mechtly et at . , 1973; Sechrist et al. 3 1969] . In 
traced line is the rocket profile obtained on January 31. 

Between 72 and 75 km, the profile on January 31 is higher than any 
one obtained previously in high absorption days. The gradient in electron 
density, however, appears between 82 and 84 km, and not below 80 km, as 
normally happens on high absorption days. Such features suggest the in- 
fluence of X-ray radiation on the ionization between 72 and 78 km. 

5.4 Measurements of Reflection Coefficients 

Since the gain of the antenna, approximately 16 dB, and the gain of the 
receiver were measured, it was possible to determine the values of the 
reflection coefficient as a function of height. The expression to determine 
the reflection coefficient is the following: 

R = 10 log 01*750) + L . * L * L. - G - P. (5.1) 

o o a p h a v 

2 

where i? Q is the reflection coefficient for the ordinary mode, (71^/50) is 
the power at the input of the receiver, I> a is the ionospheric absorption, 

1,^ is the path loss, L ^ is the ohmic loss in the cables between the trans- 
mitter and the antenna and between the antenna and the receiver, G is the 
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Figure 5.16 Electron-density profiles obtained by rocket measurement on January 31, 1972 (traced 
line), and by taking the maximum electron density at each height from previous 
measurements [Sechrist et at, 3 1969; Mechtly et at. , 1972]. 
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antenna gain and P is the transmitted power, all the above quantities 

~D 

given in dB» 

The calculated R is the average reflection coefficient over a height 
interval of 7.5 km, since 50 ysec pulses were used on the measurements. 
Figure 5,17 shows the reflection coefficient profiles for the ordinary 
mode, measured on December, 1971, in two days of high absorption (Dec, 12 
at 11:55 and Dec. 13 at 12:35), one day of low absorption (Dec. 15, 12:50) 
and one day of medium absorption (Dec, 19, 11:50). It can be observed that 
high absorption corresponds to higher reflection coefficients. The re- 
flection coefficient at 69 km, for example, is approximately four times 
greater on December 3 than December 15. Near 80 km the reflection coef- 
ficient shows a steep gradient, that corresponds to the gradient in 
electron density at the same heights. On the high absorption day this 
gradient appeared at a height 4 km lower than on low absorption day. 

Figure 5.18 shows the reflection coefficient profiles, for the 
ordinary mode, on several days of January 1972. It is observed a rather 
large variation of the reflection coefficient at a given height, for dif- 
ferent days. At 75 km, for example, the reflection coefficient on 
January 7 at 11:25 was 3.5 times greater than on January 25 at 11:35. 
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Figure 5.17 Reflection coefficient profiles for the ordinary mode, measured during 
the month of December, 1971. 
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Figure 5.18 Reflection coefficient profiles for the ordinary mode, measured during the 
month of January, 1972. 
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6. CONCLUSIONS AND SUGGESTIONS FOR FURTHER WORK 
6.1 Cone lusions 

The main objectives of the present work were to analyze the validity 
of the scattering theories and signal processing methods used in partial- 
reflection measurements, to suggest possible improvements, and to execute 
reliable measurements of the D region during the winter of 1971-1972. 

In Chapter 4 partial reflections were studied for a region where the 
scattering cross section changes rapidly with height. Such a situation 
happens in the D region near 80 km altitude, at the same heights where 
electron-density profiles measured by rockets show a steep gradient. 

Figures 4.4 to 4.6 show A q profiles measured during the month of August, 
1971, where large gradients are observed between 75 and 80 km. The 
electron-density profiles corresponding to the A profiles of Figures 4.4 
to 4.6 are shown in Figures 4.1 to 4.3. It is observed that at the heights 
where the A q profiles change rapidly with height there appear valleys in 
the electron-density profiles. Such valleys cannot be considered as real, 
but a consequence of the finite pulse width of the transmitted pulses and 
of the changes of the scattering cross section with height. In Section 4.1, 
it was shown that in a region where the electron-density profile is a 
constant, but the scattering cross section changes with height, partial- 
reflection measurements will produce fictitious valleys in the electron- 
density profile, as shown in Figures 4.10, 4.11 and 4.13. Comparison of 
the model calculations of Figures 4.10, 4.11 and 4.13 with the measured 
profiles of Figures 4.1 to 4.3 shows a great similarity in the observed 
valleys. The corresponding^ profiles show the same similarity. 

Two possible mechanisms of reflection that can be responsible for the 
gradients in scattering cross section near 80 km were analyzed in Chapter 4: 
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reflections produced by gradients in electron density and by random 
irregularities with a locally homogeneous background medium. Reflections 
produced by gradients in electron density were analyzed in Section 4.2. 
The calculated reflection coefficient for a gradient similar to that 
observed on rocket profiles was 5.4 x 10~ 6 . Measured reflection coeffi- 
cients, shown in Figures 5.25 to 5.29 are of the order of 10 -5 near 
80 km. This result indicates that at such heights the reflections are 
at least partially produced by gradients in electron density. The magni- 
tude of the reflection coefficient near 60 km indicates that at low alti- 
tudes the gradients are also an important factor in the production of 
reflections. Further evidences of reflections produced by gradients in 
electron density are also presented in Chapter 2, Table 2.1, that shows 
that the correlation coefficient between A signals, 1.5 or 3.0 km apart, 
increases with height above 75 km, and is very high above 82 km, and in 
Chapter 4, Figures 4.18 and 4.19, that shows that the number of reflec- 
tions observed at each height is a maximum near 80 km, and in Figures 4.4 
and 4.6 that show that at the same heights where gradients in electron 
density are observed there is a sharp increase in the A profiles. Re- 
flections produced by gradients in electron density can explain the exis- 
tence of preferred heights of reflections, frequently observed in previous 
measurements, and confirmed in the measurements made at the University of 
Illinois 

The dependence of the reflection coefficient on the size of the re- 
flector, shown in Figure 4.15, indicates that resonances can appear, and 
that they are not the same for the ordinary and extraordinary modes. To 
overcome the influence of such resonances, a large number of samples (at 
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least 500) must be used in the calculation of one electron-density profile, 
and all samples must be used, without rejecting signals with low signal- 
to-noise ratio or signals that reached the saturation of the receiver. In 
case of saturation, the data cannot be used. The result of rejecting part 
of the received data is shown in Figure 4,17. 

Reflections produced by random irregularities in a locally homogeneous 
background medium were discussed in Section 4.3. An expression of the 
scattering cross section for such a medium was developed and is presented 
in expression (4.74). Application of the developed expression to partial- 
reflection calculations shows that the measured electron density will 
depend on the power spectrum of the irregularities, at higher altitudes. 
Calculations of electron- density profiles, assuming that the real profile 
is that shown in Figure 4.26, for different power spectrums, are shown in 
Table 4,2, If the spectrum is a Gaussian spectrum with a correlation dis- 
tance greater than 50 m, serious errors can be committed in partial- 
reflection measurements, producing fictitious valleys. 

Possible contributions of random reflections to partial reflections 
near 80 km are discussed in Section 4,3.3. It is concluded that if the 
power spectrum is a Gaussian spectrum with a correlation distance less 
than 50 m, fluctuations of electron density of only 0.1 percent will be 
enough to produce reflections of the same order of magnitude of the 
observed ones. 

The effects discussed above will be minimized if narrow pulses are 
used in the measurements. As a consequence, pulses with a period of 
25 psec or less are recommended. Besides reducing the pulse width, two 
possible improvements in the signal processing methods are discussed in 



Section 4,5. The first one is to deconvolute the A q and A profiles. In 

the attempts made during the present work no reliable result was obtained 

from deconvolution. The second one is to use only maximums of the A and 

A^ profiles at each height in the calculations of electron density. 

Figure 4.28 shows two electron-density profiles obtained from the same 

data, by using only maximums of the A and A _ profiles, or by averaging 

0 cc 

the signals at heights separated by 1.5 km. The two methods can produce 
quite different results. 

Based on the theoretical studies described above, the partial- 
reflection system of the University of Illinois was set up to execute 
measurements of the D region during the winter 1971-1972. The results 
are presented in Chapter 5. The month of December, 1971 was characterized 
by a periodic variation of the absorption, with a period of approximately 
6 days, and by a large difference between the absorption of a low and a 
high absorption day. In the days of high absorption there was an increase 
of electron density all over the D region, down to 63 km. During the 
months of January and February, 1972 , the behavior of the D region was 
quite different. The daily variations of absorption were small, and no 
day of high absorption was observed. The variations of ionization were 
closely correlated to variations in geomagnetic activity. Comparison 
between rocket and partial-reflection measurements made on January 31, 

1973, shows good agreement between 78 and 84 km, but between 72 and 78 km 
the measurements did not agree. Possible influence of ionization by 
X-rays on the rocket measurement was suggested as a possible cause for 
this discrepancy. Measurements of the reflection coefficient, shown in 
Figures 5,17 and 5.18 indicate that its value changes on an average from 
10 ^ at 60 km to 10 ^ at 80 km. Above 80 km a strong gradient is observed. 
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Measurements of electron density during sunrise, made on October 1971, 

at Urbana, Illinois show that strong reflections, probably produced by 

gradients in electron density, are observed from the C layer. Reflections 

from the C layer permitted the measurement of collision frequency at lower 

7-1 

altitudesc At 63 km, a value of 2,0 x 10 s was obtained. 

6,2 Suggestions for Further Work 

For the improvement of the partial-reflection technique, a better 
understanding of the nature of the irregularities must be acquired. Sta- 
tistical studies of the reflections can produce useful information. Direct 
measurements by rockets carrying current probes, designed to measure both 
the background electron density and the fluctuations, is another source 
of information. 

The scattering theory developed in Chapter 4 should be extended to 

include an arbitrary geometry of the background medium, and not only a 

linear variation, as considered here. 

Measurements made simultaneously with different pulse widths should 

be made frequently, in order to verify the influence of the measured 

pulse width on the measured profiles. Pulse widths of less than 25 nsec 

should be included in the measurements. 

Techniques for deconvoluting the A and A profiles should be ini- 

proved. Comparative analysis of profiles obtained by taking only maximums 

of the A and A profiles, or by taking averages at heights separated by 

o x 

1.5 km should be made. 

Simultaneous rocket and partial-reflection measurements on days of 
high absorption, when normally valleys are observed on the partial- 
reflection measurements, would permit a better study of the causes of 
existence of such valleys. 
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APPENDIX I. SIMPLIFICATION OF THE SCATTERED ELECTRIC FIELD 


Let's take equation (4.53), that is repeated below: 
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Applying Gauss' theorem to the second term of the second member of 
equation (l.A) it results: 
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Taking n(h') - n(OJ << 1, equation (2. A) can be written as: 
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Applying Gauss' theorem to equation (5. A), it results: 
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Substituting this expression on equation (3. A) , and taking 
— -3/2 _ — 1/2 — -2 

n(?j') = n(h') • n( 0) " the following expression is obtained for E~ 
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exp(fc 0 rc(°)fl) kk n (0) exp (ik ( .n (O)fk) _ ^ ^ 
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Equation (7. A) shows that the second term of the second member is 
longitudinal (direction of m ) , with same magnitude and in opposite direction 
to the longitudinal component of the first term of the second member. 
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APPENDIX II. COMPUTER PROGRAMS AND SUB-ROUTINES 
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**********************PR0AX******************************** 
PROAX PRINTS AX/AO AND ELECTRON DENSITY PROFILES BETWEEN 
S0 AND 90 KM. VALUES OF AX AND AO ARE READ FROM DECTAPE. 
AVERAGE POWER OF ORDINARY AND EXTRAORDINARY SIGNALS AND 
OF NOISE ARE CALCULATED. THE AVERAGE POWER OF NOISE IS 
SUBTRACTED FROM THE AVERAGE POWER OF SIGNAL AT EVERY 
HEIGHT. FRAMES APE REJECTED IF THE AVERAGE NOISE FOP BOTH 
AX AND AO EXCEEDS A GIVEN NUMBER, BMXNS, IF THE DIFFERENCE 
BETWEEN A SAMPLE AT A GIVEN HEIGHT AND THE SAMPLE AT THE 
SAME HEIGHT IN THE FRAME I MED I A TELLY BEFORE EXCEEDS A VA- 
LUE RE1 OP PE2 , THE SAMPLE AT THAT HEIGHT, 1.5 KM ABOVE. 

AND 1.5 KM BELOW ARE REJECTED. 

***************** **************************** fc***^**^^^^ 

DIMENSION FNAM(2),A0(21),AX(21 ) ,A VA0(21 ) ,A VAX(21) , 

1X0(2!) ,IRJ0(21),IRJX(2I > ,DI FN0C4) ,DIFNX(4) , 

2R0(21 ) ,RX(2 1 ) ,B NO (4) ,P NXC4) ,BBN0(4) ,BRNX(4> 

CALL HEAD (0) 

C INITIAL VALUFS 

SN0=0. 

SNX-O. 

I R =P« 

IRNOrO 
I R NX -?> 

R MO = 0. 

B MX =0 . 

DO 15 I =| ,21 
XO(I ):0. 

A VA 0 ( I )rO. 

AVAX(I)rPI. 

AO(I) = 0. 

AXCI ):0. 

iRjoa >=a 
is iRJX(n = o 

DO 17 I =1 ,4 
BMOCI >=0. 

1 7 BNX(I) = fl. 

CALL VALUE 
CALL DI NIT 
WRITE( 5,20) 

2 7 F OR MA T (1 5 H WHICH DATAFILE) 

R FAD (4,70) FNA. M 
77 FORMAT (2A 5) 

CALL FSTA T (2 , F NA M , L OG ) 

I E (LOG . NE . O) GO TO 40 
VR I TF (5,75 ) FNA M 

75 FOR MAT (EH FILF ,2A5,J9H NOT FOUND ON DAT 2) 

GO TO 10 

47 VR I TF ( 5 , 4 1 ) 

•1! FORMAT (24H COLLISION FRFO . PROFILE/ 

1 25 H SUMMER, VI NTEP, OR EQUINOX) 

PEAD (4, 42)PFS? 

42 FORMAT (A 5) 

CALL SEEK (2 ,FNA M) 

UP I TE (5, 47 ) 

47 FOPMA T ( 1 4H MAXIMUM NOISE) 

READ (4,44)BMXNS 
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44 FOR MA T (FI 0.0) 

UR I TEC 6* 45) 

45 FOR MA T (2 1 H REJECTION BELOW 72KW) 

READ (4,44)RE1 

VR I TE ( 6 f 46) 

46 F0RMATC21K REJECTION ABOVE 72KM) 

PEAD(4,44)RE2 

KEOFO=0 
KEOFX =0 
ID -0 

50 DO 52 I =1 ,21 

R 0(1 )=A 0(1 ) 

52 BX(I)=AX(I) 

DO 54 I = 1,4 
BP NO (I ) =P NO (I ) 

54 BBNXd )=RNX(I ) 

0 DREAD READS 21 VALUES OF SIGNAL, BETWEEN 60 AND 90 KM, 

C A ND 4 VALUES OF NOISE TAKEN BETWEEN 45 AND 51 KM 

CALL DREAD(AO,PNO,IERR,ID,KEOFO) 

I R(KEOFO,EO .1 ) GO TO 114 

CALL DREAD (AX,BNX,IERR , ID, KEOFX) 

I E (KEOFX. EO • 1 ) GO TO 112 
C AVERAGE POWER OF NOISE OF THE FPAME 

BMEA MO r 0 . 

BMEA NX rPI, 

DO 56 I -1 ,4 

BMEANOrBMEA NO+BNO (I)**2 
6 6 P ME A NX rRMEA NX+BNX (I ) **2 

PMEANO=SOPT (BMEA. NO/4.) 

P MEA NX =SORT (BMEA NX/4 . ) 

C REJECTION OF FRAMES IF AVERAGE. .NOISE FOP BOTH AO AMD AX 

C EXCEEDS A VALUE RMXNS 

I F(PMEA NO ,GT .R MXNS . A ND .BMEA NX .GT .BMXNS) GO TO 100 
C CALCULATION OF TOTAL AVERAGE NOISE 

DO 5* 1=1,4 

D I RMO ( I ) rP NO (I ) -B B NO (I ) 

62 DIF NX (I) :R NX (I) -PB NX (I) 

DO 66 Ir| ,2 
J=l + 1 
K=I+2 

C REJECTION OF ORDINARY SAMPLE IF FADING EXCEEDS 

C R El, OP IE SATURATION IS PEACHED 

I F( DI FNOCI ) .GT.PEI .OP .D I RNO ( J) .GT .R El . OR .DI FNO(K) .GT.REI .(F . 
lBNO(J) .GT.51 0..OR.BNX(J) ,GT.5Jfl.)G0 TO 60 
BM0=RM0+RN0(J)**2 

00 TO 62 

C MUMPER OP SAMPLES OF ORDINARY NOISE REJECTED 

60 IPNOrIRNO+l 

I E (? M0( J) .GT .5 1 0 . .OR ,P NX ( J) .GT .5 1 0. ) GO TO 64 
62 I F (DI F NX ( I ) .GT.PEI .OP.DIFNX(j) .GT.PEI .OF.DI FNX(K) .GT.PEI) 

1 GO TO 64 
RMX=PMX+RNX(J)**2 
GO TO 66 

64 I R NX =IR HX+ 1 

66 CONTINUE 

62 DO P4 I =1 ,21 

0 REJECTION OF SAMPLE OF ORDINARY SIGNAL IF FADING EXCEEDS 
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C 


77 ! 

C 

7 ? 

C 

74 

76 


7 ? 


87 ! 


82 

84 


c 

1 88 

1 ! 8 


C 

I 18 
1 ! 4 


0 


116 


C 

r. 


A VALUE REI OR RE2 , OR IF SATURATION IS REACHED 
J=I + 1 
L = I -I 

IE(L.EO.0)L = 1 
I F ( J.GT .21 ) J=2 I 
DIFE=A0(1)-B0(I) 

DIFEA =A 0<L)-BOCL) 

D I FEB =AO(J)-BO(J) 

I F(I .GT .9)00 TO 78 

I FCDI FE.GT .REI .OR .DIFEA .GT .RE 1 .OP .DI EE3 .GT .REI .OR. 
JAO(I). GT «510«. OP. A XCD.GT. 510.) GO TO 74 

0 0 TO 72 

IFCDin- .GT.RE2.0R.DIEFA.GT.RE2.0R.DIEFB.GT.RE2.0P. 
1AOCI).GT.510..OR.AX(I).GT.510.)GO TO 74 
AVERAGE POWER OF ORDINARY SIGNAL 
A VA 0 ( I ) r A VA 0 ( I ) +A 0 ( I ) **2 
GO TO 7fi 

HUMBER OF SAMPLES OF ORDINARY SIGNAL RFJECTED 
IRJO(I)-IRJO(I)+! 

IF(AO(I) .GT.510..OP .AX(I> .GT.510.)GO TO 82 
D I FF=A X( I >-BX (I > 

DIFEA =AX(L>-eX(L> 

DI FEB=AX(J)-BX(J> 

1 F(I .GT.9.A ND .1 .LT.l 6)00 TO 78 

IF(DIFE.0T.PE!.0R.DIFEA.0T.RE1.0R.DIFEB.GT.RE1)G0 TO 82 
00 TO 80 

I F(DI FE.GT ,RF2. .OP ,DI FEA .GT *RE2 .OR .DI FEB.0T.RE2 ) GO TO 82 
A VA X Cl ) =A VA XCI )+AX (I >**2 
GO TO 84 

IRJX(n = IRJX(I) + l 
CONTI HUE 
GO TO 1 10 

NUMBER OF FRAMES REJECTED DUE TO EXCESSIVE NCISE 
IP: IR+1 

SHO=FMEANO**2+SNO 
S MX=BMEANX**2+SNX 
GO TO 50 

NUMBER OF SAMPLFS TAKEN 
I D =ID-I 
lD=ID/2 
BID =ID 

TOTAL NUMBER OF PFJFCTFD SAMPLFS AT EACH HEIGHT 

DO 116 I rt ,21 

IRJOCl) = IRJO(I)-fIP 

IPJXCI)=IFJXCI>+IR 

IPMO = IR H0+2+IR 

IR HX=!RNX+2*IR 

AVNO:SORT(SNO/PID> 

A VWX=SORT CSNX/P ID ) 

R NP = ID*2.-IPN0 
P MX -ID *2 -I R NX 

AVERAGE POWER OF SIGNAL SUBTRACTED FROM THE A VERA OE 
POWER OF NOISE 
DO N8 1=1,21 
RSAMOrlD-IPJOCl) 

R SA MX =ID-I R JXCI ) 

A VOC =A VA 0 ( I ) /P SA M 0 -3 MO /R NO 





200 


A VXC :A VA X ( I > /R $A MX *6 MX /R NX 
A VA 0 (I ) : (ABS (A VOC > /A VOC ) *S0R T (ABS (A VOC > ) 

A VAX (I ) = (ABS (A VXC) /A VXC)*SQRT (ABS (A VXC) ) 

UP CONTINUE 

PMOrSQRT (BMO/R NO) 

BMX:SQRT(BMX/RNX> 

CALL HEAD ( l ) 

WRITE (6, 120)AVNO,AVNX,BMO,IRNO,BMX,IRNX 
lf»fl FORMAT (2X, 1 7HA VERA GE NOISE 0RD.F8.I /2X,! 7HA VERAGE NOISE FXT,IX 
! ,F 7.1 //2X,9K NOISE OPD.FP.l ,15 » J IH REJECTI 0NS/2X, 

29HN0ISE EXT ,FP.1,I5,11H REJECTIONS) 

VRITE(6,!22)ID,IR 

122 FORMA T(//1X,I4,14H SAMPLES TAKEN, 5X, 15 , 1 SH FRAMES REJECTED// 

19H REJECTED, 2X,9H REJECTED ,2X,SHHEIGHT ,2X*CHAV. A0,2X, 

2CHAV. AX/4X,3H0RD,RX,3HEXT> 

HT =5 P.5 
DO ISC 1:1,21 
HT =HT+ 1 .5 

VRTTFCS,124)IRJO<I) ,IRJX(I> ,HT,AVAO(I) ,AVAX(I) 
t 2.4 F0PMAT(3X,I4,7X»I4,3X,F5.1 ,3X,F6.1 ,2X»F6 e l) 

12* CONTINUE 

CALL HEAD CD 
C AX/A 0 RATIOS 

DO I2P 1:1, 21 

IF(AVAO(I) ,LF .O.O.OR.AVAX(I) .LE.O.fl) 00 TO 12P 
XOCI>=AVAX(I)/AVAO(I) 

12P CONTINUE 

MI NI T=I D/1 50 
C AX/AO PROFILE 

CALL PLOTI (X0,1 ,21 ,MI NIT) 

C CALC CALCULATES AND PRINTS ELECTRON DENSITY PROFILE 

CALL CALC(X0,1 ,20,RESP) 

GO TO 10 

STOP 

END 
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**********************pp 0A XC*********************** ******** 
PROA XC PRINTS AX/AO AND ELECTRON DENSITY PROFILES AND 
NUMBER OF MAXI MUMS OBSERVED IN THE REFLECTED SIGNALS 
AT EACH HEIGHT. THE SIGNAL PROCESSING METHOD IS THE 
SAME AS USED IN THE PROGRAM PROAX. 
*********************************************************** 


I 0 

c 


16 
1 7 

2 0 
30 


53 

40 

41 


42 


45 

44 

45 


DIMENSION FNAM(2),A0(21> ,AX(2l) ,A VA 0 (21 ) ,A VAX <2! > ,1 REJ (2 1 ) , 
1X0(21 ) «HE 1 0(21 ) ,HEIX(21 ) »IRJ0C2l ) »IRJX(21 ) , 

230(21 ) ,BX(2l > ,BN0<4) ,9 NX (4) ,3BN0(4) ,BRNX(4) ,DIFN0(4) ,DIFNX(4) 
CALL HEAD (0) 

INITIAL VALUES 
SNO=0. 

S NX =0. 

IR =0 
IRNOzO 
IRNX=0 
BMO=0, 

BMX =0, 

DO 16 I =1 ,21 
HEIO(I>=0. 

HEIX(I):0. 

XO(I)r0. 

AVAO(I> = 0. 

AVAX(I) = PI. 

AO(I) = fl, 

AX(I ) = 0. 

IRJO(I) = 0 
IRJXCI ):PI 
DO 17 1=1 ,4 
9NO(I) = 0. 

8NX(I) = 0. 

CALL VALUE 
CALL DINIT 
WRI TE ( 6,20) 

FORMAT (1 5H WHICH DATAFILE) 

READ (4,30) FNA M 
FORMAT (2A5) 


CALL FSTAT(2, FNAM.LOG) 

I F(LOG . NE . O) GO TO 40 
WRITE(6,35)FNAM 

FOR MAT (6H FILE ,2A5,J9H NOT FOUND ON DAT 
GO TO 10 
WR I TE ( 6, 4 1 ) 

FORMAT (24H COLLISION FRED. PROFILE/ 

1 25H SUMMER, WINTER, OP EQUINOX) 

READ (4, 42 >RESP 
FOR MAT (A 5) 

CALL SEEK (2 ,F NA M) 

WRI TE ( 6 , 43) 

FORMATd 4H MAXIMUM NOISE) 

READ (4,44)5MXNS 
FORMAT (FJ 0.0) 

WR I TE ( 6, 45 ) 

FORMAT (2IH REJECTION BELOW 72KM) 

READ(4,44)PE1 

WR I TE ( 6, 46) 


2 ) 


original was® 

OF POOR SPAIH* 
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46 F0RMAK21H REJECTION ABOVE 72KM) 

READ (4 ,44>RE2 
KEOFO=0 
KEOFX=0 
ID=0 

50 DO 52 1=1.21 

BO<I) = AO«) 

52 BX(I)=AX(I) 

DO 54 I =! ,4 
PBNO (I ) =B NO (I ) 

54 BBNX<I)=BNX<1) 

C DREAD READS 21 VALUES OF SIGNAL. BETWEEN 60 AND 90 KM, 

C AND 4 VALUES OF NOISE TAKEN BETWEEN 45 AND 51 KM* 

CALL DREAD (AO,B N0,I ERR ,1 D ,KEOFO) 

I F (KEOFO .EQ • t ) GO TO 114 

CALL DREAD (AX ,B NX , I ERR .ID .KEOFX) 

I F (KEOFXeEG *1 > GO TO 112 
C AVERAGE POWER OF NOISE OF THE FRAME 

RMEANO=0. 

BMEA NX =0. 

DO 56 I =1 ,4 

RMEA NOrBMEA NO+9NO (I ) **2 
*5 6 B MEA NX =8MEA NX+B NX (I ) **2 

BMEA NO=SQRT (BMEA NO/4.) 

B MEA NX =S0RT (BMEA NX/4.) 

C REJECTION OF FRAMES IF AVERAGE NOISE FOR BOTH AO AND AX 

C EXCEEDS A VALUE BMXNS 

I F(BMEA NO.GT.PMXNS.A ND .BMEA NX.GT.BMXNS) GO TO 130 
C CALCULATION OF TOTAL AVERAGE NOISE 

DO 5F 1 = 1,4 

DIFN0(I)=8N0(I)-B9N0(I) 

58 DIFNX(I)=BNX(I)“BBNX(I) 

DO 66 I =1 ,2 
J=I+1 
K =1+8 

C REJECTION OF ORDINARY SAMPLE IF FADING EXCEEDS 

0 R El, OR IF SATURATION IS REACHED 

IFCDIFNO(I) .GT.RE1 .OR.DIFNO(J) .GT.RE1 .OR.DIFNOCK) .GT.RE1 .OR. 
IBNO(J) .GT.510..OR.BNXCJ) .GT.51fl).)G0 TO 60 
8M0=BM0+9N0(.J>**2 
GO TO 62 

C NUMBER OF SAMPLES OF ORDI NAPY NOISE REJECTED 

60 IPN0=IRN0+1 

I F(BNO ( J) .GT.5 1 0. .OR .BNX ( J) *GT .5 1 0* ) GO TO 64 
62 IF(DIFNX(I).GT.RE1 . OR .DI FNX (J) *GT *R El .OR.DIFNX(K) .GT.REl) 

1 GO TO 64 
PMX=RMX+B NX (J)**2 
GO TO 66 

64 I R NX =1 R NX+ 1 

66 CONTINUE 

6R DO P4 1=1,21 

H REJECTION OF SAMPLE OF ORDINARY SIGNAL IF FADING EXCEEDS 

C A VALUE RE1 OR RE2 , OR IF SATURATION IS REACHED 

J = I+t 
L=I-l 

I F (L.EQ .0)L=l 
I F ( J.GT .2 1 ) J=2 1 
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70 

C 

f> nr > I 


0 

272 

c 

72 

r 

74 

76 


78 
2 71 4 


27 f 

87 ) 

7 ? 

«4 

C 

1 77 
1 10 


112 
1 1 4 

V. 

1 1 f 


\ 


DIFE=A0(I)-B0d) 

DIFEA=A0(L)-B0CL> 

DIFEB=A O(J)-BOCJ) 

I F Cl ,GT.9)G0 TO 77 

IF (DIFE.GT.REl .OR.DIFEA.GT ♦R El .OP.DI FEB.GT.RE1 .OR. 

1A0CI) .GT.510..OR.AXCI) .GT.510.)GO TO 74 
GO TO 200 

I F(DI FE.GT.RE2 .OR ,DI FEA .GT.RE2 .OR .01 FEB .GT .RE2.0R . 
!AO(I).GT.510..OR.AX(I) .GT.5I 0.) GO TO 74 
DETERMINATION OF MAXI MUMS IN THE REFLECTED SIGNAL 
EJ:AO(I)-0.5*AO(L)-0.5*AO(J> 

IF (A 0(1 ) .GT.A 0 ( J) ,AND.AO(I).GT.AO(L).AND. 
lEJ.GT.,05*AO(I),AVD.AOd) ,GT .2 ,*BMFA NO) GO TO 202 
GO TO 72 

NUMBER OF MAXI MUMS 
HEI 0 (I ) = HEI 0(1 ) + l . 

AVERAGE POWER OF ARDINARY SIGNAL 
AVAO(I )rAVAO(I)+AO(I )**2 
GO TO 76 

NUMBER OF SAMPLES OF ORDINARY SIGNAL REJECTED 
IRJ0(I>=IPJ0(I>+1 

I F (A 0 ( I ) ,GT.5!0..OR.AX(I) *GT .510.) GO TO 82 
DIFE=AX(I>-BX(I> 

DIFEA=AX(L)-BX(L) 

DIFEB=AX(J)-BX(J> 

IF(I .GT.9.AND.I *LT , I 6) GO TO 78 

IFCDIFE.GT.PEl .0R.DIFEA.GT.RE1 .0R.DIFEB.GT.RE1 )G0 TO 82 
GO TO 204 

IF(DIFE.GT.RE2.0R.DIFFA .GT ,RE2 . OR .DI FEB.GT ,RE2 ) GO TO 82 
EJrAXd)-0.5*AXCL)-0.5*AX(J> 

IF(AX(I).GT.AX(L) .A ND ,AX(I ) .GT.AX(J) .AND .EJ.GT. ,05*AX(I ) * 
lAND.AX(I) .0T.2.*BMEANX)G0 TO 206 
GO TO 80 

HEIX(I):HEIX(I) + 1 . 

A VAX (I )=A VAXd)+AX(I)**2 
GC TO 84 

IRJX<I)=IRJXCI)+1 
CONTI NUE 
GO TO 11 0 

NUMBER OF FRAMES REJECTED DUE TO EXCESSIVE NOISE 
IR=IR + 1 

SN0=BMEAN0**2+SN0 
SNX:£MEANX**2+SNX 
GO TO 50 
ID -ID-1 


ID=ID/2 

PID=ID 

TOTAL NUMBER OF REJECTED SAMPLES AT EACH HEIGHT 


DO 116 I =1 ,21 
IRJO(I)=IRJO(I )+IR 
IP JX(I ) = IRJX d >+IF 
IRN0 = IRNCH-2*IP 
IRNXrIPNX+2*IP 
AVNO=SORT(SNO/BID) 
A VNX :SORT (SNX/B I D ) 
MR = (ID /I 00) *30 
CALL HEAD ( 1 ) 


ORIGINAL PAGE IS 
0$ POOR gUAXOT; 
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WRITE(S,208) 

ZM FORMAT <20X,25HHEIGHTS OF REFLECTION, AO) 

CALL PLOTF(60.,88.5,20,HEIO,NR) 

WRI TE<G,21 0) 

FOR MA T ( //20X,2.5HHEI GHTS OF REFLECTION, AX) 

CALL PLOTF(G0.,88.5,20,HEIX,NR) 

R N0=ID*2-IRN0 
R NX =1 0*2-1 R NX 

C AVERAGE POWER OF SIGNAL SUBTRACTED FROM AVERAGE 

C POWER OF NOISE 

DO 118 I =1 ,21 
RSAMOrlD-IRJO(I) 

RSAMX=ID-IRJX(I > 

AVCC=AVAO(I >/PSAMO-BMO/RNO 
AVXC =AVAXCI)/PSAMX-BMX/RNX 
AVAO(I)=CABS(AVOC)/AVOC)*SORT(ABS(AVOO) 
AVAX(I)=(ABS(AVXC)/AVXC)*SORT(ABS(AVXC)) 

IIP CONTINUE 

RMO=S0RT(BMO/FMO) 

PMX=S0RT (BMX/R NX) 

CALL HEAD < I ) 

VRITE(fi,12P!)AVN0,AVNX,BM0,IRN0,BMX,IRNX 
12 0 FORMAT (2 X, 1 7HA VERA GE NOISE 0RD,F8.1 /2X,1 7HAVERAGE NOISE EXT, IX 

1 ,F7,I //2X,9HN0ISE ORD ,F8. 1 ,1 5 , 1 IH REJECTI 0NS/2X, 

29HN0ISE EXTjFB.l ,I5,1IH REJECTIONS) 

WRITE(6,I22)ID,IR 

122 F0RMAT(//1X,I4,14H SAMPLES TA KEN,5X,1 5 , 1 6H FRAMES REJECTED// 

19H RE JECTED ,2X, 9H REJECTED , 2X, 6HHEI GHT,2X, 6HA V. A0,2X, 

2SHAV. AX/4X ,3H0FD, PX,3HEXT) 

HT=5P.5 
DO 12S I =1 ,21 
H T -HT+ 1 *5 

WRITE(S*124)IRJ0(I),IRJX(I),HT,AVA0<I) ,AVAXd) 

124 FORMAT <3X,I4 e 7X,I4,3X,F5.1 ,3X,F6. 1 ,2X,F6.1) 

12F CONTINUE 

CALL HEAD (1 ) 

C AX/AO RATIOS 

DO 12F I =1 ,21 

IF(AVAO(I).LE.0.0.OR.AVAX(I).LE.0,0)GO TO 12P 
XOCI >=AVAX(I)/AVAO(I) 

12 F CONTINUE 

MI NIT=ID/150 
C AX/AO PROFILE 

CALL PLOT! (X0,1 ,2! ,MI NIT) 

C CALC CALCULATES AND PRINTS ELECTRON DENSI TV PROFILE 

CALL CA LC <X0, 1 ,20 ,R ESP ) 

GO TO 10 

STOP 

END 
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**********************pp OAT******************************** 
PROAT PRINTS AX/AO AND ELECTRON DENSITY PROFILES FOR PAR- 
TIAL REFLECTION DATA TAKEN WITH A PROGRAMMED ATTENUATOR 
THAT OPERATES ON ALTERNATE FRAMES. THE AO SIGNALS USED 
ABOVE A HEIGHT AHT ARE THE ATTENUATED SIGNALS. THE SIG- 
NAL PROCESSING METHOD IS THE SAME AS USED IN PROGRAM 
PROAX 

*********************************************************** 


10 

C 


IS 

17 

2* 

30 

3? 


40 

41 


42 


43 


DIMENSION ENA M(2), A 0(21), AX (21) ,A VA 0(2 1 ) ,A VAX(21 ) , 
1X0(21),IRJ0(21),IPJX(21) ,DIFN0(4),DIFNX(4), 

2P0C2I) ,BX(21),BN0(4),B NX (4) ,BBN0(4) ,B8NX(4) 
3,IRJ0T(2I),B0T (21>,BNXN(4) ,A XNC21 ) ,A OT (21 ) ,A XT (21 ) 
4.BN0T (4) ,B NXT (4 ) ,BB NOT (4) ,DIFN0T(4) ,AVA0T(2l ) 

CALL HEAD (0) 

INITIAL VALUES 
SNO = 0. 

S NX -0. 

IR=0 
IRMO=0 
I R NX=0 
iRNOTrO 


BMO:0. 

BMOT = 0, 

BMX=0. 

DO IS I =1 ,21 
XO(I ) = 0. 

AVAO(I)=0. 

AVAX(I> = 0. 

A 0(1 >-0. 

AX(I) = 0. 

IRJO(I)=0 
IRJOT (I ) = 0 
I R JX (I ) = 0 
DO 17 1=1,4 
B NO(I ):0« 

B NX ( I ) = 0 . 

CALL VALUE 
CALL DINIT 
WRITE(S,20) 

F OR MA T ( 1 5H WHICH DATAFILE) 

READ (4 , 30) FNA M 
FORMAT (2A5) 

CALL FSTAT (2, F NAM, LOG) 

IF(LOG.NE.0)GO TO 40 
VR I TE ( 6 , 35 ) FNA M 

F OR MA T ( SH FILE ,2A5,19H NOT FOUND ON DAT 2) 
GO TO 10 


WP I TE ( 6, 4 1 ) 

FORMAT (24H COLLISION FPFO. PROFILE/ 
1 25H SUMMER, VI NTFR, OP FOUI NOX) 

PEAD (4,42)RESP 
FORMAT (A 5) 

CALL SEEK (2 ,FNA M) 

VP I TF ( S , 43 ) 

FORMAT! 1 4H MAXIMUM NOISE) 

READ (4,44)BMXNS 





206 


WRITE (6,300) 

300 F OR MA T (22 H HEIGHT STARTING ATTEN) 

R EAD (4 ,44>AHT 
IT=(AHT-58.5) /\ .5-1 . 

44 FORMAT (FI 0.0) 

WRITE (6, 45) 

45 FOR MAT (2 1H REJECTION BELOW 72KM) 

READ (4„44)RE1 

WR I TE ( 6,46) 

Af F0RMAT(2IH REJECTION ABOVE 72KM) 

READ (4,44>RE2 
KEOFO=0 
KEOFX =0 
ID =0 

CALL DREAD (A 0,B NO,I ERR, I D ,KEOFO) 

CALL DREAD (AX,B NX,I ERR , I D ,KEOFX) 

C A LL DREAD (A OT ,B NOT , I FRR , I D,KEOFO) 

CALL DREAD (A XT ,B NXT,I ERR ,1 D ,KEOFX) 

IDM=4 

IF(AOT (9). LT. A 0(9)) GO TO 50 
CALL DREAD (A OT,B NOT ,IERR ,ID,KEOFO) 

CALL DREAD (AXT,BNXT,IERR, ID, KEOFX) 

IDtf = 6 

50 DO 52 I =1 ,21 

BO(I)=AO(I) 

BOT(I)=AOT(I) 

52 9X(I)=AX<I> 

DO 54 I =1 ,4 
PBNO(I)=BNO(I) 

BB NOT (I )rRNOT (I ) 

54 BBNX(I)=BNX(I > 

C DREAD READS 21 VALUES OF SIGNAL, BETWEEN 60 AND 90 KM, 

r AND 4 VALUES OF NOISE TAKEN BETWEEN 45 AND 51 KM 

CALL DREAD (AO,BNO,I ERR,ID,KE0F0) 

I F (KEOFO.EO .1 ) GO TO 326 

CALL DREAD (AX,B NX, I ERR, ID, KEOFX) 

I F (KEOFX *EQ .1 ) GO TO 320 

CALL DREAD (AOT,B NOT, I FRR , I D ,KEOFO) 

I F ( KFOF 0 . EO . 1 ) GO TO 322 

CALL DREAD (AXT,B NXT,I ERR ,1 D, KEOFX) 

I F (KEOFX .EO ,1 ) GO TO 324 
C AVERAGE POWER OF NOISE OF THE FRAME 

B MEA NO = 0. 

B MEA NX=0. 

DO 56 I =1 ,4 

9 MEA NOrBMEA NO+BNO(I >**2 
5 6 9 ME A NX =BMEA NX+PNX (I ) **2 

B MEA NO=SQPT (BMF A NO /4 . ) 

P MEA MX =SORT (PMFA NX/4 . ) 

C REJECTION OF FRAMES IF AVERAGE NOISE FOR BOTH AO AND AX 

C EXCEEDS A VALUE BMXNS 

I F (BMEA NO .GT.BMXNS .A ND ,RMFA NX.GT .BMXNS) GO TO 100 
r. CALCULATION OF TOTAL AVERAGE NOISE 

DO 5* 1=1 ,4 

DIENO(I)=BNO(I)-BBNO(I) 

D I F NOT (I ) =B NOT (I) -BB NOT (I ) 

5? DIFNX(I) = BNX(I)-BBNX(I) 
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C 

C 


C 

60 

s?. 


64 

66 


302 

304 

6P 

C 

c 


73 

c 

7 ? 

c 

74 
76 


306 

30P 


DO 304 1=1 ,2 

J=I+1 

K=I+2 

REJECTION OF ORDINARY SAMPLE IF FADING EXCEEDS 
RE1, OR IF SATURATION IS REACHED 

IF(DIFNO(I).GT.REl .OR.DIFNO(J) .GT.RE1 .OR.DIFNOUO .GT.RE1 .OR « 
IB NO (J) ,GT ,? 1 0 , . OR .B NX (J) »GT .5 1 0.) GO TO 60 
BMO=BMO+BNO(J)**2 
GO TO 62 

NUMBER OP SAMPLES OF ORDINARY NOISE REJECTED 
IR NO=IR NO+1 

IF(BNO(J> ,GT.510..OR.PNX(J) ,GT.510.)GO TO 64 
IF(DIFNX(I> .GT.PE1 .OR.DIFNX(J) .GT .RE! .OR.DIFNX(K) .GT.RE1 ) 

I GO TO 64 
RMX=RMX+B NX ( J)**2 


GO TO 66 
I R NX =1 P NX+ 1 

I F (BNOT (J).GT.510..OR.B NXT ( J> .GT.510.5GO TO 302 

IF(DIFNOT(I) .GT.RE1 .OR.DIFNOT(J) .GT .RE 1 . OR .DI FNOT (K> . 

2GT.PEDG0 TO 302 

BM0T=BM0T+BN0T(J)**2 

GO TO 304 

IRNOT=IRNOT+l 

CONTINUE 

DO 76 1=1 ,IT 

REJECTION OF SAMPLE OF ORDINARY SIGNAL IF FADING EXCEEDS 
A VALUE RE I OR RE2, OR IF SATURATION IS REACHED 
J = I + 1 
1 = 1-1 

I F(L.EG.0)L = 1 
I FC J.GT .2 1 ) J=2 1 
D I FE=A 0(1 ) -BO (I ) 

DIFEA =A OCL)-BO(L) 

D I FEB =A 0 ( J>-BO ( J) 

I FCl .GT.9) GO TO 70 

I F(DI FE.GT.RF1 .OR. DIFEA .GT.RE1 .OR ,DI FEB.fiT.REl .OR . 
1AO(I).GT.510..OR.AX(I).GT.510.)GO TO 74 
GO TO 72 

I FCDI FE.GT.PE2 .OR .DIFEA .GT.RE2 .OP .DI FEB .GT .RE2 .OR . 

1 A 0 <1 >.GT.510..OP.AX(I) .GT .51 0, ) GO TO 74 
AVERAGE ROVER OF ORDINARY SIGNAL 
A VA 0(1 > = A VA 0( I )+A 0 (I ) **?. 

GO TO 76 

NUMBER OF SAMPLES OF ORDINARY SIGNAL REJECTED 

IRJO(I )=IPJO(n + l 

CONTI NUE 

DO B4 I =1 ,21 

I F Cl .GT .1 T) GO TO 306 

IF(AO(I ).GT.5] 0..OP .AX(I) ,GT.6!0.)GO TO P2 
GO TO 30P 

I F (AX( I ) .GT .5 1 0. ) GO TO R2 
J=I + I 
L =1 -1 


IF(L.ED.0)L=l 
I F (J.GT ,2 1 ) J=2 1 
D I FE=A X (I ) -BX(I ) 
DIFEA =AX(L)-BX(L) 





D I FEB =A X ( J> -BX f J) 

IFCI .GT.9.A ND.I ,LT.I6)G0 TO 78 

I E(DIFE.GT.RE1 .OR.DIFEA .GT.REl .OR ,D1 FEB.GT ,RE1) GO TO 82 
GO TO 80 

78 IF(DIFE.GT.RE2.0R.DIFEA.GT.RE2.0F.DIFEB,GT.RE2)G0 TO 82 

80 A VAX(I)=AVAXU)+AXa)**2 

GO TO 84 

8?. IRJX<I)=IRJX(I>+1 

84 CONTINUE 

ITM=22-IT 
DO 314 1=1 .1 TM 
IJ = I+IT-I 
J=IJ+1 
L =1 J-l 

IF(L.EQ.0)L=! 

IF(J.GT.21)J=21 
D I FF. =A OT (I J) “BOT (I J) 

Dl FFA =A OT (L) -90T (L ) 

DIFEB=A0T(J)-80T(J) 

I F (DI FFA .GT .RE t .OP .01 FEB .GT.RE ! .OR ,DI FE.GT .RE l .OR . 
2A0TCIJ) .GT.510.)GO TO 312 
AVA OT (I J) =A VA OTCIJ)+A 0T(IJ)**2 
GO TO 314 

312 IRJOT(I J)=IRJ0T(IJ) + 1 

314 CONTINUE 

GO TO 110 

C NUMBER OF FRAMES REJECTED DUE TO EXCESSIVE NOISE 

100 IR =IR+ 1 

110 SN0=8MEAN0**2+SN0 

S NX =BMEANX**2+SNX 

00 TO 326 

320 ID =ID-1 

GO TO 326 
322 I D =ID -2 

GO TO 326 
324 I D =ID -3 

326 ID=(ID-lDM>/4 

BID =ID 

C TOTAL NUMBER OF REJECTED SAMPLES AT EACH HEIGHT 

DO 116 Irl ,21 
IRJ0(I) = IRJOCn+IR 
I RJOT <1 ) = IR JOT (I)+IR 
116 IRJX(I>=IRJX(I)+IF 

I R NO=IR N0+2*IP 
IR NOT =1 R N0T+2*IR 
IRMX=IRNX+2*IR 
AVNO=SQRT(SNO/BID) 

A VNX =SQRT (SNX/B ID ) 

R N0=ID*2-IFN0 
PNX=IP*2-IPNX 
R NOT =ID*2-I P NOT 

0 AVERAGE POWER OF SIGNAL SUBTRACTED FROM THE AVERAGE 

C POWER OF NOISE 

DO 118 1=1 .21 

1 F(I .GT .1 T> GO TO 328 
RSAMO=ID-IRJO(I) 

AV0C=AVA0(I) /RSAMO-BMO/RNO 
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A VA 0 (I ) r ( ABS (A VOC ) /A VOC ) *S0RT (A BS ( A VOC ) ) 

328 RSAMX=ID-IRJX(I) 

A VXC =A VA X ( I > /R SA MX -B MX/R NX 
A VA XCI >= (ABS (AVXC) /AVXC)*SQRT (ADS (A VXC) ) 

118 CONTINUE 

DO 33fl I =1 ,ITM 
IJ=I + IT-1 

PSAMOHD-IRJOTCI J> 

AVOC =A VA OT (I J) ^ SAMO-B MOT /R NOT 

AVA 0T(IJ)=(ARS (AVOC ) /AVOC)* SORT CABS (AVOC) ) 

330 IPJO(IJ)=IRJOT (I J) 

C ORR =A VA 0 ( I T ) /A VA OT (IT) 

D COR r20.*AL0G I 0 (C ORR ) 

ITC=21-IT 
DO 332 I =1 ,ITC 
I JC=I T+I 

332 A VAO(IJC)=A VAOT (I JC)*CORR 

B MO -SORT (BMO/P NO) 

RMX=SOPT(BMX/PNX> 

CALL HEAD (1 ) 

WRITE (G, 120)AVNO,AVNX,BMO,IRNO,BMX,IRNX 
120 F0RMAT(2X, 17HAVERAGE NOISE 0RD,F8. 1 /2X, 1 7HA VERAOE NOISE EyT,lX 
1 ,F7.l//2X,RHN0ISE 0RD,F8, I ,15,1 1H REJECTI 0NS/2X, 

29HN0ISE EXT,R8.1 ,I5,HH REJECTIONS) 

WRITER, 122>ID,IR,DC0R 

122 FORMAT (//l X,I 4 , 1 4H SAMPLES TA KEN,5X ,15,1 *SH FRAMES REJECTED/ 

2 F6.2 , 1 5H DB ATTENUATION// 

3 9H REJECTED, 2X,9K REJECTED ,2 X, 6HHEI GHT ,2X, SHA V. A0,2X, 

2GHA V, AX/4X ,3 WORD, 8X,3HEXT) 

HT=5 8 ,5 
DO 128 1=1 ,21 
HT=HT+1 .5 

WRITE(F,12A)IRjn(I) ,IPJX<I) ,HT,AVAO(I),AVAX(I ) 

124 F0RMAT{3X,I4,TX,I4,3X,F‘'.!,3X,F6.I ,2X,F 8,1) 

178 CONTINUE 

CALL HEAD < 1 ) 

C AX/AO RATIOS 

DO 128 I =1 ,21 

I F (A VA 0 ( I ) ,LE ,0, P! ,0R , A VA X (I ) .LE ,0.2)00 TO 128 
XO(I)=AVAX(I)/AVAO(I) 

12? CONTINUE 

*INIT=ID/I5fl 
C AX/AO PROFILE 

CALL PLOTI (X0,1 ,21 ,MINIT) 

0 CALC CALCULATES A ND PPINTS ELECTRON DENSITY PROFILE 

CALL CALC(X0,1 ,20,PFSP> 

00 TO IP) 

STOP 

END 
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c ***************+*****SUBR0UTINE OREAD********************* 

C DREAD READS 21 SAMPLES OF SIGNAL AND A SAMPLES OF NOISE 

C FROM DECTAPE. THE OUTPUT VOLTAGES ARE TRANSFORMED INTO 

C INPUT VOLTAGES 9Y LI NAP. 

c ********************************************************** 

0 

SUP-ROUT! NE DRFAD(A ,RMEA N,IERR,ID,KEOF> 

DIMENSION A (2n,IDAT(26) ,BMEANC4> 

KEOF=P 

CALL DUMPP (IDAT, NEGF) 

C CHECK ID CONSECUTIVE 

I F (ID -I DAT ( 1 > + l ) 10,15,10 
C CHECK FOR EOF 

Ifl IF(IDAT(I).NE. 77777) GO TO I ERR 

C F.O.F. 

K F OF ; 1 
GO TO 2JO 
15 ID-IDAT Cl > 

AO DO A?. MINTS, 2 6 

mFVFtMI N-5 
A (MFVF)rlDAT(MI N) 

A (MFVE)tA (MFVF>*100./51 . 

CALL LI NAP (A (MFVE) > 

A2 CONTINUE 

C SAMPLES OF NOISE 

DO 150 J=1 ,4 
JEL=J+t 

BMEAN(J)=IDAT(JEL) 

BMEAN(J)=RMFAN(J>*t00./51 . 

CALL LI NAP (BMEA N< J) > 

150 CONTINUE 

210 CONTINUE 

p r z 'jp N 

END 


*********************SU3P OUTI NE LI NAP********************* 
LI NAP TRANSFORMS OUTPUT VOLTAGES INTO INPUT VOLTA GFS OF 
THF RECEIVER. THE CALIBRATION DATA IS CONTAINED IN SUB- 
ROUTINE VALUE. 

******************** ******* *** ****** ********** **** ****** ** 
INPUT AND OUTPUT: 

A IS THE OUTPUT VOLTAGE THAT IS TRANSFORMED INTO 
INPUT VOLTAGE 

SUBROUTINE LI NAP (A > 

COMMON / / S(29),TU(29>,TUOC50) 

IECA.GT,1000.0)A=1000.0 
DO 5 I =1 ,29 
J=I + I 

IF(A .GT.TUOCD.AND.A .LE.TUO(J) ) GOTO 10 
5 CONTINUE 

A =0. 

PETURN 

I* A=(A-TUO(I))*S(I)+TU(I) 

I F (A .LT.0.0>A =0.0 

PETURN 

END 
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*********************$UBROUTI N£ CALC********************* 
SUBROUTINE CALC CALCULATES AND PRINTS ELECTRON DENSITY 
PROFILES BETWEEN 60 AND 90KM. 
********************************************************* 

C INPUT DATA* 

C ARRAY(I): AX/AO RATIOS FROM S3 TO 90KM 

C LL=LOWER LIMIT IN THE CALCULATION (1=60.7KM) 

C LH=UPPER LIMIT IN THE CALCULATION <20=89. 2KM) 

C R ESP COLLISION FREQUENCY PROFILE < SUMMER, W I NTER OR EQUINOX) 

C 

SUBROUTINE CALC (ARRA Y,LL,LH ,R ESP) 

DIMENSION ARRAY (21 ),P(21),R (3>,CF(3) ,EL<20) 

DATA SUM, VI N,EQU/5HSUMME,5HVI NTE,5HEQUI N/ 

IF (RESP.EQ.RUM) GO TO 100 

IF (RESP.EQ.WI N) RO TO 200 

IF (RESP.EQ.SQU) RO TO 300 

0 COLLISION FREQUENCY PROFILES 

100 PCD =240* 

P (2)=182. 

P<S>=140. 

P(4) = 108. 

P<5)=83. 

P(R>rfi3. 

P(7>=49. 

P<8)=38. 

P(9>=29. 

P(10)=22. 

P (1 1>=1 7. 

P<12)=12.9 
P C!3) = I0.1 
P (14>=7.7 
P (15)=R. 

P ( ! S)=4.S 
P ( 1 7) =3 .5 
P (18>=2.7 
P(I9) = 2.02 
P (20): 1 .55 
P (2 1 ) : 1 .2 0 
GO TO 400 
2 00 P ( 1 ) = 1 B3 . 

P (2) = 1 45 . 

P (3)=1 1 0. 

P (4): 89, 5 
PC5) = <58. 

P (R):5R. 

P(7) = 44. 

P(P):35. 

P(9) = 27.5 
P ( 1 3 ) : 2 1 ,5 
P(!I):17. 

P ( 1 2 ) : 1 3 . 2. 

P(!3):10,4 
P(14):7.P 
P ( ! 5) :R. 1 
P (10=4.6 
P ( 1 7) :3 .5 
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P( 1R)=2.65 
P(1S>):2,1 
P(20):l .5 
P <2 ! ) = l .2 
GO TO 400 
300 P<!)=220. 

P(2):172. 

P (3) =132. 

P(4):|04. 

P(5> = P0. 

P(6):62. 

P <71 = 48. 

P(R) = 37. 

P(9)=2P.5 
P ( 1 0> = 22 . 

P ( 1 1 ) = I 7, 

P(12)=13. 

P (13) = ! 0. 

P ( 1 4) :7.F 
P(15)=6.0 
P (16)=4.6 
P (17>:3.S 
P (IP):2.75 
P (19>=2.12 
P(20)=l .65 
P (21 ) = ! .35 

400 DO 401 I =1 ,21 
P(I) = PCI)*(10.**5) 

401 CONTINUE 
CALL HEAD (1 ) 

K=0 

DO 20 I=LL,LH 

C CALCULATION OF ELECTRON DENSITIES 

R (1) =ARRAY(I ) 

P <2):APRAY(I+n 
K =K+ 1 

CF(1 )=P (K) 

CF(2):P (K+l ) 

I F(R ( I ) .FO.0.0.OR.R (2).EO.0.0) GO TO 20 
C FUNCTION ELDEN CALCULATES ELECTRON DENSITIES 

FL(K):ELDEN(R,CF, t .5E+3) /< 1 0.**S) 

20 CONTINUE 

C PLOT (J) PLOTS THF ELECTRON DENSITY PROFILE 

CALL PLOTJ (EL) 

RETURN 

END 
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*********************SUgR0U7I NE VALUE* ******************** 
VALUE CONTAINS THE CALIBRATION OF THE RECEIVER AND IS 
CALLED BY LI NAP, TO CORRECT THE PARTIAL REFLECTION DATA 
FROM NONLINEARITY OF THE RECEIVER. 

********** ************41 *********************************** 


S ( I ) =PA TI 0 OF INCREMENTS AT INPUT AND OUTPUT OF RECEIVER 
TU ( I ) = I NPUT OF RECEIVER 
TUO (I ) = OUTPUT OF RECEIVER 


SUBROUTINE VALUE 

COMMON / / S(29),TU(29),TUO(30> 


S (1 ) = 

4.467 

S (2 ) = 

1.722 

S(3) = 

1.247 

S (4> = 

1.115 

S (5> = 

. 988 

S (6) = 

.797 

S (7) = 

.668 

S(8> = 

.639 

S (9>r 

.569 

S (1 71) 

= .493 

sun 

= .451 

S(I2) 

= .427 

S (IS) 

= .414 

S (1 4) 

= .391 

S ( I 5) 

= .398 

S (1 6) 

= .363 

S O 7) 

= .344 

SCI8) 

= .348 

S(1 9) 

= .360 

S (20) 

= .393 

S (2 1 ) 

= .432 

S (22 ) 

= .492 

S (23 ) 

= .511 

S (24) 

= .583 

S (25) 

= .689 

S (26) 

= .676 

S (27) 

= .796 

S (28) 

= .839 

S (29) 

= 1.112 

TU( I ) 

= 0. 

TU(2) 

=4.467 

TU (3 ) 

= 6.310 

TU(4) 

=7.943 

TU(5) 

=8.913 

TU(6) 

= 1 0.0 

TU ( 7) 

=12.589 

TU(*> 

=14.125 

TU(9) 

=15.849 

TU(10>=17.783 

TU ( N 

>=19.953 

TU(12>=22.387 

TU ( 1 3 

>=25. 1 19 

TU ( 1 4 

) =28. 1 84 

TU(15>=3! .623 

TU(16)=35.481 
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TU(!7)=44.66fi 
TU<1?)=70.795 
TU(19) = 1 12.202 
TU (20) = 125. 893 
T U<2 1 ) =1 58.489 
TU (22) =1 99.526 
TU(23)=25! .! B9 
TU (24) =28! .838 
TU(25)=31 6.228 
TU(26)=354.813 
TU(27)=398.107 
TU(28)=446.6B4 
TU(29)=50l .1 87 
TU0(l>=-.4 
TUO (2) = • 6 
TU0(3)=1 .67 
TUO (4) =2 .98 
TUO (5) = 3. 85 
TUO(6)=4.95 
TU0(7)r8.2 
TUO (8)=I 0.5 
TU0(9)rl3.2 
TUO (1 0) = l 6. 6 
TUO (1 1 )=2 1 . 
TU0(12)r26.4 
TU0(13)=32.8 
TUO(14> = 40.2 
TU0(15)=49. 
TUO(l 6)=58,7 
T UO (1 7) =84 * 
TUO< 1 8) = 1 60. 
TUO ( J 9) =279 . 
TUO(20)=317» 
TUO (21 ) =480* 
TUO (22) =495 . 
TUO (23) = 600. 
TUO (24) = 660. 
TUO (25) = 71 9. 
TUO (26) = 775 • 
TUO (27) =839 . 
TUO(2R>=900. 
TUO (29) = 965 . 
TUO (30) = ! 020 . 
RETURN 
FND 
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*********************** y^LS************************* 

VALS CALCULATES ELECTRON DENSITIES AS WOULD BE MEA- 
SURED FROM PARTIAL REFLECTIONS, IN A REGION WHERE 
THE SCATTERING CROSS-SECTION PER UNIT VOLUME CHAN- 
GES IN STEPS AS A FUNCTION OF HEIGHT. 

** * * ** * * * *** * * ** ** * * £ *** ** * 4c + ** * * ** * * * 4c * *** * * * * 


INPUT: 

S! SCATTERING CROSS-SECTION BELOW A HEIGHT ZA 

52 SCATTERING CROSS-SECTION BETWEEN HEIGHTS ZA AND ZB 

53 SCATTERING CROSS-SECTION ABOVE A HEIGHT ZB 
ZA,ZB=HEIGHTS DEFINDED ABOVE 

CE=C OLLISI ON FREQUENCY, ASSUMED INDEPENDENT OF HEIGHT 
EDRrREAL ELECTRON DENSITY OF THE REGION, ASSUMED CONS- 
TANT 

W=PULSE WIDTH OF THE TRANSMITTED PULSE FOR P.R. MEASUREMENT 


1 

1 I 

1 ?. 


13 


14 

1C 

15 
1 7 

1 P 

1 9 


DIMENSION RO ( 1 00) ,R X ( 1 00) ,R A ( 1 00) 

C 52 (X>= (X*(X*(X+G.G9459>+1 G.901 00>+I . 1 G30G) /(X*(X* CX*<X*CX+ 
! G.G31 45 )+35 .3 5520 + SB. 92050) +64 .09340+4. 3605 7) 

C =2. 9 9792 ES 
VRI TE ( G» 1 0) 

FORMAT C3H sn 

READ(4 t ll)Sl 

FORMAT(F10.2) 

VR ITE ( G, l 2) 

FOR MAT (3H S2 ) 

READ (4,1 1)S2 
WRITE<G, 13) 

FORMAT (3H S3) 

READ (4,1!) S3 
WRITE(G,14> 

FORMAT (I3H CF, FORMAT E) 

R FAD (4 , 1 5) CF 
VRITE(G,1 G) 

FORMAT (13H ED, FORMAT E) 

READ(4,15)EDP 
FORMAT (El 0.2) 

WRITE(G,17) 

FOR MAT (3 H ZA) 

READ C4,l DZA 
WRITE(G, 1 P) 

FORMAT (3H ZB) 

READ (4,1 1)ZB 
WR I TE ( G, I 9) 

F OR MA T ( 1 2 H PULSE WIDTH) 

R EAD (4,1 I ) W 
C0=2.59G1 4E7/CF 


CX=7.3*FGEG/CF 

FO=((5./4.)*(3lR2.G0! /C)*C?2(C0))/C0 
FX=((5./4.)*(?1 R2.GF1 /C ) *C52 (CX) ) /CX 
A F0= p 0* EDR 
flFX=FX*EDR 
Z l=ZA-W/4.-2000. 

N= ( ZB -Z A + W/2 .+3 000 . ) *2 . F-3+2 . 

A OP =EXP (AFO* W) 
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AXP=EXP(AFX*V) 

AXN=EXP<-I S *AFX*W) 

Cl=ZA-V/4. 

C2=ZB-V/4. 

C3 -ZA+W/4 • 

C4=ZB+W/4. 

C CALCULATION OF THE REFLECTED SIGNALS, AO AND AX 

WR I TEC 6,20) 

?3 FORMAT (//7X ,2HPA , 1 9X »2HR 0,1 0X,2HRX, l 0X,2HZ 1 /) 

DO 45 1=1, N 
Z1:Z 1+500. 

IFCZI .OT.Ct .AND ,Z1 .LE.C2) GO TO 31 
IFCZI .GT.C2.AND.Z1 .LE.C3)G0 TO 32 
I FCZ I .0T.C3 .A ND .Z 1 *LE ,C4 ) GO TO 33 
IFCZI .GT.C4)Q0 TO 34 

FOCI ) = CSI/C4.*AFO))*(AOP-AON)*EXPC-4.*AFO*Zl> 

RXa) = (Sl/<4.*AFX)>*CAXP-AXN)*EXPC-4.*AFX*Zn 

00 TO 304 

31 RO(I) = CCS2-S) )/C4.*AF0))*EXPC-4.*AF0*ZA)+CCSl*A0P-S2*A0N>/ 

! C4.*AF0) )* EXP C-4 «*AF0*Z1 ) 

PXCI)=C(S2-Sl )/C4.*AFX))*EXP(-4.*AFX*ZA) + CCSI*AXP-S2*AXN)/ 

1 C4.*AFX))*EXPC-4.*AFX*Z1 ) 

00 TO 304 

3? ROCI) = CCS2-S1 ) /<4.*AF0) )* EXP C-4 ,*AFO*ZA)+ ( CS3-S2 ) / 

1 C4.*AF0))*EXPC-4.*AF0*ZB) + < (SI *A 0P-S3*A ON) /( 4.*AF0> )* EXP (-4 ,* 
2A FO*Z l ) 

RXCI>=CCS2-S1) /<4.*AFX))*EXPC-4.*AFX*ZA> + ( CS3-S2) / 
2C4.*AFX))*EXP C-4.*AEX*ZB)+( CSl *AXP -S3*AXN) /( 4 .*AFX) )*EXP (-4 .* 
2AFX*Z 1 ) 

00 TO 304 

33 R0CI> = C(S3-S2)/(4.*AF0))*EXP(-4.*AF0*ZB> + ( (S2*A0P-S3* 

1A ON)/(4.*AFO) )* EXP C-4.*AFCM<ZI) 

RXCI) = CCS3-S2)/C4.*AFX))*EXPC-4.*AFX*ZB) + (CS2*AXP-S3* 

2AXN) /C4.*AFX))*EXPC-4.*AFX*Zl > 

00 TO 304 

34 RO(I)r (S3/C4. *AFO) )* <A OP -A ON)* EXP (-4 .*AFO*Zl ) 

P X ( I ) = < S3 / ( 4 . *A FX ) ) * (A XP -A X N) * EXP C -4 ,*A FX*Z 1 > 

3 914 IFCROCI).EP.0.)GO TO 305 

RA CI)=SORTCRXCl) /ROCI)) 

00 TO 306 

305 PA (I ) = I 300. 

306 VRITE(6,307)PA (I) ,PO(I) ,PX<I) ,Z t 

3 0, 7 FORMAT (4F 12.3) 

45 CONTINUE 

C CALCULATION OF FLECTRON DENSITIES AS OBSERVED FROM 

C PARTIAL REFLECTIONS 

ZZ =ZA -W/4.-1 250. 

WR ITF ( 6 , 47) 

47 FORMAT (//7X,?.HED, 10X.2HZM/) 

J - M- 3 

DO 46 K=1 ,J 
L =K+3 
C CL =K 

ZM=ZZ+50O.*CCL 

FD=ALOO(RA (X) /RA (L) ) / (3O00.* (FX-FO)) 

WR I TE ( 6,3 I O) ED ,Z M 
31 0 F OR MA T ( 2 El 2.3) 

46 CONTINUE 
0 0 TO 1 
STOP 
END 
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******************** \jftl£X************ ************** 

VALEX CALCULATES ELECTRON DENSITIES AS WOULD BE MEA- 
SURED FROM PARTIAL REFLECTIONS, IN A REGION WHERE 
THE SCATTER I NO CROSS-SECTION PER UNIT VOLUME CHANCES 
EXPONENTIALLY 

***************************************************** 

I NPUT: 

SIrSCATTERI NO CROSS-SECTION BELOW A HEIGHT ZA 
S3 =S CA TTERI NO CROSS-SECTION ABOVE A HEIGHT ZB 
B =EXPOE NT OF THE SCATTERING CROSS-SECTI ON ,S :EXP <RZ) , 
BETEWEEN HEIGHTS ZA AND ZB 
ZA,ZB=HEIOHTS DEFINED ABOVE 

CF=C OLLI SI ON FREQUENCY, ASSUMED INDEPENDENT OF HEIGHT 
EDRrREAL ELECTRON DENSITY OF THF REGION, ASSUMED 
CONSTANT 

WrPULSE WIDTH OF THE TRANSMITTED PULSE FOP P.R. MEASUREMENT 
DIMENSION RO(!fl*),PA (1PIO) ,PX(1<W 

C52(X): (X*<X*(X+6.69459) + 16.9?!lP>e» + l .1 6306) /(X*(X*<X*(X*<X+ 
ll6.63l45)+35.3552 6)+6R.92050>+S4.09346)+4.3 6P'57) 

C ~ 7 . .99792 EP 

f WR I TE < 6, 8) 

R FORMAT (3 H sn 

READ (4, l 1 >S 1 
WRITE <6,1 0) 

1 8 FORMAT C24H B (S=EXPCRZ>), FORMAT E) 

READ<4,Ii)B 

11 FOR MAT (FI PI .2) 

WRITE (6,1 6) 

16 FORMAT <2 OH S3, IF NEG=EXP <RZR) ) 

READ (4,11) S3 
WRI TE ( 6, 1 7) 

17 FORMAT (3H ZA ) 

R EAD (4,11 )ZA 
WRITE(6,19) 

19 FOR MAT (3 H ZP ) 

R EA D ( 6 , 1 1 ) ZP 
WR I TE (6 , 1 2) 

12 FORMAT (5H CF-E) 
n EAD (4,20 CF 

2™ FOR MA T (El 0.2) 

WR I TE ( 6 , 13) 

13 F OR MA T ( 1 2 H PULSF WIDTH) 

R EAD (4 , 1 I ) W 
WR I TE ( 6 , I 4 ) 

14 FORMAT (5H ED-F) 

P EAD (4 ,2R) EDP 
IF(S3.LT.0.)S3=EXP(B*ZP) 

C CALCULATION OF REFLECTED SIGNALS, AO AMD AX 

C 0:2.5961 4E7/CF 
CX-7.3RR6F6/CF 

F0:((5./4.)*(31R2.6PI1/C)*C52(C0))/C0 
FX:((5./4.)*(31R2.6P)!/C)*C52(CX))/CX 
AFOrFO* EDR 
AFX:EX*EDR 
n 0:R -4.*AF0 
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BX=B-4.*A FX 
C l=ZA-W/4. 

C2=ZB-W/4. 

C3=ZA+W/4. 

C4-ZB+W/4. 

Z !=ZA-V/4.-20fl0. 

N= (ZB-ZA+3000.+W/2.)*2 e E-3+2 sl 
WRlTE(fi,2l) 

2.1 FORMAT(//7X,2HRA,10X,2HRO # 10X # 2HRX f ! 0X f 2HZl/) 

DO 45 I =! ,N 
Zl=ZH-50P. 

IF(ZI.GT.C1 .AMD.Z1 ,LE.C2>00 TO 51 
I FCZ ! .OT.C2.A ND .Z 1 .LF.C3) 00 TO 52 
IF(Z1 .0T.C3.AND.Z! ,LE.C4)00 TO 53 
I F (Z 1 .AT .04)00 TO 54 

PO(I) = (SI /(4,*AF0))*(EXPCAF0*V)-FXP(-l ,*AFO* W>>* EXP ( -4 .*AF0*7 I ) 
RX<I)=(S1 /(4.*AFX))*(EXP (A FX*W) -FXP C -1 ,*A FX*V) )* EXP (-4 ,*A Fy*Z 1 ) 
no TO 42 

51 P 0(1 )= (SI /( 4 »*AFO) ) * (EXP (-4 ,*AFO* (Z 1 -W/4 . ) ) -FXP (-4.*AF0*ZA)) 
l + (l ./RO)*(FXP(RO*(Zl+V/4.)>-EXP(BO*ZA)) 

RX(I)=(S1 /(4.*AFX))*(EXP (-4,*AFX*(ZI -W/4 . ) ) -EXP (-4 .*AFX*ZA ) ) 
1+(1 ,/BX)*(EXP(BX*(Zl+W/4.))-EXP(RX*ZA)) 

00 TO 42 

52 R 0(1 ) = (SI /(4.*AFO))*(EXP(-4.*AFO*(ZI-W/4.))-EXP(-4.*AFO*ZA>) 

1 + (S3/(4.*AFO))*(EXP(-4.*AFO*ZB)-EXP(-4.*AFO*(Zl+W/4.)) > 

2+ ( 1 . /P 0) * (EXP (BO*ZR) -EXP (B 0*ZA ) ) 

RX(I): (SI /(4 ,*AFX) )* (EXP (-4.*A FX* (Z I -W/4 . ) ) -EXP ( -4 ,*AFX*ZA ) ) 
l + (S3/(4.*A FX) )* (EXP (-4«*A FX*ZB> -EXP (-4 .*A FX* (Zl +W/4 ,) ) ) 

2+ ( 1 ./B X)* (EXP (BX*ZB ) -EXP (BX*ZA ) ) 

00 TO 42 

53 R0(I) = (1 ,/BO)*(EXP(BO*ZB)-EXP(BO* (Z 1 -W/4. ) > ) + (S3/( 4 e *A FO) )* ( 

! FXP(-4.*4FO*ZB)-EXP(-4.*AFO+ (Zl+V/4.))) 

PX(I)=(1 . /BX)* (EXP (RX*ZB> -EXP (BX*(Zl -W/4, )) )+(S3/( 4. *AFX>)*( 

1 EXP (-4 ,*AFX*ZB ) -EXP (-4 ,*A FX* (Z l+W/4 „ ) ) ) 

00 TO 42 

5 4 FOCI ) = (S3 /( 4. *AFO))* (FXP (A FO*W) -FXP (-1 .*4 F0*W>)* EXP (-4. *AFO*Zl ) 

RX(I) = (S3/(4.*AFX))*(EXP (AFX*W) -EXP ( - 1 .*AFX*W) )*EXP (-4.*AFX*Z 1 ) 

42 I F(RO(I> «E0 .0.) 00 TO 44 
R A (I ) = SORT (RX (I)/RO(I)) 

00 TO 43 

44 PA(I)-1000. 

43 WRITE (A, 41 )FA (I ) p P 0(1 ) ,R X(I ) ,Z 1 

41 P0PMATC4E12.3) 

45 CONTINUE 

C CALCULATION OF ELECTRON DENSITIES AS WOULD RE 

C MEASURED FPOM PARTIAL REFLECTIONS 

ZZ=ZA -W/4. -1250. 

WR I TE (C, 4 7) 

47 FORMAT (//7X,2HED,I0X,2HZM/) 

J-N-3 

DO 4S I =1 ,J 
L = I +3 
CCL = I 

7M=ZZ+500.*CCL 

FD=ALOO (RA (I) /RA (L ) ) /( (FX-FO) *3000.) 

VRITF(C,310)ED*?M 
312 FORMAT (2 El 2. 3) 

4<? CONTINUE 

00 TO 1 
STOP 
END 
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C ************************ADI st************** **** ********** 

c adist prints a histogram of the distribution of values 

C OF AO, AX, OR AX/AO, 

c ********************************************************* 

c 

DIMENSION FNAM(2),A 0(21 ) ,AX(21 ) ,BN0(4> ,BNX(4) ,1 FA (3 I ,21 ), 
2IP (2 1 ) ,AXA0(2 ! ) ,A OE (21) ,AX2(2I) 

PI CALL HEAD (0) 

IRNO=0 
I R NX =0 
CALL VALUE 
CALL DINIT 
VRITE(6,20> 

20 FORMAT ( 1 5H WHICH DATAFILE) 

READ (4,30) FNAM 
30 FORMAT (2A 5) 

CALL FSTATC2, FNAM.LOG) 

I F(LOG . NE. 0) GO TO 40 
VR I TE ( A , 35 ) FNA M 

35 FORMAT (AH FILE ,2A5,!9H NOT FOUND ON DAT 2) 

GO TO 10 

40 CALL SEEK (2,FNAM) 

51 FORMAT (F 1 0.0) 

C CV IS THE CENTRAL VALUE OF THE HISTOGRAM 

VRITE(A,52) 

52 FORMAT ( 1 4H CENTRAL VALUE) 

R EAD (4,51 )CV 

C STEP IS THE INTERVAL IN VALUES OF AO, AX OR AX/AO 

C TO RE USED IN THE HISTOGRAM 

WRITE(S,54) 

54 EORMATdlH WHICH STEP) 

R EAD (4 ,5 1 )STEP 

C SAMPLES WILL BE REJECTED IF NOISE IN THE ORDINARY 

C OR EXTRAORDINARY FRAME EXCEEDS A VALUE BMXNS 

WRITE (6,5 6) 

5 6 FORMATd 0H MAX NOISE) 

READ (4,51 )BMXNS 

C SAMPLES WILL BE REJECTED IF THE SIGNAL TO NOISE PATIO 

C IS BELOW A VALUE SNR 

WRITE (6,57) 

E OR MA T (I 9H SIGNAL-NOISE RATIO) 

PEAD(4,51 )SNR 

C SAMPLES WILL BE PFJEOTED IF THEY INCREASE BY MORE 

C THAN A UALUE RE, IN RELATION TO THE SAMPLE TAKEN 

C AT THE SAME HEIGHT, I MMEDI ATELLY BEFORE 

WRITE(6,5P) 

5F FORMATd 1H MA X FADING) 

R FAD (4 ,5 1 )PE 
WRITE (6, 62) 

62 FORMAT (24H AO, AX OR AX/AO (1,2,3)) 

R EAD (4, 63 >1 S 

63 FORMATd 1) 

KEOFO=0 
KEOFX=0 
I D =0 

DO 64 I N: 1 ,2 1 
IR (I N):0 
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A 0 < I N> 

R4 AXdN>r0, 

DO €5 I =1 ,31 
DO «5 Js 1,21 
*5 IFA(I,J)2f 

CVS=CV-STEP*1S. 

CMI MsCV-1 5 **STEP 
CMAX=CV+15.*STFP 
4P DO EG I =1 ,21 

A02(I):A0(I) 

AX2(I)=AXa> 

CALL DREAD (A 0,B NO t I ERR , ID,KEOFO) 

I FCKEOFO.EO.I ) GO TO 5 PI 

CALL DREAD (AX,B NX,I ERR »I D ,KEOFX> 

I F (KEOFX *EO .1)00 TO 50 
B"!EANO = 0. 

PMEANX=PI. 

DO El I =1 ,4 

BMEAN0=BMEAN0+BN0d)**2 
*1 BMEANX=BMEANX+BNX(I>**2 

pnIEA NOrSQRT (BMEA NO/4 . ) 

P ME A NX ^OR T (BMEA NX /A . ) 

IF(IS*EQ.l>BMFANX=fl. 

IF(IS.EO.2)BMEAN0:P). 

PLO:S NR*B*1EA NO 
BLX=SNR*BMEANX 
DO ^ I N= 1 ,21 
I F (IS .EO .1 )AXd N)t5 00)» 

IFCIS.EO .2)A0<IN)=5P!P!. 

I FCPMEA NO .GT .BD1XNS .OR ,B MEA NX .GT «B MXNS .OR , A 0 d N) .LT.BLO.OR. 
1AX(IN) .LT.BLX) GO TO FI 
D I FE=A 0 (I N) -A 02 <1 N) 

lF(DIFE.GT.PE.OR.AO(lN).GT.510..OR.AXdN> .GT.5I0.)GO TO FI 
79 DIFE=AXCIN)-AX2(IN> 

IF(DIEE.GT .RE) GO TO FI 
GO TO F2 

«1 IRdN)=IRdN) + l 

GO TO 95 

92 A XA 0 (I N) =A X <1 N) /A 0 (1 N) 

I F(IS «EQ »1 )AXA 0 (I N) =A 0(1 N) 

IF(IS.EQ.2)AXA0(I N)rAX(I N) 

IFfAXAO(lN) .LE.CMI N) GO TO 9! 

IFCAXAOCIN) .GF.CMAX) GO TO 92 
DO P4 I ,31 
Cl =1 

C =CVS+CI*STFP 
Ct*=C-STEP/2. 
pp STEP /2 . 

1F(AXA0(IN) !gF.CM*AND.AXAOCIN).LT.CP)GO TO 90 
F4 CONTINUE 

9F IFA (1 ,1 N) = IFA (I ,1 N) + l 

GO TO 95 

91 IFA(1,IN):IEA(1,IN) + ! 

GO TO 95 

92 I FA (31 ,1 N)=I FA (31 ,1 N) + l 

95 CONTINUE 

GO TO 4F 
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5 fl I D =1 D /2 

HT=60. 

CALL HEAD (I ) 

DO 120 1:1,21 

WRITE (6, 1 1 fl)HT,ID,IR (I ) 

110 FORMAT (///2X, GHHEIGHT,F4, 1 /2X»I4, !X,13HSAMPLES TAKEN, 3X,I4, 
21 X, FHREJECTED /) 

CS=CMI N 
WRITER, 1 12) 

I 12 FORMA T <4X , 5HVALUF ,3X» 4HFP EQ , 4X,5HVA LUE,3X, 4HFREO » AX , 

25H VA LUE, 3X, 4HFREO, 4X, 5HVALUE ,3X, 4HFREQ /) 

DO 1 IP J:1 ,29,4 
CS1=CS+STEP 
CS2:CSl + STFP 
CS3 :CS2 + STEP 
I F (J*EQ ,29) 00 TO 113 

WRITE (S, 1 14)CS,IFA<J,I ),CS1 ,IFA(J+1 ,1 ) ,CS2 ,1 FA (J+2 ,1) 

2, CS3, 1 FA (J+3, 1) 

00 TO HR 

113 WF I TF (6, 1 14)CS,IFA (J, I ) ,CS1 ,1 FA (.1+ I ,1 ) ,CS2 ,1 FA (J+2 ,1 ) 

114 F0PMAT(4(F9.2,I7)) 

HR CS=CS+STEP*4. 

120 KT=HT+1.5 

00 TO 10 
STOP 
END 


ORIGINAL! PAGE IS 
OF POOR $UALOT 
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C ak***********************ADI STT**************************** 

C ADISTT PRINTS A HISTOGRAM OF THE DISTRIBUTION OF VALUES 

C OF AO, AX, OR AX/AO, FOR DATA TAKEN WITH A PROGRAMED ATTE- 

C NUATOR, THAT OPERATES ON ALTERNATE FRAMES, 

C ********************************************************* 

c 

DIMENSION FNAM(2> ,A0(21 ) ,AX(2I ) ,B NO (4) ,B NX (4) ,IFA (31 ,21) , 
2IRC2I) ,AXA0(2l),A0T(2l),AXTC21 ) ,B NOT (4) ,B NXT ( 4 ) ,A 02 (21 ) , 
3A X2 (2 l > 

IP! CALL HEAD (0) 

IRNOrfl 
IRNX=0 
CORR -! . 

DO 15 1=1 ,31 
DO 15 J = 1 ,21 
15 IFA<I,J) = fl 

CALL VALUE 
CALL DINIT 
WRITE(5,2P!) 

2P! FORMAT (15H WHICH DATAFILE) 

READ (4,30) FNAM 
30 FORMAT (2A5) 

CALL FSTAT (2, FNAM, LOG) 

I F (LOG . NE . 0) GO TO 4P 
WR I TE ( 5 , 35 ) FNA M 

35 FORMAT (5H FILE ,2A5,19H NOT FOUND ON DAT 2) 

GO TO 10 

40 CALL SEEK<2 ,F NA M) 

51 FORMAT(F 1 0.0) 

C CV IS THE CENTRAL VALUE OF THE HISTOGRAM 

WRITER, 52) 

52 FORMAT (14H CENTRAL VALUE) 

R EAD (4,5 l ) CV 

C STEP IS THE INTERVAL IN VALUES OE AO, AX OR AX/AO 

C TO BE USED IN THE HISTOGRAM 

WRI TE ( 5,54 ) 

54 FORMAT (1 1H WHICH STEP) 

P EAD (4 ,5 1 ) STEP 

C SAMPLES WILL BE REJECTED IF NOISE IN THE ORDINARY 

c OR extraordinary frame fxceeds a value BKXNS 

WRITE (5, 55) 

5 5 FORMAT (1 OH MAX NOISF) 

REA0(4,5t >9MXNS 

C SAMPLES WILL BE REJECTED IF THE SIGNAL TO NOISE RATIO 

C IS BELOW A. VALUE SNR 

WPI TE (5,57) 

57 FORMAT (1 9H SIGNAL-NOISE PATIO) 

READ (4,5 l)S NR 

C SAMPLES WILL BE RFJECTFD IE THEY INCREASE BY MOPE 

C THAN A VALUE PE, IN RELATION TO THE SAMPLE TAKEN 

C AT THE SAMK HEIGHT, IMMFDI ATELLY BEFORE 

WRITE (5 t 5 F) 

5 P FORMAT ( 1 1 H MAX FADING) 

READ (4,51 )PE 
VR I TE ( 5 , 4 1 ) 

41 FORMAT (23 H HEIGHT STARTING ATT AO) 

REA.D(4,5I )AHO 
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42 


4 3 


42 


43 


21 * 
44 


48 

44 


4 ! 


WF I TE (4 42 ) 

FORMAT (23H HEIGHT STARTING ATT AX) 

R EAD (4 , 5 l )A HX 
WRITE (4,43) 

FORMAT (22 H CORRECTION ATT SIGNAL) 

R EAD (4 ,5 1 )CORR 
WR I TE ( 4 » 42 ) 

F OR MA T ( 2 4H AO, AX OR AX/AO (1,2,3)) 
READ (4, 43) IS 
FORMAT (1 1 ) 

IHO: (AHO-58.5) /I .3 

I HX = (AHX-5 8.5) /I .5 

KEOFO=0 

KEOFX=0 

ID =0 

IDM=4 

CALL DRFAD(AO,RNO,IERR,ID,KEOFO) 
CALL DREAD(AX,BNX,IERR, ID, KEOFX) 
CALL DREAD ( A OT,P NOT , I ERR , I D ,KEOFO> 
CALL DREAD (A XT, 8 NXT,I ERR ,ID,KFOFX) 
IFCAOT (9) ,LT.A0(9))G0 TO 210 
IDM = 4 

CALL DREAD (AOT,PNOT,IERR,ID,KEOFO) 
CALL DREAD (AXT,P NXT , I ERR , I D ,KEOFX) 
DO 44 I N=t ,21 
I R ( I N) = 0 
A 0 (I N) : 0. 

AX (I N> = 0. 

CVS^CV-STEP*1 4. 

C Ml N-C V- 1 5 .*STEP 
CMAX=CV+15.*STEP 
DO 44 I = 1 ,21 
A 02 ( I ) =A 0 ( I ) 

AX2CI)=AX(I) 

CALL DREAD (A 0,R NO ,1 ERR ,ID,KEOFO) 

I F (KEOFO.EQ , 1 ) GO TO 50 

CALL DRFA D (A X,P NX, I ERR ,ID, KEOFX) 

I F (KEOFX »EQ . ! ) GO TO 50 

CALL DRFAD(AOT,B NOT , I ERR ,ID,KEOFO) 

I F (KEOFO.EO .1 ) GO TD 50 

CALL DREAD (AXT,P NXT , I FRR , I D , KEOFX) 

IF(KEOEX .EQ .1 ) GO TO 50 

PMEA NO = 0. 

BMEANX:0, 

DO 4| I ; 1 ,4 

PMEA NO=BMFA NO+?NO(I )**2 
PMFA NX -PMEA NX+RNX ( I )**2 
PMEA NG=SORT (BMFA NO/4. ) 

B MEA NX =SORT (PMEA NX/4 . ) 

I F(I S .EG . ! ) PMEA NX "0 . 
IF(IS.EO.2)9MFANO = 0. 

P L 0 =S NR *P MF A NO 
PLX=?NR*FMEANX 
DO 95 IN-1,21 

IF(IN.GE.I H0)A0(IN)=A0T (IN)/CORR 
IF (I N.GE.IHX)AX(IN) =AXT (I N) /C ORR 
SAT=51 0./COPR 


ORIGINAL PAGE IS 
OF POOR ^UALETYi 
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79 

Rl 


92 


*4 

9* 

9! 


92 

95 

58 


111* 


1 12 


113 
! 14 
IIP 
12 P 


IF(lSoEQol)AX(IM>-30fl o 
1 FdS.EQ .2)AO(IN>=300 o 
I F(BMEAN0 o GT.BMXNS o 0R.RMEANX 
lAXdN) e LT.BLX)GO TO 81 
D I FE=A 0 (I N) “A 02 (I N) 

I F(DI FEeOT.RE.OR .A0(1 N) .GT.SAT 


,QT 0 8 MX NS o OR D A 0(1 N) oLT e BLO c 


, OR .AXCI N> .GT .SAT) GO TO 1 


OIFE=AX(I N)“AX2(IN) 
IF(DIFE.GT.PE)GO TO SI 
GO TO S2 
IR (I N)=IR(IN)U 


GO TO 95 

AXAO(IN)=AX<IN) /AO(IN) 

I FCIS.EQ .1 )AXAO(I N) = A O(IN) 

I F(IS »FQ »2) AXAO<X N)=AX(IN) 
IF(AXAO(IN) .l-E.CMl N> GO TO 91 
I F(AXAO (1 N) .GE.CMAX) GO TO 92 


DO P4 1=1 ,3! 

Cl =1 

C ^VS+CI*STEP 
CM=C ->STEP/2. 

C P STEP /2 © 

I F(AXA 0 (I N) IgE.CM.A ND .A FA 0 (I N) „LT.CP ) GO TO 98 


CONTINUE 

IFA (1,1 N>rl FA (I »I N)+i 
GO TO 95 

IFA(1*IN) = IFAU*IN)+1 
GO TO 95 

IFA(31,IN) = IFAC3I,IN) + 1 

CONTI NUF 

GO TO 4S 

ID = (ID -I DM) /A 

HT = KPi. 

CALL HEAD ( 1 ) 

DO 120 I =1 ,21 

URITEC6, nPl)HT,ID,IR(I) _ 

FORMAT (2 X» GHHEI GHT , F4 . 1 /2 X,I 4 , 1 X» I 3HSAMPLES TAKEN,3X,I 
2 SHR EJECTED) 

CS=CMI N 


WRITECM 12) 

F0RMAT(4X,5HVALUE»3X,4HFPEQ»4X fl 5HVALUE,3X»4HFREQ»4X » 

25H VALUE » 3X, 4HFREQ , 4X , 5 RVALUE »3X#4HFRE@ /) 

DO IIS J = 1 ,31 ,4 
CSI =CS+STEP 
CS2=CS1 + STEP 
CS3=CS24-STEP 
I F(J„EO *29) GO TO l 13 

WRITE <6, !!4)CS,IFA (J»I ) ,CS1 ,IFA CJ+ 1 ,1 ) ,CS2 ,1 FA (J+2,1) 
2,CS3,IFA(J+3,I> 

GO TO UR 

WRITE IS, 1 !4)CS,I FA (J f I > ,CSI ,IFA (J+l ,1 ) ,CS2 ,1 FA (J+2,I ) 
FORMAT (4 (F9.2,I 7 )) 

CS=CS+FTEP*4. 

HT =HT+1 .5 
GO TO 18 
STOP 
END 


OR . 

1 


4 , IX e 
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**********************rexard******************************* 

RETARD CALCULATES ELECTRON DENSITIES FROM 72 TO 90 
KM, CORRECTI NR THE VALUES OF AX, DUE TO THE DIFFE- 
RENT GROUP VELOCITIES OF THE ORDI NARY AND EXTRAOR- 
DINARY MODES. 

*********************************************************** 


I 

C 


c 


c 


DIMENSION CF( 1 3) ,FD < 12) ,A 0 (13) ,A X( 1 3) 

C3m = (X*(X*(X*<X+2.4653l 1 El 5 + 1. 1 394 1 6E2)+1 .128751 El > + 
22.39R347E-2) /<X* (X* (X* (X* (X* (X+2 . 465 681 El ) + l .20495122) + 
32.89580RE2 )+l . 4921 25E2 )+9. 387737) + l . 8064 1 2 E-2) 

C5CX)= (X* (X*(X+6. 69459) + ! 6 ,901 00) + l . 1 5306) /CX* (X* (X* (X* CX+ 
26.631 45 )+35 .3552 6) + 68* 92 050) + 64. 093 4 6) + 4. 3 605 7) 

CALL HEAD (0) 

COLLISION FREQUENCY MODEL 

CF(!)=27.5E5 

CF C2) -2 1 .5E5 

CFC3)=17.E5 

CF(4) = 13.2E5 

CF(5)=I 0.4E5 


CF(6) = 7.RE5 
CF(7> = 6.1 F5 
CF(R)=4.6E5 
CF(9) = 3.5E5 
CF(10)=2.65E5 
CF<II)=2.1F5 
CFC12) = 1 .5E5 
CF(13> = 1 .2E5 

INITIAL CONDITIONS AND CONSTANTS 

HT0=72.E3 

XORrfl. 


TXP=0. 

C =2.99792528 
TP =6. 2 83 18 
V=2.66E6*TP 
F=2.66E6 
W0=2 .5 961 E7 
Wy=7.3RR6E6 
A4=2.392347F-2 
A3=l .128751 El 
A2 = 1 .1 3941 6E2 
A 1=2.46531 1 FI 
9 6= I .80641 2 E-2 
9 5=9.3 87737 
9 4=1 .4 92I25E2 
83 =2. 895808F2 
82= 1 .204951 E2 
9 1 =2 .465681 El 
D3= 1 .163064 
D2.= 1 .690I00E1 
D 1 = 6, 6945 93 
F5 = 4 . 3605 73 
F4=6.409346E1 
E3 = 6 . R92050E 1 
E2 = 3 .535525E1 


ORIGINAL PAGE IS 
,03! POOR gUALTO 
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E! =6.63! 449 

READ VALUES OF A0,AX AND ELECTRON DENSITIES CALCULATED 
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C WITHOUT CORRECTION FOR THE RETARDATION OF THE X MODE, 

C BETWEEN 72 AND 90 KM. 

WRITE(6»l !) 

I! FORMAT (20H ED, 72.7 TO 89,2 KM) 

READC4,l2)CED«),I=!,t2) 

12 FORMAT (El 0,2) 

WR I TE ( 6 , 1 3 > 

13 FORMAT (15H AO, 72 TO 90 KM) 

READ (4, 1 5) (A 0 (I ) ,1 =1,13) 

15 FORMAKFl 0.2) 

WR I TE ( 6 , ! 4 ) 

14 FORMAT (3H AX) 

READ (4,15) (AX(I) ,1=1 ,13) 

CALL HEAD ( 1 ) 

WRITER, !0> 

10 FORMAT (4X,2HED, 8X,2HHT , RX,3HMU0,7X,3HMUX,7X,2HT0, 

2BX,2HTX,7X,5HAX/A0) 

DO 30 I =1 ,1 1 
HT-HT0+1 .5E3 
CFM= (CF(I)+CFCI+l))/2. 

0=W0/CFM 

XrWX/CFM 

C DERIVATIVES OF THE FUNCTIONS C3/2CX) AND C5/2CX) 

A N0=O (0* CO* (0+A 1 )+A2)+A3)+A4 

D0=0+ (0+ (0* (0* <0* (0+5 1 ) +52 ) +B3 >+B4 ) +B5 ) +3 6 

A NOD =4 .* 0**3+3 ,*A 1* 0**2+2 ,*A2* 0+A3 

D OD =6 .* 0* *5+5 .* 0* *4*B 1 +4 .+32* 0**3+3 . *B3* O* *2+2 . *B 4* 0+35 
ANX=X*(X*CX*(X+Al)+A2)+A3)+A4 
DX=X* CX*CX*(X*(X*(X+31 )+B2)+B3)+B4)+B5)+B6 
A NXD=4.*X**3+3.*AI*X**2+2.*A2*X+A3 

DXP =S.*X**5+5 .*5 1 *X** 4+4 ,*32*X**3+3 .*B3*X**2+2.*B4*X+B5 
A M05 = 0**3+Dl*C**2+D2*0+D3 
D 05= 0**5+El * 0*4+E2*0**3+E3*0**2+E4*0+E5 
A N05D =5 .* 0**4+4 ,* El *C**3+3.*H2*(X<*2+2.*E3*0+E4 
D 05D =5 .* Qfc*4+4 .* El *0* *3+3 .*E2* 0**2+2»* E3*0+E4 
A NX5=X**3+D1*X**2+D2*X+D3 
D X5=X**5+ El *X**4+E2 *X**3+E3 *X**2+E4*X+E5 
A MX 50 =3.*X**2+2.*D1*X+D2 

DX5D =5 ,*X**4+4.*E1*X**3+3.*E2*X**2+2.*E3*X+E4 
C3 OD =TP* (DO* A NOD -A N0*D0D )/(CFM*D0**2) 

C3XD =TP* (DX*A NXD-A NX*DXD ) /(CFM* DX**2) 

C 50D =TP* CD05*A N05D -A N05*D05D ) /(CFM* D05**2) 

C 5XD =TP*CDX5*A NX5D -A NX5*0X5D ) /(CFM* DX5**2) 

C CALCULATION OF RX/RO 

01 =WO/CF(I ) 

X1=VX/CF(I) 

02 =WO/CF (1+ 1 ) 

X2=WX/CF(I+1 ) 

PX1= (X1*C3 (X1))**2+(2.5*C5 (Xl))**2 
R 0 1 : (01 *C3 (01 ))**2+(2.5*C5 (0I))**2 
PX2: (X2*C3 (X2) )**2+(2.5*C5 (X2) )**2 
R02r ( 02 *C3 ( 02 ) )**2+ (2 .5 * 05 ( 02))**2 
ABB = ( 1 .25*31 82 . 6* (C5 (X)-C 5 ( 0) ) ) /(C*CFM) 

DLOOR =A LOG (RX2/R 02 ) /2 .-A LOG (RXl /R 0 1 ) /2 , 
AXA01=AX(I)/A0(l) 

18 VN2=3 182.6* ED (I) 

C GROUP VELOCITIES FOR THE ORDINARY AND EXTRAORDINARY 



227 


C MODES 

REO=-l .* F*WN2*0*C30D /(TP*CEM) -WN2*C3 (0>*(W0+V)/<TP*CFM**2) 
REXr-1 .*F*WN2*X*C3XD/(TP*CFM) -WN2*C3 (X)*(WX+W)/(TP*CFM**2> 

A I MO= (5 */2 . )*WN2*C5 (O) / (TP*CFM) + (3 ,/2 »)*P* WN2*C50D/(TP*CFM) 
AIMX = (5. /2 . )*WN2*C?(X) /(TP*CFM) + (5./2 »)*P* WN2*C3XD/(TP*CFM) 
A N2 CR=1 .-WN2*W0*C3 ( 0) /(W*CFM* *2 ) 

AN2 01 r(5./2.)*WN2*C5C0)/(W*CFM) 

AN2XRrl .-WN2*WX*C3 (X) /<W*CFM**2> 

A N2XI =(5./2.)*WN2*C5(X)/<W*CFM> 

T 002 = A N2 01 /A N2 OR 
TGX2=AN2XI /AN2XR 
C OSOE =SQRT ( 1 ./(I ,+TG02**2 ) ) 

C OSXSrSQRT ( 1 . / (1 ,+TGX2**2 ) ) 

SINOrSORTCd .-C0S02)/2.) 

cosn=soRT((i .+cosoe)/ 2 .) 

SI NX=SORT ( ( 1 .-C 0SX2)/2.) 

C OSX =SORT ( ( I «+C0SX2) /2 , ) 

AMOD OrSORT (A N2 0I**2+A N2 CR**2> 

A MODX rSPPT (A M2 XI **2+A N2 XR**2 > 

A NOR =A MOD 0* COSO 
A NO! rA MOD 0* SI NO 
A NXR rAMODX*COSX 
ANXI =AMODX*SINX 

A MU Or (A NOR* (PF0+2.*F)+AIM0*AN0I ) /(2.*F* <A N0R**2+A N0I**2) ) 
AMUXr<ANXR*CRFX+2.*F)+AIMX*ANXI)/(2.*F* (ANXR**2+A NXI**2>) 
VOrC /AMUO 
VXrC/AMUX 

C CORRECTION OF AX/AO, PY USING A LINEAR INTERPOLATION 

DH = CHT-HT0) 

T1 OrDH/VO 
T I X =DH/VX 
T0=T0R+T1 0 
TX=TXR+T1X 
DT=TX-TO 

DAX=AX(I+2)-AX(I+l> 

AXC=AX(I+! )+DAX*DT/(DH/C) 

AXA02=A XC/AO (1+1) 

c CORRECTED ELECTRON DENSITIES 

D LOGA rA LOG (A XA 02 ) -A LOG (A XA 0 1 ) 

ED I = (DLOGP -DLOOA ) / (ABB*DH*2. ) 

D ED rABS ( (ED ! -ED (I > ) /FD (I > ) 

I F (DED .LE . 3» 1 5) GO TO 210 
I F (DAX.LF.0.) GO TO 120 
ED(I) = (EDt + ED(n>/2. 

WRITE (6,900) EDI ,FD (I) 

9*0 FORMA T (3 H 1,2E12.3> 

IF(ED (I ) .LT.5.F7.0R.ED (I) .GT.I .E10)GO TO 210 
GO TO 19 

12 0 IF(ED1.GT.ED (I))ED (I)rl ,1*ED(I) 

IF (ED 1 .LT.ED (I ) ) ED (I ) r . 9* ED <1 ) 
IF(ED(I).LT.5.E7.OR.ED(I).GT.l.E10)GO TO 213 
WRITE (6, 901) EDI ,ED (I) 

921 FORMAT (3H 2,2E12.3) 

GO TO 19 
210 TOR:TO 

TXP rTX 
HT0=HT 
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WRITE(6»2R) ED (I) »HT»AMUO»A MUX,TO,TX,AXAOe 
2B FOR MA T ( 1 P 7E 1 0 .2 ) 

30 CONTINUE 

00 TO 1 
STOP 
END 



c *********************** DEC ONM************************** 

C DECONM CALCULATES DECONVOLUTION OF AO AND AX, AND 

C ELECTRON DENSITIES, USING A LEAST SQUARE APPROXI MATI ON 

C ***************************************************** 

C 

DOUBLE PRECISION A 

DIMENSION A 0(26), AX (26) ,G (21 ) ,C0(2t ) ,CX(21) ,A (21 ,21) , 
2S0S (2 ! ) , SXS (21) ,A OD (21 ) ,AXD (21 ) ,X0(2l ) 

CALL HEAD(O) 

C G (I ) APE SAMPLES OF THE TRANSMITED POWER TAKEN EVERY 

C 1.5KM 

WRITE (6, 8) 

F FORMAT (1 7H PULSE, 5 SAMPLES) 

R EAD (4,15) (G(I),I=1 ,5) 

DO 10 I =1 ,15 
K =1+5 

10 G (10 = 0. 

V!R I TE ( 6, 12) 

12 FORMA T ( 1 4K AO, 21 VALUES) 

READ (4,15) (AO(I) ,1=1 ,21) 

15 FORMAT (FI 0.2) 

WRITE (6, 17) 

17 FORMAT (3H AX) 

READ (4, 15) (AX(I ) ,1=1 ,21 ) 

WRITE (6, IP) 

IF FORMAT (24H COLLISION FREQ. PROFILE/ 

2 25H SUMMER, WINTER, OR EQUINOX) 

R EAD (4,1 9)RESP 

19 FORMAT (A 5) 

WR I TE ( 6, 20) 

20 FOR MAT ( 1 5H GAMMA, FORMAT E) ' 

READ (4,22) GAM 

22 FOR MA T(E 1 0.2) 

DO 25 1=1 ,5 
K=I+2t 
A 0(K) = 0 . 

25 A X( K) = 0 . 

C CALCULATION OF EQ (A (N,M) +GA M*B ( N, M) )*S < N) =C ( M) 

0 DEFINE A ( N,M) 

1 DO 35 N=1 ,21 
DO 35 M=1 ,21 

35 A ( N , M ) = 0 . 

DO 40 N= 1 , 1 7 
DO 40 M= I ,5 
Kl = Nf M-l 
K = N- M+ 1 
DO 35 L= 1 ,5 
LM=L+M-1 

35 A(N, K1)=A(N, K1) + G(L)*G(LM) 

I F (K .LE .0) GO TO 40 

A (N, K)=A <N,K! ) 

40 CONTINUE 

DO 43 NM = 1 ,4 
N= NM+ 1 7 
MM 4= N-4 

A (N,NM4)=G(t)*G(5) 

NM3= N-3 
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ACM. NM3)=6(1 )*G(4)+G(2)*G(5) 
IF(N.RT.20)A (N,NM3)=G(l)*G(4> 

NM2 = N-2 

A (N, NM2)sG(l >*G(3)+G(2)*G<4)+G(3>*G(5> 
IF(N,EQ.20)A CN,NM2 > = G(1 )*G(3)+G(2)*G(4> 
IFCN.EQ ,21)A (N»NM2)=G<1)*G<3) 

NM 1 - N-l 
LNr22-N 
DO 42 L = 1,LN 
Ll-L+I 

42 A(N, NM1)=A (N, MM I )+G(L)*G(Ll ) 

DO 43 MM=l ,LN 

M = MM+ M- 1 
KL =LN-MM+1 
DO 43 L=l ,KL 
LMtL+MM-1 

43 A(N,M)=A (N,M) + G(L)*G(LM) 

C ADD GAM*B(N,M) 

DO 5 0 N=l,15 

NM3=M-3 

MM2rM-2 

NMlrM-l 
N0=M 
Nt = M+t 
N2=m-2 
N3=M+-3 

A(M,NM3)=GAM*(-1 ,/2.5+A <M,NM3) 

A (M, NM2) :0AM* (13 * /4 . )+A (M, NM2 ) 

A (M,NMI )=GAM*(-I7./2#)+A(M,NMl) 

A (M, N0)=GAM*49./4.+A (M,N0) 

A (M»N1 )=GAM* (-1 7./2.)+A (M,Nl > 

A (M, N2 ) = GA M* 1 3 * /4 »+A (M,N2 ) 

50 A(M,N3) = GAM*(-I ./2 0 )+A(M,N3> 

A ( 1 , 1 )=GAM* 1 */4.+A (l ,1 > 

A ( 1 ,2) = GAM* (-1 »)+A ( 1 ,2) 

A(l,3) = GAM*5./4.+A(! ,3) 

A (I ,4> = GAM*(-I ,/2.)+A (J ,4) 

A (2,1)=GAM* (-l.)+A (2,!) 

A(2,2)=GAM*17./4.+A(2,2) 

A(2,3)=GAW*(-S.)+A(2,3) 

A (2,4> = GAM*13./4.+A (2,4) 

A (2,5) - GA M* (-1 ./2,)+A (2,5) 

A(3,I ) = GAM*5./4.+A (3,1) 

A (3, 2): GAM* <-6.)+A (3,2) 

A (3,3)= GAM* 42 ,/4,+A (3,3) 

A(3,4) = GAM*(«17./2.)+A(3,4) 

A(3,5) = GAM*13./4.+A (3,5) 

A(3,$)=GAM*(-t ,/2.)+A (3,5) 

A (t 9,1 6) =A (19,1 5) + GAM*(-l./2.) 

.4(19,1 7>tA (I 9,!7>+GAM*13./4. 

A(19,IF)=A (l9,l£)+GAM*(-17./2.) 

A (19, 19) =A (19,19) +GA M* 4 5 • /4 . 

A ( 1 9 ,20) =A ( l 9 ,2 0) + GAM* (- 1 5 . /2.) 

A ( 1 9 ,2 1 ) =A ( 1 9 ,2 1 )+GAM*2. 

A (20,1 7) =A (20, I 7) + GAM* (-1./2.) 

A (20, 1 2) rA (2S , 1 8)+GAM* 1 3 . /4 . 
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A <20,1 9) =A (20, 1 9)+GA W* (-15./2.) 

A (20,20)=A (20, 20)+ GAM* (29. /4.) 

A(20,21)=A(20,21) + GAP*(-5./2.) 

A (21 ,!8)=A(21 ,1 8) + GAM*(-l ,/2.) 

A (2 1 , 1 9) =A (21,1 9) + GAM*2. 

A <21 ,20>=A (21 ,20)+GAM*<-5./2.) 

A (2 l ,21 )=A (21 ,2 1 )+GA M 
DO 52 1=1,21 
CO(I> = 0. 

52 CX(I>=0. 

DO 58 N= I ,21 
DO 55 K=1 ,5 
L=K+ N-I 

C 0 ( N) =C 0 ( N)+G(K) *A 0(L)**2 
55 CX(N)rCX(N)+G(K)*AX(L)**2 

58 CONTINUE 

C INVERSION OF EQUATIONS 

CALL MATI NV(A,21 ,DET) 

DO 61 N=!,2l 
S OS ( N) = 0, 

61 SXS ( N) = 0 . 

DO 62 N=1 ,21 
DO 62 M=1 ,21 

SO$(N)=A (N,M)*CO(M)+SOS(N> 

62 SXS(N)=A <N,M)*CX(M)+SXS(N) 

CALL HEAD (1 > 

DO 68 1=1 ,21 
A OD (I ) = 0« 

68 AXD(I) = 0. 

DO 70 N= 1 ,21 
DO 69 M=l,5 
M N= N-M*- 1 

I F (MN.LE. 0) GO TO 70 
A OD (N)=A OD (N) + G(M)*SOS(MN) 

69 AXD(N)=AXD(N)+G(M)*SXS(MN) 

75 CONTINUE 

DO 72 N=1 ,21 

S OS ( N) : (ABS (SOS ( N) ) /SOS ( N) ) *SQRT (ABS (SOS(N) )) 
SXS(N)=(A3S(SXS(N) )/SXS(N) )*S0RT (ABS (SXS( N) ) ) 

A OD (N) = (ABS (A OD (N) ) /A OD (N) ) *SQRT (ABS (A OD ( N) > ) 

72 AXD (N) r (ABS (AXD (N) ) /A XD ( N) )*SPRT (ABS (AXD ( N) ) ) 

DO 74 1=1,19 
N=22-I 

S OS ( N) =S OS ( N-2 ) 

SXS (N) =SXS ( N-2 ) 

74 XO(N)=A3S(SXS(N)/SOS(N)> 

WRI TE (6,75) GAM 

75 F OR MA T ( 6H GAMMA , 1 PE 1 2 .2 ) 

WRITE (6, 80) (SOS (I ) ,1 =1 ,2 1 ) , (SXS (I ) ,1 = 1 ,21 ) , (AOD (I ) ,1 = 1 ,21 ) , 
2 (AXD (I ) ,1 =1 ,21) 

*0 FORMAT (1P7E1 0,2/1 P7E1 0.2/1 P7E1 0.2//) 

WRI TE (6, 80) (XO (I ) ,1 =1,21) 

CALL CALC(XO, 1 ,20,RESP) 

GA M= GA M* 1 000. 

GO TO 1 

STOP 

END 
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C ******************* SUBR OUTI NE MATI NV* ********************* 

C MATI NV INVERTS A SYMMETRIC MATRIX AND CALCULATES ITS DE- 

C TERMlNANTo 

C ********************************************************** 

C 

C USAGE 

C CALL MATI NV (ARRAY* NORDER, DET) 

C 

C DESCRIPTION OF PARAMETERS 

C ARRAY - INPUT MATRIX WHICH IS REPLACED BY ITS INVERSE 

n NORDER - DEGREE OF MATRIX (ORDER OF DETERMINANT) 

C DET - DETERMINANT OF INPUT MATRIX 

C 

C SUBROUTINES AND FUNCTION SUBPROGRAMS REQUIRED 

C NONE 

C 

C MODIFICATIONS FOP FORTRAN II 

C OMIT DOUBLE PRECISION SPECIFICATIONS 

C CHANGE DABS TO ARSF IN STATEMENT 23 

C 

C COMMENTS 

C DIMENSION STATEMENT VALID FOR NORDER UP TO 21 

SUBROUTINE MA TI NV (ARRAY* NORDER, DET) 

DOUBLE PRECISION ARRAY, AMAX, SAVE 
DIMENSION ARRAY (21*21), IK(21), JK(21) 

10 DET = l. 

11 DO 100 K=l , NORDER 

FIND LARGEST ELEMENT ARRAY(I»J> IN REST OF MATRIX 
AMAX = 0. 

21 DO 3 2 I=K, NORDER 
DO 30 J -K , NORDER 

23 IF CDABS(AMAX) - DABS (ARRAY ( I ,J )) ) 24, 24, 30 
2.4 AMAX = ARRA Y ( I * J) 

IK (K) = I 
JK (K) = J 
3 0 CONTINUE 

INTERCHANGE ROWS AND COLUMNS TO PUT AMAX IN ARRAY(K.K) 

31 IF (AMAX) 41, 32, 41 

32 DET = 0. 

GO TO 143 

4! I r IK (K) 

IF (I-K) 21, 51* 43 
43 DO 50 J:1 , NOPDFR 
SAVE = ARRAY (K,J) 

ARRA Y (K, J ) = ARP AY ( I ,J> 

S3 ARRAY (I , J) = -SAVE 
51 J 3 JK ( K) 

IF CJ-K) 21, SI* 53 
53 DO S3 1=1, NORDER 
SAVE 3 APPAY (I ,K> 

ARRAY (I , K) 3 APPAY(I ,J) 

S0 ARRAY (I , J) = -SAVE 
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ACCUMULATE ELEMENTS OF INVERSE MATRIX 

61 DO 70 1st, NORDER 
IF (I -JO 63, 70, 63 
63 ARRAY(I ,K) = -ARRAY(I,K) / AMAX 

70 CONTINUE 

71 DO 80 Ir|, NORDER 
DO 80 J=1 , NORDER 
IF (I-K) 74, 80, 74 

74 IF <J-K) 75, 80, 75 

75 APRAYCI.J) = ARRA Y(I ,J) + ARRAY (I ,K)*ARRAY<K,J) 

80 CONTINUE 

81 DO 90 J=l, NORDER 
IF (J-K) 83, 90, 83 

83 ARRAY(K,J) = ARRAY(K,J) / AMAX 
90 CONTINUE 

ARRA Y(K, K) = ! , / AMAX 

100 DFT = DET * AMAX 

RESTORE OR DERI NO OF MATRIX 

101 DO 130 L-l , NORDER 
K - NORDER - L + 1 
J = I K (K) 

IF (J-K) 111, 111, 105 
105 DO 110 1=1 , NORDER 
SAVE = ARP AY (I, K) 

ARRAY (I , K) = -A RRA Y (I ,J ) 

110 ARRAY (I , J) = SAVE 
1111= JK (K) 

IF (I-K) 130, 130, 113 
113 DO 120 J=l, NORDER 
SAVE = ARRAY(K,J) 

ARRAYS, J) = -ARRAY(I,J) 

120 ARRA Y(I , J) = SAVE 
130 CONTINUE 
140 RETURN 
END 



